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ANALYSIS  OF  A  LIMITING-AMPLITUDE  PROBLEM 
IN  ACOUSTO-ELASTIC  INTERACTIONS,  I 

1.  INTRODUCTION 

We  are  concerned  here  with  the  correct  formulation  and  the  solution  of  a  boundary  and  interface 
problem  in  elliptic  partial  differential  equations  that  models  the  physical  setting  of  certain  steady-state, 
or  time-harmonic,  fluid-elastic  interactions.  Specifically,  we  consider  an  elastic  shell-like  body  immersed 
in  a  homogeneous,  inviscid  fluid  of  unbounded  extent  in  all  directions.  The  entire  continuum  is  assumed 
to  be  driven  by  the  action  of  certain  known,  externally  applied  sources  of  energy  generating  an  acoustic 
disturbance  propagating  in  the  fluid  and  forces  applied  to  the  surface(s)  and  the  bulk  of  the  elastic  body. 
These  sources  sustaining  the  motion  are  supposed  to  begin  at  some  instant  and  thereafter  approach  harmonic 
time-dependence,  with  a  common  angular  frequency,  and  we  assume  that  the  induced  fields  of  response  in 
the  fluid  and  elastic  media  exhibit  the  same  qualitative  sort  of  temporal  behavior  as  the  process  evolves.  We 
set  the  problem  of  determining  the  time-harmonic  parts  of  these  resultant  fields,  paying  particular  attention 
to  the  solution  for  the  field  in  the  fluid  and  its  far-field  pattern.  To  this  end,  we  must  first  ensure  that 
we  formulate  a  mathematical  model  problem  having  as  its  unique  solution  precisely  the  desired  limiting- 
amplitude  fields  and  then  show  how  to  set  up  a  convergent  scheme  for  their  numerical  approximation. 

Thus,  we  are  examining  a  simplified  version  of  the  problem  that  must  be  addressed  whenever  one  wishes 
to  predict  the  sound  field  that  is  both  radiated  and  scattered  by,  say,  a  vibrating  metal  structure  submerged 
in  an  ocean  in  which  there  is  also  propagating  an  acoustic  wave  originating  from  other  sources.  That  is, 
the  setting  described  is  essentially  a  first  approximation  to  that  prevailing  in  a  broad  range  of  technical 
questions  of  fundamental  significance  to  the  Navy,  including  direct  problems  of  design  and  inverse  problems 
of  detection.  It  is  evident  that  a  thorough  investigation  of  the  solution  of  the  simplified  problem  posed  will 
be  invaluable  as  a  guide  in  attacking  the  study  of  more  realistic  models  incorporating  the  effects  of  spatially 
varying  fluid  properties  and  the  presence  of  air-fluid  and  earth-fluid  interfaces. 

While  the  developments  of  the  subsequent  sections  are  conducted  in  a  systematic  manner,  to  avoid  ob¬ 
scuring  the  basic  issues  in  this  introductory  section  on  motivation  and  orientation  we  shall  be  intentionally 
rather  cavalier  vis-a-vis  matters  of  precision  and  complete  speci '  ■  'on  of  hypotheses.  Consider  a  homoge¬ 
neous  and  isotropic  elcistic  medium,  with  Lame  parameters  A  >  (^  ni.d  /i  >  0  and  density  Qa,  occupying  the 
closure  Qg  of  a  bounded  and  (for  simplicity)  connected  open  set  fio  in  R^.  The  restrictions  specified  for 
the  Lame  parameters  are  equivalent  to  the  reasonable  requirements  that  Young’s  modulus  be  positive  and 
Poisson’s  ratio  be  nonnegative  and  less  than  one  half.  We  suppose  for  now  that  the  boundary  Pj,  :=  dQg  of 
fio  is  “smooth,”  and  denote  by  the  normal  field  of  unit  magnitude  on  To  that  is  “directed  into  the  exterior 
of  Qg."  The  complement  R‘^\f2(,,  which  need  not  be  connected,  is  decomposed  into  its  open  connected  com¬ 
ponents  and  written  as  17+  Ufi-,  with  0+  representing  the  single  unbounded  component  and  denoting 
the  union  of  'h*'  (finite  number  of)  bounded  components.  In  case  n_  is  the  empty  set,  or,  as  wc  prefer  to 
say  here  (for  a  reason  that  will  soon  be  apparent),  “absent,”  one  should  make  the  obvious  adjustments  in 
the  reaisoning.  We  write  r_  :=  dCl-  and  r+  ;=  dQ+,  so  that  I j,  ^  r_  i  i  r+.  and  indicate  by  n_  and  n+  the 
restrictions  of  n,,  to  and  r+,  respectively.  A  homogeneous  fluid  is  supposed  to  fill  f2+;  while  this  fluid 
is  modeled  as  inviscid,  it  is  assumed  that  it  may  possess  some  mechanism  for  internal  damping,  so  that  the 
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propagation  of  small  disturbances  in  the  fluid  is  governed  by  the  “damped”  wave  equation.  The  density  of 
the  quiescent  fluid  is  designated  by  ,  its  characteristic  phase  speed  of  small  disturbances  by  c+ ,  and  the 
damping  constant  by  y+,  with  y+  >  0.  Actually,  we  shall  carry  along  the  hypothesis  that  y+  can  take  on 
certain  complex  values  by  supposing  that 

7+  G  {  2  G  C  1  either  z  =  0orRe2>0}, 

although  an  appropriate  physical  interpretation  is  not  clear  when  Imy^  ^  0.  The  set  f2_,  accounting  for 
any  “cavities”  within  the  elastic  medium,  is  taken  as  evacuated  (i.e.,  containing  no  material);  if  some  other 
continuum  were  postulated  to  lie  in  fi- ,  there  would  be  introduced  certain  additional  complexities  that  we 
wish  to  avoid  in  the  present  exposition. 

Prior  to  some  time  <o,  say,  <o  =  0,  we  suppose  that  the  system  rests  in  a  state  of  equilibrium  in  which 
the  pressure  in  the  fluid  has  the  constant  value  po  throughout  and  no  body-force  field  is  applied  to  any 
portion  of  either  medium.  Then,  reckoned  relative  to  the  in  vacuo  position  of  the  elastic  body,  the  static, 
or  quiescent,  elastic  displacement  field  Uq  :  flo  — ^  satisfies 


Al^Uo  =  0, 

in  fto. 

(1.1) 

T“+[Uo]  =  — pon+, 

on  r+. 

(1.2) 

T“-[Uo]  =  0, 

on  r_ 

(1.3) 

(conditions  serving  to  characterize  Uq  only  to  within  an  added  rigid-displacement  field).  Here,  the  operator 
^  is  given  by 

A)(  ^  :=  -/icurlcurl  (A  -f  2/i)graddiv ; 

in  terms  of  the  Cartesian  components  (Uj)  of  a  vector-valued  field  U,  the  Cartesian  components  of  A* 
then  appear  eis 

(Al^V)j=pAUj  +  (X  +  fi)Ui,tj, 

with  A  ;=  divgrad  denoting  the  Laplacian  (acting  on  scalar  fields)  and  a  comma  signifying  partial  differ¬ 
entiation  with  respect  to  the  Cartesian  coordinates  indicated  by  the  following  subscripts.  Here  and  in  the 
sequel,  the  conventions  of  the  index  notation  are  in  force,  with  the  range  of  all  free  and  summed  indices  to 
be  understood  ais  {1,2,3}.  By  introducing  the  (Cartesian  components  of  the)  strain  tensor  corresponding  to 
a  field  U  in  Q,,  according  to 

+  (1.4) 

and,  in  turn,  the  (Cartesian  components  of  the)  stress  tensor  generated  by  U  as 

(T,i[U]  ;=  Ae„[U]f),*  -I-  2/ie;*[U],  (1.5) 

Sjt  being  the  Kronecker  symbol,  the  traction  field  T"”[U]  induced  by  U  on  To  is  defined  as  having  the 
Cartesian  components  given  by 

(T"»[U])^  :=u-^t[Uj|r„n„fc,  (1.6) 

(uoi ,  no2,  lios)  denoting  Cartesian  components  of  iXo.  when  the  traces  o'jit[U]|r„  exist  in  an  appropriate  sense. 
Recall  that  T"°[U]  represents  physically  the  tractions  exerted  on  the  eleistic  material  in  fio  by  agents  lying 
outside  that  set.  In  (1.2)  and  (1.3)  we  have  indicated  by  T"+(Uo]  and  T"- [Uq]  the  restrictions  of  T"‘’[Uo] 
lo  r+  and  to  r_,  respectively. 
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Subsequent  to  the  time  <o  =  0,  we  assume  that  various  externally  impressed  disturbing  agents  begin 
smoothly  to  act  on,  and  in  a  neighborhood  of,  the  elastic  body,  viz.,  tractions  applied  to  both  the  vacuum- 
elastic  boundary  r_  and  the  fluid-elastic  interface  r+,  a  body-force  field  within  flo,  and  an  “incident  field” 
propagating  in  the  fluid  contained  in  Henceforth,  we  shall  consider  all  scalar  fields  as  C-valued  (C 
denoting  the  complex  numbers),  all  vector  fields  as  C^- valued;  the  corresponding  fields  of  physical  significance 
are  then  to  be  obtained  by  taking  the  real  parts  (when  7^.  is  real).  Points  of  R"*  we  write  in  the  form  (x,  *4), 
with  X  =  (xi,a;2ia;3)  G  R^  and  X4  €  R;  the  operators  curl,  grad,  and  div,  when  applied  to  appropriate 
fields  defined  on  a  cylinder  such  as  fio  ^  R  in  R'*,  shall  be  interpreted  as  acting  in  the  “spatial  variables” 
only.  Thus,  we  are  given  vector- valued  fields  To  :  To  x  R  — ►  and  Fo  :  fio  x  R  — ►  C^,  smooth  in  their 
fourth  arguments,  and  satisfying  To(-,t)  =  0  on  To  and  Fo(-,i)  =  0  in  Qo  for  t  <  0.  Consistently,  as  we 
have  done  for  no  and  T"°[Uo],  when  /„  is  some  (scalar-  or  vector-valued)  field  on  To,  we  shall  denote  by 
/_  and  its  restrictions  to  r_  and  F^.,  respectively,  i.e.,  f±  /o|r±,  the  notation  g\F  being  reserved 
throughout  to  indicate  the  restriction  of  the  function  g  to  the  subset  F  of  its  domain.  Thus,  we  write,  t.g., 
T_(-,t)  ;=  To(-,t)|r-  ^•nd  T+(-,t)  ;=  To(',<)lr4..  To  describe  the  fluid  incident  field,  we  suppose  that 
there  is  specified  an  open  subset  of  R^  containing  the  elastic  body  along  with  its  cavities,  i.e.,  containing 
fi_  U  f2o,  and  a  smooth  scalar  field  :  Q'  x  R  — ►  C,  the  velocity  potential  of  a  disturbance  that  would 
propagate  in  if  no  elastic  medium  were  present  and  the  fluid  were  to  fill  all  of  Q.‘ .  Then  is  to  satisfy 
$'(.,<)  =  0  in  whenever  t  <  0  and 


-A^‘ -f  ^$‘,4 -t- 4-^‘,44  =  0,  infi'xR.  (1.7) 

C4.  C4. 

We  are  implying  here  that  any  “sources”  sustaining  the  disturbance  represented  by  <^'  are  situated  outside 
(and  so,  in  particular,  outside  the  elastic  medium).  In  terms  of  this  velocity  potential,  the  incident 
particle-velocity  field  V‘  and  the  incident  acoustic,  or  excess,  pressure  field  P‘  are  defined  in  x  R  by 

V‘  := — ^grad-J*  (1.8) 

e+ 

(again,  with  the  gradient  operator  acting  in  the  spatial  variables  alone),  and 

P‘ :=  $‘,4-1-7+ (1.9) 

we  shall  use  the  corresponding  rules  of  calculation  when  defining  the  velocity  and  pressure  fields  induced  in 
the  fluid  by  any  other  velocity  potential,  as  well. 

In  the  full  space-time  formulation,  for  the  determination  of  the  resultant  motion  of  the  entire  continuum 
in  response  to  these  specified  disturbances,  we  would  then  wish  to  find  a  vector  (elastic-displacement)  field 
U  :  Qo  X  R  — >■  and  a  scalar  (fluid-velocity-potential)  field  $  :  Q4.  x  R  — ►  C  satisfying,  in  an  appropriate 

sense,  the  hyperbolic  partial  differential  equations 

—  A$  +-  —^$,4  -| - 2'^>44  ~  9, 

c\  ’  c\ 

and 

-b  ?oU,44  =  Fo, 

the  initial  conditions 

$(•,0  =  0 

U(.,<)  =  0 


in  Q+  X  R, 

(1.10) 

in  Qo  X  R, 

(1.11) 

in  Q+  1 

>  for  <  <  0, 

(1.12) 

in  Qo  I 
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the  interface  conditions 


d"  '11+  — 

T“+[U]+  ($,4|r+  +7+4^|r+)n+  =  T+  -  ($U|r+  +  t+*‘lr+)n+, 


and  the  boundary  condition 


T“-[U]  =  T_,  onr_xR. 


on  r+  X  R,  (1.13) 


(1.14) 


Here,  (•)lr+  denotes  a  trace,  and  (•),„^  denotes  a  normal  derivative  on  r+  (taken  with  respect  to  the  spatial 
variables  at  a  fixed  time),  while  a  ■  b  :=  at6t  for  elements  a  =  (01,02,03)  and  b  =  (&i, 62,^3)  of  C^.  The 
equalities  in  (1.10)  through  (1.14)  result  upon  first  setting  down  the  conditions  to  be  fulfilled  by  the  state  of 
the  medium  when  the  elastic-displacement  field  is  calculated  relative  to  the  in  vacuo  position  of  the  elastic 
body  and  then  employing  (1.1)  through  (1.3)  to  remove  the  appearance  of  the  state  of  prestress  induced 
by  the  constant  pressure  po.  Thus,  U  represents  the  perturbation-displacement  field  of  the  elastic  medium, 
taken  relative  to  its  quiescent  position  of  equilibrium  at  times  preceding  to,  the  displacements  relative  to 
its  in  vacuo  position  being  given  by  U  -(-  Uq,  with  Uo  as  described  above.  The  velocity  potential  $  has  no 
intrinsic  physical  interpretation  when  is  nonzero.  Rather,  it  is  the  sum  in  (Q+  n  fi‘)  x  R  that 

describes  the  disturbance  propagating  in  the  fluid  (the  state  of  the  fluid  being  obtained  by  superposition  of 
this  perturbation  onto  the  quiescent  state  with  constant  pressure  po);  when  there  are  neither  tractions  To 
impressed  upon  To  nor  a  body  force  Fo  acting  in  Qoi  it  is  customary  to  refer  to  $  as  the  velocity  potential 
of  the  “scattered  field”  corresponding  to  the  incident-field  velocity  potential  <I>‘.  The  fulfilment  of  the  first 
of  conditions  (1.13)  ensures  the  continuity  of  the  normal  component  of  the  continuum  velocity  on  r+  (which 
is  the  most  that  can  be  demanded,  in  view  of  the  inviscid  nature  of  the  fluid),  while  the  satisfaction  of  the 
second  condition  in  (1.13)  provides  for  the  continuity  of  the  continuum  tractions  on  r+.  Of  course,  with 
a  X  b  denoting  the  vector  product  of  a  and  b  in  C^,  from  (113)  we  have 

n+  X  T'*+[U]  =  n+  X  T+,  on  r+  x  R, 


so  that  the  tangential  component  of  T"+  [U]  on  r+  vanishes  if  T+  is  a  normal  field  on  r+;  this  result  is  also 
due  to  the  inviscidity  of  the  fluid. 


Now,  to  pose  the  problem  of  present  interest,  let  us  consider  the  special  circumstance  in  which  the 
impressed  traction  and  body-force  fields  asymptotically  approach  harmonic  time-dependence,  of  angular 
frequency  w  >  0,  on  the  respective  spatial  parts  of  their  domains  of  definition  as  t  00,  while  the  fluid- 
incident-field  velocity  potential  displays  such  behavior  in  a  neighborhood  of  n_  Uflo-  In  such  a  case,  we  can 
write 


To(',<)  =  T:(.,f)  +  to(-)e--‘ 

on  r„  for  <  6  R, 

(1.15) 

Fo(',<)  =  F:(',<)-hfo()e--‘ 

in  fio  for  <  6  R, 

(1.16) 

$‘(.,t)-<I>‘^(-,0  +  ^‘(-)e-''‘"‘ 

in  for  <  6  R, 

(1.17) 

wherein  fi-Uflo  C  Cl‘  C  fl',  and  the  functions  F^(',<),  and  <^‘’■(■,0,  along  with  sufficiently  many  of 

their  derivatives,  vanish  on  their  respective  domains  in  the  limit  as  /  — ►  oo,  while  to,  fo,  and  ip'  are  certain 
known  (spatially  dependent)  complex  amplitudes.  Again,  we  write  T+(',t)  '•=  TJ(-,t)lr±  and  t±  ;=  to|r±. 
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Then,  from  (1.7),  under  appropriate  conditions  of  transience  of  this  sort  for  and  of  smoothness  for 
the  latter  function  must  satisfy 


or 


Avc>‘  +  k\ip^  =0  in  f2‘ 


(1.18) 


in  which 


K+  =  K+(w;c+,')+)  := 


(w  +  i7+) 


(1.19)- 


Here  and  throughout,  (•)^^^  =  v/T)  denotes  the  principal  branch  of  the  square-root  function:  when  z  is  real 
and  nonnegative,  is  the  nonnegative  square  root  of  z,  and  then,  in  general,  2^/^  ;=  exp(t  arg  2/2), 
with  arg  2  €  (— ir,  jr]  indicating  the  principal  argument  of  z  e  C.  Thus,  /c+(w;c+,0)  =  w/c^..  If  Re 7+  >  0, 
then  the  principal  argument  of  uj{ui  4-  lies  in  (0,  ir),  since 


w(w  -t- 17+)  =  u{lj  —  Im74.)  -I-  iwRe74., 


and  so  its  principal  square  root  k;+(w, c+,7+)  has  positive  real  and  imaginary  parts.  Explicitly,  setting 

R+  =  R+(a>,7+)  :=  if  (w  -  Im7+)  7^  0, 


we  find  that 


and 


K^(u;;  c+ ,  7+)  =  {{1  +  [1  + 


+ 


{l  +  [H-R5.]l/2}l/2i 

if  (w  -  Im7+)  >  0, 


1  /w(Im7+  -w)  / 

'"+(‘^,^+,7+)  -  c+ V  2  \ 


/  ^  (1-1-i)  /a) Re  7+ 

«,(u,^c+,r+)  =  — ^-5- 


(1.20)i 


\  {1  4-  [1  -I-  R2  ]i/2}i/2  +*{!  +  [!  + 

if  (ui  -  Im7+)  <  0,  (1.20)2 

if  (w  -  Im7^.)  =  0.  (1.20)3 


In  view  of  (1.8)  and  (1.9),  V‘  and  P‘  then  have  in  x  R  the  same  form  of  a  time-transient  superimposed 
upon  a  time-harmonic  contribution,  the  complex  amplitude  of  the  latter  deriving  from  for  every  t,  in  fl‘ 
we  have 

V‘(.,0  =  grad<I>*^(-,<)  -  — grad¥>‘(-)e-‘'"‘ 

e+  e+ 


and 


P‘(.,  <)  =  <!>*,;(.,  <)  4-  7+ <!>'"(■.  t)  4-  (7.  -  -‘(Oe—'. 
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Pursuant  to  these  hypotheses,  and  supposing  that  there  exists  a  corresponding  unique  pair  (U,<&)  satisfying 
(1.10)  through  (1.14),  we  further  make  the  fundamental  assumption  that  the  latter  functions  have  the  forms 


U(-,<)  =  U"(.,0  +  u()e--‘ 

in  Qo  for  <  e  R 

(1.21) 

$(•,<)  =  $"(•,  t)-|-^(.)c-“" 

in  for  t  6  R, 

(1.22) 

wherein  and  are  “sufficiently  weak  transients,”  i.e.,  vanish,  along  with  certain  of  their  derivatives, 
in  a  sufficiently  strong  sense  as  t  — ►  oo.  In  this  setting,  our  objective  is  the  determination  of  the  eissumed- 
to-exist  “complex  limiting  amplitudes”  u  in  fio  and  (p  in  in  particular,  we  wish  first  to  formulate  a 
boundary  and  interface  problem  for  elliptic  partial  differential  equations  in  Clg  U  that  completely  and 
directly  characterizes  the  latter  spatially  dependent  fields,  to  determine  when  this  problem  is  well-posed  in  a 
reasonable  sense,  and  subsequently  to  show  how  one  can  construct  sequences  of  approximations  converging 
in  a  useful  manner  to  the  solution  (when  the  problem  is  well-posed).  As  a  preview,  we  note  at  this  point 
that  the  desire  to  account  for  the  effect  of  driving  fields  such  as  To  and  Fo  introduces  into  the  correct 
formulation  of  the  problem  governing  the  steady-state  behavior  peculiar  features  that  are  not  of  concern 
when  it  is  supposed  that  only  a  fluid-incident-field  velocity  potential  is  present  to  sustain  the  motion. 
For  example,  one  would  hope  th^t  the  data  of  the  correct  steady-state  problem  comprise  merely  the  limiting 
amplitudes  to,  fo,  and  along  with  the  angular  frequency  u>,  the  underlying  geometry,  and  the  material 
(constitutive)  parameters,  i.e.,  that  the  limiting  amplitudes  u  and  p  can  be  characterized  without  knowledge 
of  the  transients  TJ ,  FJ ,  and  acting  during  the  evolution  from  the  quiescent  state.  Perhaps  surprisingly, 
this  is  evidently  not  the  case;  in  general,  one  must  know  TJ  and  FJ,  although  is  not  required.  We  shall 
return  to  this  matter  shortly. 

It  is  necessary  to  proceed  through  the  examination  of  the  original  space-time  formulation  to  identify 
the  correct  time-independent  problem  to  be  studied  for  the  determination  of  the  limiting  amplitudes  of  the 
response,  essentially  because  there  is  a  possible  loss  of  the  unique-solution  property  in  a  too-na'ive  passage 
from  the  evolution  problem  to  the  time-harmonic  setting.  One  must  retain  the  unique-solution  property  in 
such  a  way  eis  to  force  the  solution  of  the  steady-state  problem  to  provide  precisely  the  limiting  amplitudes 
sought;  when  this  is  done,  we  naturally  refer  to  the  resultant  formulation  as  a  limiting- amplitude  problem. 
But  this  difficulty  is  not  of  an  unfamiliar  sort.  Indeed,  the  situation  is  analogous  to  that  obtaining  in  the 
simpler  case  of  acoustic  radiation  or  scattering  in  the  exterior  of,  say,  a  rigid  body,  wherein  the  Sommerfeld 
radiation  condition  (cf.  (1.34),  infra)  is  the  appropriate  requirement  to  be  imposed  for  the  purpose  of  picking 
out,  from  among  all  solutions  of  the  Helmholtz  equation  in  that  possess  the  required  Neumann  (normal- 
derivative)  data  on  r+,  precisely  that  function  that  represents  the  amplitude  of  the  time-harmonic  motion 
asymptotically  approached  from  the  quiescent  state.  Thus,  the  Sommerfeld  condition  can  be  referred  to  the 
“(exterior)  limiting-amplitude  condition”  for  that  case.  One  should  consult  in  this  regard  the  fundamental 
work  of  Wilcox  [1]  for  the  developments  pertinent  to  that  more  familiar  setting.  The  present  problem  is 
more  complex,  requiring  an  “interior  limiting-amplitude  condition,”  as  well  as  an  exterior  one.  Indeed,  by 


merely  setting  down  the  necessary  conditions  on  u  and  p  that  evidently  result  from  (1.10),  (1.11),  (1.13), 
and  (1.14)  in  consequence  of  the  (sufficiently  regular)  assumed  forms,  one  arrives  at  the  following  collection 
of  requirements  (cf.,  also.  Theorem  5.1  and  the  remarks  following  the  proof  of  Theorem  5.2,  infra): 

Aip  -(-  =  0, 

in 

(1.23) 

-f  QoUl^U  =  -fo. 

in  fio. 

(1.24) 

^.n+  -  »f>+wu|r+  •  n+  = 

T"+[u]  (7+  -  »w)^|r+n+  =  t+  -  (7+ 

>  onr+, 

-iu;)v?‘lr+n+,  J 

(1.25) 
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and 

T“-[u]  =  t_,  onr_.  (1.26) 

However,  one  can  already  see  that  this  problem  need  not  be  well  posed,  in  particular,  need  not  possess  at 
most  one  solution.  For  example,  suppose  that  there  is  a  nontrivial  field  such  that 


with 

and 


Uu,  -h  QgLj'^Uu,  =  0, 

in  Clg, 

(1.27) 

Ua,lr+  •»+  =  0, 

on  r4.. 

(1.28) 

T“»K]  =  0, 

on  Fj,. 

(1.29) 

Then,  obviously,  the  nontrivial  pair  (u  =  =  0)  satisfies  (123)  and  the  homogeneous  versions  of  (124) 

through  (1.26).  In  fact,  for  certain  domains  f2o  and  frequencies  ui,  there  do  exist  nontrivial  satisfying 
(1.27)  through  (1.29).  When  fio  is  a  ball  (so  fl_  is  absent)  or  a  spherical  annulus  (so  is  a  ball),  such  fields 
can  be  constructed  by  means  of  separation-of-variables;  cf.  Love  [2]  for  the  former  case.  The  eigenfunctions 
Ut,j  in  question  for  either  of  these  spherical  geometries  correspond  to  vibrations  in  which  the  particles  of 
the  body  move  on  spheres  concentric  with  the  boundary.  Some  of  these  motions  are  “purely  rotatory,” 
with  all  particles  oscillating  along  arcs  of  circles  lying  in  planes  normal  to,  and  having  centers  on,  a  single 
axis  passing  through  the  center  of  the  body;  it  is  reasonable  to  anticipate  that  modes  of  this  simple  type 
can  be  supported  by  any  body  of  revolution,  although  we  have  not  yet  worked  through  the  analysis  of  the 
pertinent  eigenvalue  problem.  For  other  classically  familiar  geometries  in  which  separation-of-variables-type 
arguments  are  applicable,  see  Ref.  3,  from  which  the  appropriate  computations  can  be  constructed  on  the 
basis  of  a  reformulation  of  the  eigenvalue  problem  governing  the  purely  rotatory  modes  for  a  certain  class  of 
shapes  of  revolution.  For  example.  Ref.  3  provides  some  of  the  details  necessary  in  the  study  of  this  problem 
for  ellipsoids  of  revolution.  It  would  be  very  useful  to  prove  (as  intuition  suggests)  that  the  class  of  bodies 
of  revolution  is  precisely  the  collection  of  domains  iio  for  which  such  nontrivial  eigenfunctions  can  exist, 
since  the  entire  development  that  we  shall  give  is  much  simplified  when  it  is  known  that  fields  of  this  type 
cannot  be  supported  by  the  particular  flo  under  consideration.  In  the  absence  of  such  a  proof,  and  since,  in 
any  event,  we  wish  to  establish  results  that  shall  be  valid  for  the  important  class  of  bodies  of  revolution,  we 
shall  proceed  under  the  assumption  that  there  may  be  nontrivial  Uu,  satisfying  (1.27)  through  (1.29)  for  the 
chosen  flo  and  w.  Such  a  field  xi^  can  be  termed  the  “complex  amplitude  of  a  nonradiating  mode,”  since  the 
corresponding  elastic-displacement  field  given  by 


Uu,(-,<)  ;=u,^()e-‘“'‘  infi^forteR  (1.30) 

then  churns  about  in  the  elastic  medium  in  but  has  no  effect  on  the  inviscid  fluid;  it  is  completely 
uncoupled  from  the  fluid  at  the  interface  r+,  and  so,  in  particular,  in^'aces  no  outward  radiation  of  acoustic 
energy.  No  measurement  made  in  the  fluid  will  serve  to  detect  its  presence  in  the  elastic  body. 

Conceivably,  there  are  more  complex  situations  in  which  ( 1 .23)  through  (1  26)  do  not  suffice  to  determine 
at  most  one  pair  (u,v?).  To  see  how  such  a  situation  might  come  about,  let  us  suppose  that  Qg  is,  say,  a 
Lipschitz  domain  with,  for  simplicity,  fi_  absent  (so  r„  =  r+),  and  the  positive  w  is  such  that  the  second 
fundamental  boundary-value  problem  of  steady-state  elastic  vibrations  of  angular  frequency  u)  for  Qg,  A, 
and  Qo  is  uniquely  solvable  in  a  strong  sense  for  any  boundary  data  chosen  from  some  linear  space  H(r4.)  of 
C^-valued  functions  defined  almost  everywhere  (a.e.)  on  r+.  Thus,  we  assume  that  whenever  t  6  H(r.(.), 
there  exists  a  unique  corresponding  field  Ut  in  the  Sobolev  space  H^(flo)  ;=  H^{Qg)^  satisfying 

A);  ^Ut  +  =  0  a.e.  in  flo  (1-31) 
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and 

T“‘’[ut]  =  t  a.e.  onr+,  (1-32) 

the  latter  condition  being  fulfilled  in  the  Sobolev-trace  sense;  for  more  on  the  notations  etc., 

and  terminology  being  used  here,  cf.  §2,  infra.  Note  that  H(r+)  must  lie  in  the  fractional-order  Sobolev 
space  H3(r+)  ;=  /f3(r+)^  associated  with  the  boundary  r+.  Parenthetically  we  note  that,  by  interior 
elliptic-regularity  results,  each  Ut  must  have  real-analytic  real  and  imaginary  parts,  and  so  (1-31)  is  actually 
satisfied  at  each  point  of  fi;,.  The  trace  Utlr+  is  then  automatically  determined  as  an  element  of  H2(r+)  by 
t,  so  that  there  is  defined,  as  a  result  of  the  assumed  strong  unique  solvability,  an  injective  linear  operator 

S,,  :H(r+)-Hi(r+): 

B,^t  :=  Ut|p^  (=  B,^T"+[ut])  whenever  t  G  H(r+). 

Now  suppose  that  there  exists  a  nontrivial  function  tp  G  that  possesses,  in  some  reasonable  sense,  a 

trace  and  a  normal  derivative  in  the  Lebesgue  space  Z,2(r^.)  {e.g.,  suppose  that  p  is  in 
and  interpret  the  trace  and  normal  derivative  in  the  Sobolev  sense),  with  i^|r+n+  G  H(r+),  and  is  a  solution 
of  a  nonlinear  and  nonlocal  eigenvalue  problem  for  —A  in  0+,  specifically,  that  it  satisfies 


A<^ -h /c^<^  =  0,  in 


along  with  the  homogeneous  (nonlocal)  boundary  condition 

<^5.0+  +  ■  Bw(^lr+n+)  =  0.  (1.33) 

Then,  taking  u  :=  uj,  wherein  t  ;=  —(7+  —  tw)9plr+n+,  we  obtain  a  nontrivial  pair  (u,  1^)  satisfying  (1-23) 
and  the  homogeneous  versions  of  (1-24)  through  (1.26),  as  can  be  checked.  One  would  expect  that  this  sort 
of  eventuality  could  be  ruled  out  with  the  additional  imposition  of  the  Sommerfeld  condition.  We  remark 
that  a  generalized  exterior  Robin  problem  for  (1.23)  and  a  boundary  condition  of  a  form  corresponding 
to  the  nonhomogeneous  version  of  (1.33)  is  studied  in  Refs.  4  and  5  (Ref.  5  being  a  revised  and  somewhat 
amplified  version  of  Ref.  4),  under  certain  restrictions  on  the  operator  appearing  in  the  boundary  condition. 
In  that  case,  it  is  found  that  the  Sommerfeld  condition  does  serve  to  prohibit  the  existence  of  a  nontrivial 
eigenfunction  of  the  sort  that  we  have  just  described. 


At  any  rate,  it  is  evident  that  some  conditions  augmenting  (1.23)  through  (1.26)  are  wanted  before  there 
is  a  chance  that  the  resultant  problem  will  be  well  posed.  Keeping  in  mind  the  hypothesized  origin  of  the 
fields  being  sought  here,  we  should  demand  that  these  additional  requirements  also  follow  from  the  original 
initial- value  problem,  as  conditions  necessarily  fulfilled  by  assumed-to-exist  complex  limiting  amplitudes  of 
the  elastic  and  fluid  perturbation  fields;  cf.  the  assumptions  on  the  forms  in  (1.21)  and  (1.22).  To  see  more 
specifically  what  is  involved,  let  us  formulate  two  requirements  that  we  might  consider  imposing  in  addition 
to  (1.23)  through  (1.26),  viz.,  for  the  exterior  field,  the  Sommerfeld  radiation  condition. 


lim  r{e  ■  grad  (^(re)  —  j/<^.yj(re)}  =  0  uniformly  for  e  G  ,  (1-34) 

r— *00 

with  Si'  denoting  the  boundary  of  the  unit  ball  in  R^,  and,  for  the  interior  field,  the  orthogonality  condition 
u  u,^dA3  =  0  for  every  Uu,  satisfying  (1.27)  through  (1.29),  (1.35) 


L 


an  overbar  signifying  complex  conjugation  and  A3  denoting  the  Lebesgue  measure  on  R^;  the  precise  sense(s) 
in  which  (1.27)  through  (1.29)  are  to  be  satisfied  in  (1.35)  (and  in  (1.36),  infra)  shall  be  specified  later.  Now, 
one  can  show,  as  we  shall  in  §2,  that  the  problem  of  finding  a  pair  {u,<p)  satisfying  (1.23)  through  (1.26), 
(1.34),  and  (1.35)  can  have  at  most  one  solution  in  a  reasonably  large  class  of  pairs,  each  comprising  an 
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interior  vector  field  and  an  exterior  scalar  field.  But  this  is  not  sufficient  to  fulfill  our  stated  objective; 
we  must  also  decide  whether  (1-34)  and  (1.35)  are  limiting- amplitude  conditions,  i.e.,  whether  they  are 
necessarily  satisfied  by  our  hypothesized  complex  limiting  amplitudes  figuring  in  the  lime-harmonic  states 
approached  by  fields  evolving  from  a  quiescent  state  under  the  influence  of  driving  agents  as  in  (115)  through 
(1.17).  In  fact,  (1.35)  does  p  .  in  general  meet  this  criterion.  Rather,  evidently  one  should  require  that 

/  u-u,,dA3=  r(  /  Ta-,s)u,.|r^dAr„+  /  F;(-,  s)  ■  u„  dAal  ds 

Jcio  Jo  lJr„  Jci„  J 

for  every  satisfying  (1.27)  through  (1.29).  (1.36) 

One  of  the  goals  of  this  report  is  to  supply  motivation  for  the  claim  that  (1.34)  and  (1.36)  are,  respectively, 
the  correct  exterior  and  interior  limiting-amplitude  conditions  to  be  imposed,  by  proving  that  they  must 
obtain  at  least  under  sufficiently  strong  hypotheses  concerning  the  transients  in  (1.15)  through  (1.17),  (1-21), 
and  (1.22).  Once  this  has  been  established,  (1.36)  will  imply  that  the  interior  limiting  amplitude  u  is  not 
independent  of  the  particular  manners  by  which  asymptotic  approaches  to  steady-state  evolve  under  the 
influence  of  various  To  and  Fo  differing  only  in  their  transient  parts.  The  reason  for  this  circumstance 
resides  in  the  fact  that  the  transient  parts  of  To  and  Fo,  while  dying  away,  nevertheless  may  leave  a  vestige 
of  their  action  by  inducing  a  nonzero  L2(Oo)^-orthogonal  projection  of  the  resultant  limiting  amplitude 
u  onto  the  subspace  of  solutions  of  (1.27)  through  (1.29),  and  different  transient  components  may  induce 
different  projections.  It  is  trivial  to  see  that  this  is  so  for  the  case  in  which  To  =  0  and  Fo  =  FJ  ^  0.  While 
the  proof  of  the  assertion  in  the  case  To  =  TJ  ^  0  and  Fo  =  0  requires  much  more  work,  we  choose  not  to 
provide  the  reasoning  in  the  present  report. 

There  is  an  additional  consideration  necessitated  by  the  presence  of  the  forcing  tractions  To  and  body 
force  Fo.  We  describe  this  by  making  some  remarks  that  shall  remain  for  the  present  in  the  nature  of  a 
plausible  conjecture,  since  we  shall  not  prove  them  here.  Consider  for  a  moment  the  driving  of  the  elastic 
body  in  fio,  but  in  vacuo,  from  a  quiescent  state  by  means  of  (bounded!)  applied  tractions  and  body 
force  as  in  (1.15)  and  (1.16):  in  this  case,  when  w  has  any  one  of  an  infinite,  discrete  set  of  positive  values 
corresponding  to  the  eigenvalues  of  the  traction  problem  for  —  in  a  “resonance”  must  be  expected, 
I.e.,  one  or  more  of  the  pertinent  eigenmodes  will  be  excited  and  the  displacement  response  of  the  body  will 
not  be  bounded  in  Ug  x  (0,oo),  so  certainly  a  time-harmonic  state  will  not  be  approached  asymptotically. 
Now,  when  the  elaistic  body  is  surrounded  by  fluid,  evidently  the  situation  with  regard  to  possible  resonances 
is  significantly  altered,  for  one  would  expect  that  the  fluid  acts  as  an  energy  sink  into  which  the  body  suffers 
losses  by  working  at  the  interface.  Consequently,  even  when  w  corresponds  to  one  of  tne  eigenvalues  just 
specified,  it  can  be  anticipated  that  the  tendency  to  resonate  owing  to  excitement  of  a  traction  eigenmode 
will  be  suppressed  to  the  extent  that  the  response  will  remain  bounded  unless  the  body  possesses  such  a 
mode  that  can  exchange  no  energy  with  the  fluid,  w  corresponds  to  the  eigenvalue  for  that  mode,  and 
and  fo  fail  to  satisfy  a  certain  orthogonality  condition  relative  to  that  mode.  When,  as  in  our  case,  the 
fluid  is  inviscid,  the  special  eigenfunctions  in  question  are  just  those  (nontrivial)  u^,  satisfying  (1.27)  through 
(1.29)  for  some  w  >  0.  Formal  computations  based  on  Laplace  transformation  indicate  that  a  condition  on 
to  and  fo  that  is  necessary  and  sufficient  to  ensure  that  no  resonance  be  induced  by  excitement  of  any  such 
(nonradiating)  mode  appears  as 

I  to  •  +  /  fo-Uu>dA3  =  0  for  every  satisfying  (1.27)  through  (1.29).  (1-37) 

On  the  other  hand,  it  is  al.so  easy  to  use  certain  integral  identities  to  show  that  (1.37)  is  necessary  for  the 
existence  of  a  solution  to  the  problem  generated  by  (1-23)  through  (1.26),  (1.34),  and  (1.36),  while  the  same 
condition  will  appear  naturally  in  the  development  of  the  existence  assertion  of  §4.  Observe  that  (1-37)  is 
automatically  fulfilled  (and  (1.36)  reduces  to  (1.35))  when  the  only  forcing  agent  in  the  problem  is  the  fluid 


9 


ALLAN  G.  DALLAS 


incident  field.  While  the  conjectures  in  this  discussion  probably  cannot  be  substantiated  unless  the  domain 
fio  satisfies  some  additional  requirements  on  its  shape,  they  provide  an  heuristic  basis  for  anticipating  the 
need  for  (1.37). 

As  we  proceed,  it  becomes  clear  that  it  would  be  of  great  assistance  to  have  knowledge  of  conditions 
on  the  shape  of  sufficient  to  guarantee  that,  for  the  given  u;  >  0,  there  can  exist  no  nontrivial  fields  of 
the  sort  that  we  have  termed  “complex  amplitudes  of  nonradiating  modes”  for  and  w,  those  very  special 
(eigen )functions  satisfying  (1.27)  through  (1.29).  For,  when  we  are  certain  that  such  fields  do  exist  or 
when  we  lack  any  separate  assertion  of  their  nonexistence,  we  must  deal  at  every  step  of  the  analysis  with 
their  presence  or,  respectively,  possible  presence,  an  aspect  significantly  complicating  matters.  In  particular, 
this  is  true  in  the  numerical  work,  since  it  turns  out  that  if  one  wants  to  approximate  both  the  elastic  field 
u  and  the  acoustic  field  <fi  satisfying  (1-23)  through  (1.26),  (1.34),  and  (1.36)  (when  such  a  pair  exists),  then 
complete  information  concerning  the  solutions  of  the  problem  embodied  in  (1.27)  through  (1.29)  is  required. 
However,  we  also  wish  to  show  that  if  one  is  content  with  the  approximate  calculation  of  the  fluid  field  alone, 
then  one  can  make  do  with  very  meager  information  about  the  latter  problem. 

To  sum  up,  our  aim  is  threefold: 

1.  We  wish  to  establish  the  unique  solvability  of  the  problem  generated  by  (1.23)  through  (1.26),  (1.34), 

(1.36) ,  and  (1.37),  in  a  weak  sense  to  be  specified. 

2.  We  intend  to  show  how  the  Galerkin  method  can  be  used  to  construct  a  convergent  sequence  of  approx¬ 
imations  for  the  unique  solution  of  the  problem  just  cited  above,  provided  that  one  “knows  all  about” 
both  the  exterior  Neumann  problem  for  the  Helmholtz  equation  (1.23)  with  the  Sommerfeld  condition 
(1.34)  and  the  collection  of  solutions  of  the  problem  (1.27)  through  (1.29).  Also,  we  want  to  verify  that 
the  acoustic  part  of  the  solution  can  be  approximated  without  explicit  knowledge  of  the  latter  family  of 
functions. 

3.  We  wish  to  prove  that,  under  sufficiently  strong  hypotheses  concerning  the  assumed  asymptotic  approach 
to  steady-state,  (1.34)  and  (1.36)  are  properties  of  the  complex  amplitudes  u  and  ip  in  the  Eissumed  forms 
(1.21)  and  (1.22)  of  the  solution  (U,$)  of  (1.10)  through  (1.14),  with  forcing  data  having  the  forms 
appearing  in  (1.15)  through  (117).  The  results  here  are  less  than  satisfactory,  in  that  they  provide  only 
a  “plausibility  argument,”  on  the  basis  of  which  we  are  motivated  to  suspect  strongly  that  (1.34)  and 

(1.36)  are  genuine  limiting-amplitude  conditions. 

We  shall  carry  out  the  analysis  in  this  same  order.  Thus,  in  §2  we  shall  reduce  the  study  of  the  original 
interior-exterior  problem  to  the  examination  of  a  purely  interior  problem  (in  Qo),  while  giving  strong  and 
weak  formulations  for  both  the  original  and  the  purely  interior  problems.  Subsequently,  we  shall  restrict 
our  attention  to  the  weak  problems,  deferring  to  a  future  report  a  discussion  of  regularity  results  for  the 
weak  solutions,  in  particular,  postponing  the  specification  of  conditions  sufficient  to  guarantee  that  a  weak 
solution  is  a  strong,  or  <•  .  a  claissical,  solution.  The  abstract  results  that  we  need  for  the  analysis  of 

the  weak  purely  interim:  j>  Mem  are  presented  in  §3.  There,  we  base  our  approach  on  developments  of 
Hildebrandt  and  Wienl.on./.  [61  generalizing  the  Lax-Milgram  Lemma,  that  we  either  take  over  entirely  or 
modify  to  suit  our  own  purp^ —  ‘  In  particular,  an  evidently  new  result  is  given  in  Theorem  3.3,  identifying 
a  simple  sufficient  co  lition  inder  which  the  Galerkin  procedure  can  be  directly  applied  to  a  variational 
problem  with  an  indefinite  sesquilinear  form  appearing  as  a  perturbation  of  a  definite  form  by  a  compact 
one.  This  machinery  is  applied  in  §4  to  accomplish  the  first  two  aims  listed  above.  In  §5  we  undertake  to 
provide  the  promised  motivation  for  accepting  (1.34)  and  (1.36)  as  limiting-amplitude  conditions.  Finally, 
the  proofs  of  two  auxiliary  results  used  in  §5  are  given  in  the  Appendix. 
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2.  STRONG  FORMULATION;  DERIVATION  OF  THE  PURELY  INTERIOR  PROBLEM 
AND  WEAK  FORMULATIONS 

Following  a  description  of  the  basic  notations  to  be  used  throughout,  we  proceed  directly  in  this  section 
to  the  formulations  and  reformulations  of  the  problem  represented  by  (1.23)  through  (1.26),  (1.34),  and 
(1.36).  The  plan  consists  in  the  use  of  certain  facts  about  the  Neumann  problem  for  the  Helmholtz  equation 
in  0+  to  reduce  the  original  interface  problem  to  one  involving  only  the  eleistic  field  u  in  flo,  with  subsequent 
construction  of  the  fluid  field  tp  in  from  u. 


We  adhere  to  the  standard  notations  for  the  usual  spaces  of  complex  functions  on  subsets  of  R^ 
(N  >  2).  Thus,  C°(F)  =  C(F)  indicates  the  complex-linear  space  of  continuous  C-valued  maps  defined  on 
F  C  R^.  Let  Q  be  open  in  R^;  for  a  positive  integer  k,  C'‘{Cl)  is  the  linear  manifold  of  all  elements  of  C(n) 
having  all  partial  derivatives  of  orders  less  than  or  equal  to  k  also  belonging  to  C(r2).  When  /  G  (7*^(0), 
in  addition  to  the  notation  /,j,  we  may  also  indicate  partial  derivatives  of  /  in  an  obvious  manner  as  /,<,, 
wherein  a  =  (o^ , . . . ,  a-y^)  is  an  V-index  of  order  |a|  :=  aj  -t- . . .  -f  <  it.  The  elements  of  C^(Q)  having 
supports  compact  and  contained  in  Q  is  the  linear  manifold  denoted  by  Co(n).  The  members  of  C*(f2)  are 
those  functions  that  lie  in  C*(0)  and  are,  along  with  all  of  their  partial  derivatives,  bounded  and  uniformly 
continuous  (and  so  possess,  with  all  of  their  partial  derivatives,  continuous  extensions  to  Q);  C*(fl)  is  a 
Banach  space  under  the  usual  norm,  given  by 


In  turn,  for  a  number  i/  with  0  <  t/  <  1,  C^'‘'(Q)  is  the  linear  manifold  in  C*(D)  comprising  those  functions 


that  are,  along  with  all  of  their  partial  derivatives,  uniformly  j/-H61der  continuous  in  ft;  provided  with  the 
norm  given  by 


'C*(0)  0<|cr|<t  x.y^n 


l/,a(x)  -  /,c»(y)| 
Ix-yl" 


C*’‘'(n)  is  a  Banach  space.  The  spaces  C°°(Q),  Co°(Q),  and  are  defined  as  the  intersections  of, 

respectively,  all  C'‘(Q),  Cq{Q),  and  C*(n),  for  ^  =  1, 2, 3, . . .  . 


Whenever  /F  is  a  measure  space,  with  positive  measure  m,  L2('^)  >s  the  collection  of  all  (equivalence 
classes  of)  complex  m-measurable  functions  /  defined  m-a.e.  on  X  with  dm  <  oo;  in  this  setting, 

i2(T)  is  equipped  with  the  inner  product  (•, 


{f’9)L3{X) 


for  f,g  £  L2{X), 


under  which  it  is  a  Hilbert  space.  The  Lebesgue  measure  on  R^  shall  be  denoted  by 

Concerning  regularity  hypotheses  for  a  bounded  and  connected  open  subset  Q  of  R^,  we  use  the  defini¬ 
tions  of  Necas  [7].  Corresponding  to  a  positive  a,  let  Sa  denote  the  open  square  (—a,  a)  x  (—a,  a)  in  R^.  By 
an  affine  isometry  of  R®,  we  mean  a  function  x  ►  Lx  =  x,,  +  L„x  on  R^,  with  Xo  G  R^  and  L<,  :  R^  — »  R^ 
a  linear  isometry.  Now,  let  k  be  either  a  nonnegative  integer  or  oo;  Q  is  said  to  be  of  type  provided  there 
exist  positive  numbers  a  and  6,  a  positive  integer  n,  a  collection  of  real-valued  functions  contained 

in  C^(Sa),  and  a  family  {L/}"_,  of  surjective  affine  isometries  of  R^  such  that 
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(a)  for  /  =  1, . . . ,  n,  L(  maps  the  graph  C{/|}  of  /,  into  d€i,  the  set 

{(a;i,X2,X3)  1  (xi.Xj)  €  5a, //(Xl,  *2)  <  X3  <  //(xi,X2)  +  6} 

into  n,  and  the  set 

{(xi,X2,X3)  1  (Xi,X2)  €  5a,/,(xi,X2)  -  6  <  X3  <  /,(xi,X2)  } 

into  \  n, 

and 

(b)  5n  =  ur=iLf(^{/,})- 

Further,  if  these  conditions  are  fulfilled  and  the  family  is  in  C'‘’‘'{Sa)  for  some  number  //  6  (0, 1], 

then  n  is  said  to  be  of  type  3?*'".  Bounded  and  connected  open  sets  of  type  are  also  termed  Lipschitz 
domains.  For  the  definition  of  the  Lebesgue  measure  (induced  by  A3)  on  the  boundary  90  of  a  Lipschitz 
domain  O,  cf.  Ref.  7,  §3.1.1  or  Ref.  8,  §1.1;  we  shall  denote  this  meeisure  by  A^^j,  and  its  restriction  to  the 
measurable  subsets  of  a  measurable  subset  F  of  90  by  Ap.  For  such  O,  the  unit  exterior  normal  field  n 
exists  A^fj-a.e.  on  90,  with  Cartesian  components  A^j^-measurable  (and  bounded);  by  “exterior,”  of  course 
we  mean  here  that,  whenever  a  normal  is  defined  at  x  €  90,  the  inclusion  x  +  sn(x)  €  R^  \  0  obtains  for 
all  sufficiently  small  positive  s.  Explicitly,  suppose  that  x  £  90  is  such  that  for  some  I  £  {1, . .  .,n}  and 
(xi,X2)  £Sa  we  have  x=  L, (xj,  X2, /,(xi,  X2)):  then,  provided  that  /,,i(xi,X2)  and  /,,2(xi,X2)  exist,  n(x) 
exists  and  can  be  computed  from 

n(x)  =  L,„(n,(xi,X2)), 

wherein  L,„  is  the  surjective  linear  isometry  associated  with  L,,  and  (in  Cartesian  components) 

n(x  i  1  •= 

'  ■  {l  +  [/m(^i.i2)P  +  [/n2(^i.^2)n''' 

is  the  appropriate  unit  normal  to  the  graph  of  //  at  the  point  (xi,X2,/((xi,X2)). 

For  the  definitions  and  properties  of  the  various  Sobolev  spaces,  we  rely  on  Necas  [7];  cf.,  also,  Adams  [9], 
wherein  certain  of  the  results  are  to  be  found.  Returning  to  the  case  in  which  fi  is  open  in  R^,  for  a  positive 
integer  k  the  corresponding  Sobolev  space  (=  lF^*\n)  in  Ref.  7,  =  VF*'^(n)  in  Ref.  9)  is  the  set  of 

all  elements  of  H°{Q)  :=  L2{^)  for  which  all  weak  (distributional)  derivatives  cf  orders  less  than  or  equal  to 
k  exist  and  belong  to  Z.2(f^);  we  extend  the  notations  f,j,  f,a  to  indicate  the  weak  derivatives  of  an  element 
/  £  Defining  an  inner  product  on  i/*(D)  by 

(/’<»>9''a)L3(n)  for /,  y  G /f*(n), 

0<|a|<lfc 


i/*(0)  becomes  a  Hilbert  space.  The  definition  of  ff*(f2)  is  extended  to  encompass  nonintegral  values 
of  A:  >  0  as  in  Ref.  7,  §2.3.8,  in  a  well-known  manner,  by  using  the  approach  of  S.  L.  Slobodetskii  and 
N.  Aronszajn;  such  spaces  are  needed  here  for  the  definition  of  the  fractional-order  Sobolev  spaces  on  9f2, 
infra.  Thus,  when  s  ir  positive  d  nonintegral,  with  [s]  denoting  the  greatest  integer  not  exceeding  s,  one 
sets 

f  ^  2 


|a|  =  [«l 


for  those  u  £  //t’l(f2)  for  which  the  integrals  appearing  on  the  right  are  all  finite;  the  resultant  complex-linear 
space  is  equipped  with  the  obvious  inner  product  generating  {|  ■  forming  a  Hilbert  space  denoted  by 

H’(Q).  Finally,  when  Q  is  not  bounded,  by  we  shall  mean  the  collection  of  all  complex  measurable 
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functions  defined  a.e.  in  Q  and  having  restrictions  in  for  every  bounded  open  subset  ft'  C  in  the 

standard  manner,  is  made  into  a  locally  convex  linear  topological  space  through  the  introduction  of 

the  seminorms  u  i-* 

In  this  paragraph,  il  C  R.^  shall  be  (at  least)  a  Lipschitz  domain,  whence  a  good  deal  more  is  known 
about  the  structure  of  than  in  the  case  of  a  general  open  set.  We  shall  write  T  :=  dQ.  Thus,  C*(fl) 

is  dense  in  H^{Q)  (cf.  Ref. 7,  Theoreme  2.3.1  or  Ref. 9,  Theorem  3.18),  and  there  is  defined  the  Sobolev 
trace  operator  u  i— ►  ujp  from  H^(Q)  into  f/2(r),  as  the  unique  bounded  extension  to  H^{Q)  of  the  “classical” 
trace  map  /  i— >•  f\p  ;=  /jF  for  /  G  C'(fl),  regarded  m  densely  defined  in  ff*(f2)  and  mapping  into  7/2(1’) 
(Ref.  7,  Theoreme  2.4.2).  Moreover,  both  the  natural  injection  of  into  L^iQ)  and  the  trace  operator 

carrying  H^{Q)  into  i2(r')  are  compact  (Ref.  7,  Theoreme  2.6.1  and  Theoreme  2.6.2).  At  various  points, 
we  shall  make  reference  to  Sobolev  spaces  7f’(r)  associated  with  F,  for  certain  real  (nonintegral  as  well 
as  integral)  values  of  s.  Following  §2.4.3  and  §2.5.2  of  Ref.  7,  let  /b  be  a  positive  integer  and  fi  be  of  type 
0  <  s  <  I:,  the  corresponding  7f*(r)  is  defined  to  comprise  those  ff  in  L2(r)  for  which  the 
function  g/  defined  A2-a.e.  on  Sa  by  ^,(11, X2)  :=  g  (l>;(ii,  X2, //(xi,  X2)))  lies  in  i/'(5o)  for  /  =  1, . . . ,  n  (cf 
the  notations  employed  in  the  definitions  of  the  regularity  classes  of  open  sets,  supra)-,  when  provided  with 
the  inner  product  ensuing  from  this  construction. 


n 

1=1 

/f*(r)  forms  a  Hilbert  space.  With  this  definition,  the  range  of  the  trace  operator  on  7/' (ft)  is  characterized 
as  precisely  the  collection  of  elements  in  and  the  resultant  map,  when  regarded  on  ff*(ft)  onto 

/f  ^(F),  is  bounded  (Ref  7,  Theoreme  2.5.5  and  Theoreme  2.5.7).  Finally,  provided  that  s  >  0  and  .H'(F)  is 
defined,  denotes  the  anti-dual  of  //’(F)  (the  collection  of  all  bounded  conjugate-linear  functionals 

on  H‘{r))  and  (•,  •),  the  resultant  duality  pairing  on  //“’(F)  x  //’(F).  With  the  usual  identification,  we  can 
regard  H°{T)  :=  /-2(F)  as  a  linear  manifold  in  //~’(F),  and  then  can  write 

(/,g).  =  J^fgdXr  =  (/,s)L,(r)  for  /  €  H°{r)  and  g  G  //’(F). 

The  extension  of  the  formula  for  integration  by  parts  (“formule  de  Green”),  classically  a  corollary  of  the 
Divergence  Theorem,  to  the  Sobolev-space  setting  for  a  Lipschitz  domain  is  effected  in  Theoreme  3.1.1  of 
Ref.  7.  Thus,  we  have 

/  u,jvdXz—  I  u|ru|pnj  dAp-  /  uv,j  d A3  for  u,  u  G //^f^),  j  =  1,2,3, 

Jn  Jr  Jn 

(ni,  02,03)  denoting  Cartesian  components  of  the  unit  exterior  normal  field  n,  defined  Ap-a.e.  on  F. 

For  each  of  the  spaces  of  complex  functions  defined  in  the  preceding,  we  introduce  the  corresponding 
space  of  C^- valued  functions,  identified  as  the  Cartesian  product  of  three  copies  of  the  original  collection, 
and  denote  the  new  space  with  the  corresponding  boldface  symbol,  C(F'),  C*'(ft),  L2(/V),  H*'(ft),  etc.  Ij2(X) 
is  provided  with  the  Hilbert-space  structure  induced  by  the  inner  product  given  as 

(f.g)^.^)  =  g dm  for  f,gG  L2(A’); 

H^(ft)  is  equipped  with  the  inner  product  defined  by 

(f>a'g-<»)Mn)  for  f,gG  H*'(ft), 

0<|a|<t 
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under  which  it  also  becomes  a  Hilbert  space.  Similarly,  the  other  products  of  Hilbert  spaces  are  given  their 
natural  Hilbert-space  structures. 

Returning  to  the  developments  begun  in  §1,  we  retain  the  notations  already  introduced  there.  We 
now  suppose  that  Qo  is  a  Lipschitz  domain,  but  we  impose  the  restriction  that  (the  fluid-elastic  interface) 
r+  be  a  manifold  of  class  C^,  just  so  that  we  can  use  directly  the  results  of  Refs. 4  and  5,  wherein  this 
much  smoothness  v/as  demanded.  Still,  the  components  of  (the  vacuum-elastic  boundary)  r_  are  allowed 
to  be  nonsmooth,  with  corners  and  edges  of  a  severity  permitted  within  the  Lipschitz-domain  setting.  More 
precisely,  then,  the  (bounded  and  connected)  open  set  Ufto  shall  be  taken  throughout  to  be  of  type  3?^. 
Alternately,  the  regularity  requirements  on  r+  can  be  phrased  in  terms  of  the  notations  already  introduced: 
let  a  >  0,  6  >  0,  the  integer  n,  and  the  families  and  {LjjJLj  be  associated  with  Qo  as  a  result  of 

its  being  of  type  If  n_  is  absent,  let  ;=  n;  otherwise,  we  may  suppose  that  matters  are  arranged 

so  that  {1, . . . ,  n}  can  be  partitioned  into  {1, . . . ,  n.j.}  and  {n++  l,...,n},  with  r+  =  (J"^i  and 

r-  =  L,  (^{/,})  Then  our  restriction  is  expressed  by  the  inclusion  C  C^{Sa)-  Observe  that 

the  spaces  /f*(r+)  and  H'’(r+)  are  defined  for  |s(  <  2  and  can  be  described  more  explicitly  with  the  help  of 
a,  6,  71^.,  and  the  /,  and  L,  for  /  =  1, . .  .,n+.  For  brevity,  we  write  (with  H°(no)  -■=  I'2(Oo))  H*  ;=  H*(no) 
for  t  =  0, 1, 2, . . .  .  If  u  €  by  u(p_  and  u|p^  we  mean  the  restrictions  to  F-  and  r+,  respectively,  of  the 
trace  u|p^  of  u  on  F,,,  i.e., 

u|r_  -  (u|ro)|r-  and  u|p^  :=  (u|pj|r+  whenever  u  e  H*. 

With  this  understanding,  we  want  to  make  sure  of  some  simple  facts  concerning  the  linear  maps  u 
and  u  I— ►  u|p^  •  n+  on  ,  which  will  be  of  importance  in  the  later  developments. 

Lemma  2.1.  Recall  the  regularity  conditions  imposed  on  Qp. 

(i)  The  linear  operations  u  ^  ujp^  and  u  *-♦  ujp^  •  n+  on  are  compact  when  regarded  as  mapping  into 
L2(r+)  and  L2(r+),  respectively. 

(it)  The  linear  operations  u  u|p^  and  u  u|p^  11+  carry  H*  onto  H » (r+ )  and  into  H » (r+),  respectively, 
and  are  bounded  when  so  regarded. 

Proof  :  The  assertions  in  (i)  follow  immediately  from  the  compactness  of  the  trace  map  u  i-->  u|p^  on  into 

L2(ro)  and  the  obvious  inequalities  ll(gir+)  n+||^^(P^)  <  llgir+l|i,^(P^)  <  ||glli,,(p^),  holding  for  g  e  L2(r„). 

Turning  to  statement  (ii),  let  us  first  show  that  H^r+)  =  {  u|p^  1  u  6  }.  The  inclusion  “D”  here  is 

easily  seen  to  hold  by  virtue  of  the  inclusion  u|p^  €  H^r<,)  and  the  remarks  made  concerning  the  generation 

of  H3(r+).  The  proof  of  the  opposite  inclusion  requires  not  much  more  work;  taking  any  g  6  H^r+), 

extend  its  definition  (if  necessary)  to  get  an  element  g  of  H^(r<,)  by  setting,  say,  g(x)  :=  0  for  x  G  r_ ;  there 

exists  some  u  G  such  that  ujp^  =  g,  whence  we  have  ujp^  =  g.  Thus,  u  1— ►  u|p^  carries  H’  ontoHi(r+); 

the  boundedness  of  this  “restricted-trace”  map  then  follows  from  the  boundedness  of  the  trace  operator  from 

H'  onto  H3(ro)  and  the  evident  inequality  ||glr+||  ^  <  l|g||  j.  ,  holding  for  g  G  H2(ro).  Finally, 

H3(r+)  H2(ro) 

let  u  G  H*:  we  show  that  u|p^  ■n+  lies  in  and  that  there  exists  a  positive  C],  independent  of  u,  such 

that  llulr  n+ll  ,  <  cillulr  II  1  ,  which  then  effectively  completes  the  proof  of  (ti).  To  this  end, 

by  using  the  notation  introduced  above,  select  any  I  G  {!,...,«+}  and  consider  the  function  (u|p^  '  n+), 
defined  for  A2-almost  all  (xi,ir2)  G  Sa  by 

(u|r+  ■  ri+),(xi,X2}  :=  (u|r^  -n+j  (L,(ari,a-2,//(^i, 

=  u(r^  (L,(xi,X2,/,(xi,X2)))  ■L,„(n,(xi,X2))  =  (u|p_^),(xi , X2)  ■  L,„(n,(xi, X2)). 
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Clearly,  the  inclusion  n,  G  C^(5o)  holds  by  our  restriction  on  flo,  so  the  function  n,  is  certainly  uniformly 
Lipschitz  continuous  on  So-  Consequently,  we  get 

f  f  l(«lr+ •n+)i(^i.^2)-(u!r+ •n+),(yi,y2)P  ,,  ,  , 

JsJs. - (».,») 

<  2{/ 


y2)P 


£iA2(a:i,X2)dA2(yi 


-1/2)  I 


in  which  M,  and  Ml  depend  upon  flo  and  I  alone.  Moreover,  obviously 
ll(^lr+  ■  “+)/ll/f<>(5.)  ^  ll(“lr+)jllHO(s.)  ^ 

From  these  facts,  we  conclude  that  (u|p^  '^+)t  Jl^s  in  H3(Sa),  with  norm  not  greater  than  a  positive  multiple 
)’  multiple  depending  upon  only  Oo  and  1.  In  turn,  this  implies  that  u|p^  11+  belongs 

to  H^{r+),  with  norm  not  exceeding  a  positive  multiple  of  ||u|p^||^j^^  the  multiple  depending  on  fl,, 
alone.  ■ 


For  the  various  formulations  of  the  problem  loosely  described  in  §1,  and  for  its  reformulation  as  a  purely 
interior  problem  in  we  need  some  preliminary  results  on  the  (purely  exterior)  Neumann  problem  for  the 
Helmholtz  equation  and  the  Sommerfeld  condition  in  0+,  with  boundary  data  lying  in  L2(r+).  For  this, 
it  is  most  convenient  to  rely  on  the  developments  already  in  place  in  Refs.  4  and  5.  There,  a  generalized 
Robin  problem,  subsuming  the  Neumann  problem,  is  studied  by  using  potential-theoretic  methods,  although 
the  analysis  is  not  carried  out  within  the  Sobolev-space  setting,  but  completely  in  £.2(r+);  the  later  article 
Ref.  10  is  concerned  with  the  connections  between  the  two  approaches.  In  particular,  in  Refs.  4  and  5,  traces 
and  normal  derivatives  on  r+  are  taken  in  the  normal-L2  sense.  To  recall  the  pertinent  definitions,  let 
N,  :  r+  be  given  for  any  s  >  0  by 

N,(x)  ;=  X sn+(x)  for  each  xGr+; 

because  of  the  regularity  imposed  upon  Q<,,  there  is  some  s+  >  0  such  that  N,(r4.)  C  n+  whenever 
0  <  s  <  s+.  Thus,  if  tp  is  any  function  defined  in  Q+,  then  o  N,  is  defined  on  r+  for  all  sufficiently  small 
positive  s.  Now  suppose,  for  simplicity,  that  tp  G  C^(0+)-  Then  we  say  that  ip  has  a  trace  on  r+  in  the 
normal-L2  sense  iff  lim,_o+  po'N,  exists  in  the  L2(r+)-sense,  i.e.,  iff 

^liin  (1^0  N,  -  v2|r+llMr+)  =  0 

for  some  G  L2(r+),  in  which  case  the  latter  function  is  termed  the  normal-L2  trace  of  p  on  r+. 
Similarly,  we  say  that  p  has  a  normal  derivative  on  r+  in  the  normal- L2  sense  iff  lim,_o+  n+  •(gradi^)oN, 
exists  in  the  L2(r+)-sense,  i.e.,  iff 

lim  |jn.+  •  (grad  p)oN,-  ||t,(r+)  =  0 
«-*0+  ^ 
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for  some  <p,a^  €  I/2(r+),  in  which  case  the  latter  function  is  called  the  normal-L^  normal  derivative  of  (p  on 
r+.  When  both  <p\p^  and  exist,  we  say  that  tp  is  L2-regu/ar  at  r+.  With  the  understanding  that  the 
trace  and  normed  derivative  of  an  appropriate  function  (p  in  Cl+  are  always  to  be  interpreted  in  this  sense, 
henceforth  we  omit  the  superscript  “2”  and  write,  respectively,  simply  p\r+  and  p,^  for  these  elements  of 

i2(r+). 

Now,  we  introduce  the  collection  W(n+;  /c+)  of  radiating  solutions  of  the  Helmholtz  equation  (1.23)  in 
that  are  also  L2-regular  at  r+; 

/c+)  :=  {p  £  C^(Q+)  |  (1.23)  and  (1-34)  hold,  and  p  is  Z2-regular  at  r+  }. 

For  p  6  W(n+;/c+),  we  refer  to  </>lr+  and  respectively,  as  the  Dirichlet  data  and  the  Neumann  data  of 
f- 

The  basic  facts  that  we  shall  use  concerning  the  Neumann  problem  for  (123)  and  (134)  in  are 
collected  in  Theorem  2.1.  To  state  the  results,  we  need  to  introduce  an  “outgoing”  fundamental  solution  for 
the  operator  A  +  we  choose  the  particular  one  given  by,  for  each  x  G  R^, 

foryGR3\{x}. 

Whenever  x  G  R^,  we  write  £^xtn+  •=  ^*‘'’)|(r+ \  {x})  for  the  normal  derivative  of  Ex*  on 

r+\{x}. 

Theorem  2.1.  Recall  the  regularity  conditions  imposed  upon  Q+  and  the  inequalities  Re/c+  >0  and 
Im  K+  >  0  fulfilled  by  k+. 

(i)  The  linear  map  p  p,^^  is  a  bijection  oflV(Q+;K+)  onto  I.2(r+). 

(»i)  By  (i),  the  linear  operator  ■  L2(r+)  — ►  L2(r+)  given  by 

Ax+g  ■■=  <Pg\r^,  wherein  pg  G  W(n+;K+)  with  Pg,n+  =  9,  for  each  g  G  i2(r+) 

is  well-defined.  This  operator  is  compact,  injective,  and  has  dense  range  in  L2(r+),  while  is 

“strictly  dissipative,”  i.e., 

Im(«^/lK^/,/)£,(r+)  <  0  whenever  /  G  T2(r+)  and  /  5^  0.  (2.1) 

(jh)  For  any  g  G  L2{T+),  the  corresponding  unique  pg  G  W(n+;«+)  such  that  V’y>ii+  =  g,  i.e.,  the  corre¬ 
sponding  unique  solution  in  W(fi+; «+)  of  (1.23),  (1-34)  with  Neumann  data  g,  is  given  by 

=  \  {K*9- E^'tn^.AK+g}d>ir^  forx£fl+.  (2.2) 

Proof  :  This  is  proven  in  Ref. 4,  §6  and  §7;  cf.,  also,  the  remarks  there  in  §2  following  the  definition  in  (2.2)i. 
Actually,  in  Ref.  4  it  is  assumed  that  k+  ^  0,  with  ImK+  >  0,  and  «:+  >  0  if  Im/c+  =  0.  It  is  shown  there 
that  the  statements  (i),  (i«),  and  («i)  hold  on  that  larger  set  of  K:+-values,  with  the  exception  of  (2.1). 
Rather,  there  is  a  generalization  of  (2.1)  asserting  that,  for  any  such  the  operator  — is  strictly 
dissipative,  i.e., 

Im(CAK^./,/)i,(r+)  <  0  whenever  /  G  Z.2(r+)  and  /  /  0, 


16 


NRL  REPORT  9172 


for  any  C  lying  in  the  complex  set  defined  by 


:= 


J  {ceRlC>0}  =  (0,oo), 

\  {c  e  c  I  imC  >  0,  <  o, 


if  Im/c+  =  0, 


<  0,  and  +  [Im(C«:+^)]^  >  0  },  ifIm/c+>0. 


One  can  check  that,  for  the  larger  set  of  «+-values  specified,  lies  in  iff  Re/c+  >  0,  whence  (2.1) 
results.  ■ 


The  operator  Ak^.  of  Theorem  2.1  can  be  referred  to  as  the  (L2{T+)-)  boundary-data  operator  for 
VT(n+ ;«:+),  since  it  maps  the  Neumann  data  to  the  corresponding  Dirichlet  data  <p\r+  for  each 

ip  G  K+).  This  operator  has  a  number  of  interesting  properties  in  addition  to  those  just  cited  in 

Theorem  2.1,  some  of  which  are  given  in  Refs.  4,  5,  and  10.  Henceforth,  we  shall  suppose  that  we  have  at  our 
disposal  some  means  for  constructing  the  image  A^_^g  for  each  g  G  L2(r+),  which  is  essentially  equivalent 
to  assuming  that  we  can  construct  the  solution  of  the  Neumann  problem  for  (1.23)  and  (1.34)  in  0+  for  any 
boundary  data  in  L2(r+).  Various  schemes  are  available  for  accomplishing  this.  For  example,  an  explicit 
representation  of  the  operator  A^^  can  be  obtained  if  one  is  willing  to  carry  out  the  Gram-Schmidt  or¬ 
thonormalization  procedure  for  an  appropriately  chosen  complete  family  in  Z,2(r+);  cf.  Lemma  6.3  of  Ref.  4. 
If  is  not  in  the  countably  infinite  collection  of  (positive)  Dirichlet  eigenvalues  for  the  operator  —A  in 
the  interior  domain  Q_  U  flo  (=  \  f^+),  then  is  given  by  (I  -f  )  '  ■5*+ 1  wherein  and  Dk+ 

are  the  compact  operators  in  L2(r+)  constructed  from  the  “direct  values”  of,  respectively,  the  single-  and 
double-layer  potentials  based  upon  the  fundamental  solution  £*+:  for  /  G  L2(r+),  S*+/  and  are  the 

elements  of  12(1+)  defined  by  setting,  for  Ar+-a.a.  x  G  r+, 

S.^f(K)  :=  f  E^-^fdXr^ 

and 

However,  unless  Im/c^.  >  0,  one  does  not  usually  know  whether  k\  is  not  amongst  the  eigenvalues  just 
specified,  i.e.,  whether  the  operator  (/  +  D^^.)  is  injective,  and  so,  in  general,  a  more  difficult  operator 
inversion  must  be  effected  to  capture  ^1^+.  For  example,  it  is  shown  in  Theorem  8.4  of  Ref.  4  that  one  can 
always  compute  on  the  basis  of  the  representation 

-^(/ -£>.+  )}{/ + 

wherein  D*^  is  the  compact  operator  in  Z,2(r+)  obtained  by  conjugating  the  kernel  of  the  (L2(r4.)-adjoint) 
integral  operator  ^  is  any  conveniently  chosen  complex  number  with  positive  imaginary  part  (or,  in  the 
more  general  case  identified  in  the  proof  of  Theorem  2.1,  with  Im^  ^  0  and  Im^  ReK+  >  0);  and 
is  the  densely  defined  and  unbounded  operator  in  L2{T+)  induced  by  the  normal-L2  normal  derivative  of 
the  double-layer  potential  in  0+  (or  in  n_).  That  is,  defining  the  (restriction  to  G+  of  the)  double-layer 
potential  with  density  /  G  L2(r+)  by 

I  for  X  G  , 

then  WK^.,n+f  is  defined  to  be  the  normal-L2  normal  derivative  lV,^{/},n+,  on  the  linear  manifold  of  all 
/  G  L2(r+)  for  which  this  normal  derivative  exists  (which  turns  out  to  be  precisely  7i(vl*^),  the  range  of  the 
operator  A^^;  cf.  Ref. 4,  'Corollary  8.3).  Another  scheme  for  calculating  the  action  of  A^^  is  implicit  in  the 
results  of  Ref.  1 1  and  stated  more  explicitly  in  Theorem  8.1  of  Ref.  4.  Although  it  has  not  been  formulated  in 
this  manner  elsewhere,  the  problem  of  deriving  representations  for  that  are  valid  even  when  happens 
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to  be  one  of  the  Dirichlet  eigenvalues  of  —A  in  n_  Ufto  has  been  the  object  of  much  study;  cf.,  e.g.,  Refs.  11, 
12,  and  the  references  given  therein. 

By  relying  on  properties  of  the  single-  and  double-layer  potentials  that  are  proven  by  Kirsch  [13],  some 
connections  are  m.ade  in  Ref.  10  between  the  Z,2-type  results  just  cited  and  those  that  can  be  established 
within  the  Sobolev-space  setting.  For  example,  with  f2_  Uflo  of  type  as  we  are  supposing  here,  it  can  be 
shown  that  'JZ{Ak^)  lies  in  /f^(r+).  For  a  given  real  number  s,  with  more  smoothness  hypothesized  for  F^. 
(depending  upon  s),  it  turns  out  that  the  restriction  forms  a  bounded  bijection  of  //’(F+j  onto 

^j+i(r^)  when  s  >  0,  while  ^4,^^  can  be  extended  to  a  bounded  bijection  of  //’(F+)  onto  //’■*"^(F+)  if  s  <  0. 
By  using  results  of  this  sort,  we  find,  for  example,  if  fi.UQo  is  of  type  and  g  is  a  member  of  /f  2  (r+),  then 
the  function  (fg  G  lV(f24.;  k^)  defined  by  (2.2)  lies  in  and  so  possesses  both  a  Sobolev  trace  and  a 

Sobolev  normal  derivative  on  F+,  elements  of  and  respectively;  moreover,  these  coincide 

with  its  trace  and  normal  derivative  on  F^.  in  the  normal-L2  sense,  given  by  AK^.g  and  g,  respectively  (cf., 
also.  Proposition  2.2,  infra).  Presently,  we  shall  see  that  we  need  to  construct  the  solution  ipg  of  the  exterior 
Neumann  problem  (with  (1.34))  here  only  when  the  boundary  data  g  £  // 3(F+);  the  preceding  remarks  will 
then  supply  the  Sobolev-space  properties  of  that  solution  when  F+  has  the  additional  smoothness  noted. 

To  concisely  formulate  the  conditions  (1.35)  and  (1.36),  it  is  convenient  to  introduce  the  linear  space 
complex  amplitudes  of  strong  nonradiating  modes  for  Qo  and  w: 

:=  {u,^  e  1  (1.27)  through  (1.29)  hold).  (2.3) 

Here,  we  mean  that  (1.27)  through  (1.29)  are  to  hold  a.e.  on  the  respective  sets  indicated,  with  respect  to 
the  appropriate  measure  (A3  or  Ar„),  the  conditions  in  (1.28)  and  (1.29)  being  interpreted  in  the  sense  of 
the  Sobolev  trace  operator  on  H^;  cf.  (1.4)  through  (1.6)  for  the  calculation  of  the  action  of  T"°[  ]  on  an 
element  of  Now,  with  (A/n,  denoting  the  orthogonal  complement  of  w  taken  with  respect  to 

H°,  we  shall  write  ^  for  the  linear  manifold  of  elements  of  this  orthogonal  complement  that  also  lie  in 

(2.4) 

Since  is  closed  in  H°,  it  is  easy  to  see  that  ^  is  closed  in  H^.  Eventually,  we  shall  find  that 

finite-dimensional  (a  fact  that  is  not  surprising,  since,  when  it  is  nontrivial,  A/^  subspace 

of  an  eigenspace  for  the  traction  problem  for  in  f2<,),  and  so  also  closed  in  each  H*,  fc  =  0,2.  In 

particular,  once  this  is  known,  we  shall  have  the  orthogonal  direct-sum  decomposition 

(2.5) 

whence  it  shall  follow  readily  that 

H'=Afl.w  +  Ad?j,,w  with  =  {0},  (2.6) 

j.e.,  we  will  have  the  (nonorthogonal)  direct-sum  decomposition  of  into  the  subspaces  indicated.  In  this 
case,  we  shall  reduce  the  problem  corresponding  to  (1.23)  through  (1.26),  (1.34),  and  (1.36)  to  that  in  which 
(1.36)  is  replaced  by  (1.35),  essentially  by  seeking  separately  the  projections  of  the  solution  onto  the  subspace 
of  amplitudes  of  nonradiating  modes  and  onto  the  H°-orthogonal  complement  of  that  subspace  (and  we  shall 
effect  this  reduction  in  both  the  strong  and  the  weak  formulations).  Since  we  must  frequently  refer  to  these 
two  sets  of  requirements,  it  is  convenient  to  introduce  some  abbreviated  references. 

Terminology  2.1.  Let  G  H°  and  to  G  L2(Fo).  Let  fl‘  be  open  in  R®  and  contain  Q_  UOo,  and  suppose 
that  <p‘  G  C^(n‘)  satisfies  (1.18).  Then  by  (Po)  =  (Po(foitoi  <p‘))  we  refer  to  the  problem  of  determining 
an  appropriate  pair  of  functions  (u,^)  satisfying,  in  a  sense  that  is  always  to  be  specified  in  the  context, 
the  conditions  (1.23)  through  (1.26),  (1.34),  and  (1.35).  Suppose  also  that  TJ  on  Fo  x  (0,oo)  and  FJ  on 
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fio  X  (0,  oo)  are  measurable  C®-valued  functions  such  that  TJ(-,  t)  G  L2(ro)  and  FJ(-,  <)  G  H'’  for  each  <  >  0 
and 

=  T;(-,s)  •  u^|r„  dAr.  +  F^O,  s)  •  O,  dAgj  e*- ds  (2.7) 

exists  whenever  Ui^  G  satisfies  (1-27)  through  (1-29)  in  a  specific  sense  to  be  designated  in  the  context; 
for  example,  for  the  existence  of  the  integrals  in  (2.7),  it  is  sufficient  to  know  that  the  functions  t  i—> 
l|To(',t)|li„(r„)  and  t  ||FJ(-, f)ll„o  are  in  Li(0,oo).  Then,  by  (P)  =  (P(fo,  t„,  V?*;  FJ, TJ))  we  shall  refer 
to  the  problem  of  determining  an  appropriate  pair  of  functions  (u,  (p)  satisfying,  in  a  sense  to  be  designated 
in  the  context,  the  conditions  (1.23)  through  (1.26),  (1-34),  and  (1.36). 

Of  course,  it  is  not  really  necessary  to  retain  the  explicit  appearance  of  the  complex  amplitude  <p‘  in 
the  conditions  in  (1.25),  since  we  could  just  replace  the  inhomogeneities  there  with  some  given  elements 
of  L2(r+)  and  L2(r+),  respectively.  However,  we  wish  to  maintain  contact  with  the  underlying  physical 
setting,  and  so  choose  to  continue  writing  the  interface  conditions  in  their  original  form  (1.25). 


We  shall  first  define  “strong  solutions”  of  (Po)  and  of  (P),  since  one  has  for  such  solutions  a  better  feel 
for  the  sense  in  which  the  various  requirements  are  fulfilled.  Eventually,  however,  we  restrict  our  attention  in 
the  present  report  to  the  study  of  weak  formulations  of  these  problems.  Whatever  else  is  demanded  of  them, 
all  solutions  of  (1.24)  (or  (1-27))  that  we  consider  shall  be  H^-functions  u  that  ate,  in  particular,  solutions 
in  the  distributional  sense,  t.e.,  such  that 

/  u  ■  {Ax  dXs  =  -  f  fo  -vdAa  for  every  v  G  C§°(fto), 

’  Jcio 

whence  their  (local)  interior  regularity  properties  will  be  known  from  the  well-developed  theory  of  regularity 
for  strongly  elliptic  systems.  For  example.  Theorem  3.1  of  Ref.  14  implies  that,  if  fo  G  H*  for  some  integer 
A:  >  0,  then  any  u  G  H°  satisfying  (1-24)  in  the  distributional  sense  will  have  restrictions  in  for 

any  open  f2  with  f)  C  flo-  The  Sobolev  Imbedding  Theorems  then  show  that,  if  G  H^,  such  a  u  will 
be  in  C^(S7<,).  and  so  also  a  classical  solution  of  (1-24)  (in  fact,  the  same  conclusion  can  be  drawn  even  if 
it  is  known  only  that  f,  lies  in  C'’'*'(no));  if  f<,  G  C“’(Qo),  then  u  will  satisfy  that  inclusion  as  well  (for 
these  assertions,  cf.  Theorem  3.II  of  Ref.  14).  If  fo  has  real-analytic  real  and  imaginary  parts,  the  same  will 
be  true  of  u  (cf.  Ref.  15).  Thus,  the  elements  of  Wq  ^  enjoy  the  latter  property  of  real-analyticity  (as  do 
those  of  A/'^^  to  be  defined  presently).  All  of  our  solutions  of  (1-23)  will  be  classical,  in  W(n+;  K+),  and 
consequently  will  also  have  real-analytic  real  and  imaginary  parts.  Global  regularity  (“regularity  up  to  the 
boundary”)  for  the  solutions  of  our  boundary  and  interface  problems  is  another  matter,  that  we  shall  touch 
upon  in  this  report  only  in  Proposition  2.2,  in/ra,  and  there  we  shall  give  a  limited  result  only  for  the  scalar 
part  of  a  solution. 


Definition  2.1.  Recall  Terminology  2.1.  The  pair  (u,y>)  is  a  strong  solution  of  (Po(fo,  t<,, <p‘))  df  u  G 
•^n„  u)>  f  €  W(n+;  K+),  (1.24)  holds  a.e.  in  flo,  (1-25)  holds  a.e.  on  r+,  and  (1.26)  is  true  a.e.  on  r_,  with 
u|r^,  T"+[u],  and  T"-[u]  interpreted  in  the  sense  of  the  Sobolev  trace  operator,  and  and  y?,  n+  taken 
in  the  normal-L2  sense  on  r+.  Suppose  it  is  known  that  ^  is  finite-dimensional,  with  dimension  >  0; 
if  Uu)  >  1,  let  constitute  an  H°-orthonormal  basis  for  Suppose  that  the  integral  in  (2.7) 

exists  for  every  u^,  €  Then  the  pair  {u,ip)  is  a  strong  solution  of  (P(f<,,  to,  FJ,  TJ))  iff 


U  =  Uo  -1- 


0  if  Uu,  =  0 


and  (uo,<^)  is  a  strong  solution  of  (Po(fo,to,¥’'))  (ef-  (2-7)  for  the  coefficients  in  (2.8)  when  >  1). 


(2.8) 
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The  requirement  that  the  vector  part  of  a  strong  solution  of  (Pq)  lie  in  ^  is  the  strong  form  of  the 
condition  (1.35);  the  requirement  that  the  vector  part  of  a  strong  solution  of  (P)  be  of  the  form  (2.8)  with 
(in  particular)  uq  €  strong  form  of  condition  (1.36),  ensuring  that  the  latter  holds  for  every 

As  a  consequence  of  Definition  2.1,  we  essentially  need  consider  only  the  study  of  the  existence, 
uniqueness,  and  construction  of  strong  solutions  of  the  problems  (Pq).  We  will  find  that,  even  in  the  best 
of  circumstances,  one  must  identify  approximate  the  elastic-field  portion  of  a  solution  of  the  latter 

problem,  so  that  there  will  be  but  little  additional  work  involved  in  using  (2.8)  if  the  ultimate  goal  should 
in  fact  be  the  solution  of  a  problem  (P). 


Obviously,  under  Definition2.1  (P(fo,  t^,  y>‘;  F^,  TJ))  can  have  at  most  one  strong  solution  if  the  problem 
(Po(fo,to,y>‘))  enjoys  that  same  property.  To  emphasize  the  role  played  by  (1.34)  and  (1.35),  let  us  prove 
that  there  can  exist  at  most  one  strong  solution  of  the  boundary  and  interface  problem  (Po).  In  preparation, 
we  need  a  familiar  integral  relation  (which  shall  be  used  later,  as  well),  viz., 

I  {v  - At(  ^u-l-fx,;i(u,v)}clA3=  /  vlr„  •T“»[u]dAr„  for  u  6  and  v  e  (2.9) 

in  which 

^^A,;i(u,v)  :=  A£jj[u]eti[v]  +  2^£:ji[u]£jt[v]  =  <rjfc[u]£jj:[v]  =  t^jit[v]£jt[u]  for  u,  v  G  H'.  (2.10) 

This  identity  is  easily  derived  with  an  application  of  the  formula  for  integration  by  parts,  cited  previously. 

Theorem  2.2.  There  can  exist  at  most  one  strong  solution  of  the  problem  (Po(f<,,  t^,  i^*)),  in  the  sense  of 
Definition  2.1. 

Proof  -.  Let  the  pair  (u,  ^)  be  a  strong  solution  of  (Po(0,0,0)),  t.e.,  with  fo  =  0,  t^  =  0,  and  =  0.  Then, 
by  using  (1.25)  and  the  mapping  property  of  following  from  its  definition  in  Theorem  2.1(m),  and  noting 
that  T"+[u]  is  a  normal  field  on  r+,  we  have  the  equalities 


=  (-»f>+wu|r+  -n^.)  ^--^T"+[u]  •  n+^ 


tg+u 


T“+[u]-u|r+ 


7+  -  jw 

between  elements  of  Li(r+),  whence  a  little  rearrangement  and  integration  over  r+  produces 

f 

)t,(r+)  = -e+  /  T"+ [u]  ■  ulr+ dAr+ . 

"  ./r+ 

On  the  other  hand,  from  (2.9)  and  by  taking  into  account  (126)  and  (124),  we  get 

/  T"+  [u]  u|r^  dXr^  =  f  T"”  [u]  •  u|r„  dAr, 

Jr„ 

=  /  {u  •  -|-fA,^(u,u)}dA3  =  /  {-£ioW^u  •  u  +  ^A,/i(u,  u)} 


(2.11) 


dAa 


clearly  implying  that  the  integral  on  the  right  in  (2.11)  is  real.  Thus,  («:^^K+V?,n+ > ‘^m+)i2(r+)  i®  ®o 
that  the  vanishing  of  follows  from  (2.1).  Statement  (i)  of  Theorem  2.1  then  allows  us  to  eissert  that  (p 
vanishes  in  fl^..  Now,  (1.24)  through  (126)  and  the  inclusion  u  G  show  that  u  G  since  we  also 

know  that  u  G  (A^n„  Definition  2.1),  it  must  be  that  u  =  0.  This  completes  the  proof.  I 

It  is  interesting  to  note  that,  without  the  condition  that  u  lie  in  we  could  still  have  deduced 

that  =  0  by  reasoning  as  in  the  proof,  but  then  would  have  been  able  to  conclude  only  that  u  G  -Vq  ui- 


20 


NRL  REPORT  9172 


Now,  we  turn  to  the  reduction  of  the  strongly  formulated  boundary  and  interface  problem  (Pq)  to  a 
purely  interior  problem  (in  flo)  with  a  nonlocal  boundary  condition,  basing  the  replacement  on  the  properties 
of  the  operator  Suppose  that  (u,<p)  is  a  strong  solution  (and  so  also  the  unique  strong  solution)  of 

(Po).  Then,  by  the  boundary-data-mapping  property  of  and  the  first  equality  in  (1.25),  we  have 


<p\r+  =  =  ^<c+(»^+wu|r+  •  n+)  - 

and  so  the  second  equality  in  (1-25)  becomes 


T"+[u]4-(7+  -»w){ip+a/^*^(u|r+  ■  n+)  -  =  t+  -(7+  -iw)^‘|r+n+. 


or 

T"+[u]  +  p+c5.»c^{>lK+(u(r+  •n4.)}n+  =  t+  -  (7+  -  iu;){<p‘|r+  -  }n+,  (2.12) 

a  boundary  condition  on  u  alone.  This  motivates 

Definition  2.2.  Recall  Terminology  2.1.  The  function  u  is  a  strong  solution  of  the  purely  interior  problem 
corresponding  to  (Po(f£)i  ta,  (p‘))  iff  u  G  (1.24)  holds  a.e.  in  Q^,,  (2.12)  is  true  a.e.  on  r+,  and  (1.26) 

holds  a.e.  on  r_. 

To  find  a  strong  solution  of  the  original  problem  (Pq),  it  is  easy  to  see  that  it  suffices  to  generate  a 
strong  solution  of  the  corresponding  purely  interior  problem.  For,  if  u  fulfills  all  requirements  of  Definition 
2.2,  let  us  construct  the  unique  <p  G  lT(f2+;«+)  having  the  Neumann  data  given  by 

=  *?+‘*'u|r+  •  n+  -  (2.13) 

observe  that  the  function  on  the  right  here  lies  in  ff  3(r+)  (cf.  our  previous  remark).  Explicitly,  from  the 
representation  (2.2)  of  Theorem  2.1,  this  is  given  by 

¥>(x)  =  -  I  J 

^  dr+ 

+  {■£’x‘^u|r+ -£'*t„^>l«+(u(r+ •n+)}dAr+  for  each  xGn+.  (2.14) 

With  (2.13),  we  have  automatically  ensured  that  the  first  equality  in  (1.25)  is  fulfilled  a.e.  on  r+,  while  the 
validity  of  the  second  requirement  there  follows  in  the  same  sense  from  (2.13)  and  (2.12): 

T"+[u]  +  (7+  -  iw)<p|r+n+  =  T"+[u]  +  (7+  -  iw)(ylK^(p,„^)n+ 


=  -p+w(w  +  i7+){^,^(ulr+  ■n+)}n4.  +  t+  -  (7+  -  iw){v?‘lr+  -  ^K+<p‘n+}n+ 

+(7+  -  *w){ie+u;^;5^(ulr+  •  n+)  - 

=  t+  -  (7+  -  iw)<p‘|r+n+. 

Consequently,  (u,  ip)  satisfies  all  of  the  requirements  placed  in  Definition  2.1  on  a  strong  solution  of  (Po).  It  is 
important  to  note  here  that  only  the  normal  component  u|r^  11+  of  the  trace  of  u  restricted  to  r+  is  needed 
for  the  construction  of  ip  by  use  of  (2.14),  i.e.,  to  calculate  the  fluid  field  once  the  strong  solution  of  the 
purely  interior  problem  for  the  elastic  field  has  been  obtained;  the  complete  implication  of  this  observation 
for  the  numerical  work  will  emerge  in  §4.  We  might  also  remark  that  the  representation  displayed  for  p  in 
(2.14)  has  split  naturally  into  a  sum  of  two  integrals  taken  over  r+.  Of  these,  Theorem  2.1  makes  it  clear 
that  the  first  is  a  representation  for  the  unique  element  of  W(n+;  /c+)  with  Neumann  data  which  is 


21 


ALLAN  G.  DALLAS 


nothing  else  but  the  complex  amplitude  of  the  time-harmonic  scattered-field  velocity  potential  that  would  be 
produced  by  the  interaction  of  the  time-harmonic  incident-field  velocity  potential  with  complex  amplitude 
with  a  rigid  obstacle  occupying  n_  Ufioi  under  the  usual  steady-state  assumption.  The  second  term  on  the 
right  in  (2.14)  must  then  encompass  all  contributions  from  the  elastic  nature  of  the  obstacle,  the  impressed 
tractions,  and  the  imposed  body  force.  By  using  Green’s  Theorem,  the  first  term  can  be  recast  into  the 
form  —(1/2)  ■E'xtn+{¥’‘lr+  —  dAr^,  in  which  there  appears  the  same  field  v^‘|r+  -  "^<5+ V’>n+  on 

r+  as  in  the  new  boundary  condition  (2.12)  of  the  purely  interior  problem;  this  is  the  sum  of  the  incident 
and  the  scattered  fields,  or  “total  field,”  on  r+  in  the  rigid-acoustic  steady-state  scattering  interaction  just 
mentioned. 


It  is  appropriate  to  include  here  an  observation  pertinent  to  the  eventual  numerical  approximation  of 
^  and  its  far-field  pattern  from  numerical  approximations  generated  for  u  by  computations  based  on  the 
interior  problem.  We  recall  that  to  each  element  ii>  of  lT(n+;K.).)  there  is  cissociated  a  unique  complex 
function  V^oo  on  the  unit  sphere  Si  of  with  the  property  that 


i/;(re) 


giic  +  r 

- V'oo(e)  -f-  O 

r 


as  r  —*  oo,  for  each  e  £  Si . 


This  function  i/jqo  is  variously  termed  the  far-field  pattern,  or  radiation  pattern,  of  rp-,  cf.,  e.g..  Ref.  12.  Since 
V'oo(oj  =  /r)~^ip(re)  for  each  e  £  Si,  by  using  the  representation  for  elements  of  lT(fl+;/c+) 

that  is  obtained  from  (2.2),  it  is  easy  to  show  that 


0oo(e)  = ^  {e^+V',n+ foreache£Si,  (2.15) 

whenever  £  W(f2+;/c+),  with  e^‘''(y)  :=  for  y  £  R^.  Now,  we  can  formulate  the  desired 

statement  in  the  form  of 

Proposition  2.1.  Let  u  be  an  element  vfH‘,  and  let  ^  be  the  unique  element  of  lV(n+,  k+)  satisfying 
(2.13),  so  that  ‘f  has  the  representation  displayed  in  (2.14);  of  course,  if  u  happens  to  be  a  strong  solution  of 
the  purely  interior  problem  corresponding  to  (Pq).  then  (u,v?)  is  the  unique  strong  solution  of  the  problem 
(Po).  Suppose  that  (Un)^_;,r  is  a  sequence  in  H'  converging  weakly  to  u  in  H' .  Construct  the  corresponding 
sequence  (v’nl^rr/v  l^(^  +  t'«+)  according  to,  for  n>  N , 


<fn(x) 


(  {£:;;+u„|r+ ■  n+ -  >4,t^.(u„|r+ ■  n+)}  dAr+  for  each  x  £  (2.16) 


Then,  for  each  3-index  o.  the  sequence  (ip„'a)^=/v  converges  to  y?,„  uniformly  on  each  subset  ofn.^.  that  is 
closed  in  R'^.  Further,  the  sequence  (‘^noo)^=,v  of  far-held  patterns  converges  to  the  far-held  pattern  of 
ip  uniformly  on  Si . 

Proof  :  We  have,  for  each  x  £  Q+  and  n>  N, 

p{x)  -  pu{x)  =  [  {Ex+(u  -  u„)ir+ ■  n+ -  -  u„)|r.^  ■  11+)}  dAr+.  (2.17) 

Jr+ 


Now,  choose  any  set  F  that  is  closed  in  R^  and  contained  in  0+.  Since  ImK+  >  0  and  the  distance  between 
F  and  the  compact  set  r+  :=  50+  is  positive,  it  is  clear  that  E^‘''|r+  and  are  bounded  on  1'+. 

uniformly  for  x  £  F,  Thus,  by  applying  the  Cauchy-Schwarz  inequality  to  (2.17)  and  using  the  boundedness 
of  in  h2(f+).  WT  find  that 

sup  |v5(x)  -  i^„(x)|  <  ^^Ar+(r+){c,  +  c2||.4„,^|(}|1(u  -  u„)|r+|lLj(r+)  for  ”  >  E, 
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wherein  ci  and  C2  depend  upon  only  F  and  k^.  The  compactness  of  the  Sobolev  restricted-trace  map 
V  from  into  L2(r+)  and  the  weak  convergence  of  (u„)“_jy  to  u  in  then  imply  that  the  first 

assertion  of  the  Proposition  is  true  for  a  =  (0,0,0).  The  proof  for  any  other  3-index  a  can  obviously  be 
carried  through  in  an  entirely  analogous  manner,  since  the  corresponding  partial  derivatives  can  be  computed 
by  differentiation  under  the  integral  in  (2.14)  and  (2.16).  The  final  assertion,  concerning  the  convergence  of 
the  far-field  patterns,  is  established  just  ais  eeisily,  by  similar  reasoning  based  on  (2.15).  ■ 

In  the  applications  of  Proposition  2.1  that  are  pertinent  to  the  present  setting,  we  shall  even  have 
available  the  norm-convergence  of  (u„)^j^  to  u  in  H’  (cf.  Theorem  4.1,  infra),  not  merely  weak  convergence; 
assuming  some  additional  regularity  for  r+,  we  can  sketch  the  proof  of  a  correspondingly  stronger  result  for 
the  convergence  of  the  sequence  (v?n)“=A/ 

Proposition  2.2.  Retaining  the  setting  and  notation  of  Proposition  2.1,  suppose  now  also  that  fi-  Uf2o  is 
of  type  while  (  converges  to  u  in  the  norm  ofH^.  Then  ip  and  the  sequence  (<Pn)^-fv  he  in 

^ioc(^+).  '^hh  lirrv,  — *00  pn  —  F  i”  the  locally  convex  topology  of  that  space.  Moreover,  ip  and  the  sequence 
he  in  C°'2(fl+),  and  lim„_oo(^„  =  p  in  that  Banach  space  (so,  in  particular,  the  convergence  is 
uniform  on  Q+  ). 

Proof  :  Now,  we  know,  by  Lemma  2.1,  that  (u„lp^  ■»+)"=jv  converges  to  u|p^  ■  n+  in  the  norm  of  //§(r+). 
We  recall  the  definition  of  the  double-layer  potential-operator  (restricted  to  Q+)  f  i— ►  W^^{f}  on  L2(r+) 
(supra),  and  now  define  the  single-layer  potential-operator  (restricted  to  f2+)  /  on  L2(r+)  by 

V^4{/}W-=/  Ei*fdXr^  forxen+. 

We  have  sufficient  regularity  of  r+  to  ensure  that,  according  to  Ref.  10,  AK,^|//s(r+)  maps  H^{T+)  bound- 
edly  into  //» (r+),  and,  by  Theorem  2.16  of  Ref.  13,  that  the  double-layer  potential-operator  is  bounded  from 
//3(r+)  into  H/l^(Q+),  while  the  single-layer  potential-operator  is  bounded  from  //a(r4.)  into  Hjl^(Q+). 
Combining  all  of  these  facts,  in  view  of  (2.14)  and  (2.16)  it  follows  that  p  and  (i^n)^=Ar  are  contained  in 
H(g^(Q+),  with  lim„_oo  V’n  =  ¥’>!'  fhat  space  (cf.  (2.17)).  Directly,  the  Sobolev  Imbedding  Theorem  (The¬ 
orem  5.4,  Part  II,  of  Ref.  9)  allows  us  to  assert  that,  for  any  bounded  open  set  Q  C  0+,  p\^  and  the 
are  elements  of  C°'3(n),  with  the  convergence  lim„_oo¥’n|D  =  'P\^  the  norm  of  the  latter  Banach  space. 
Taking  D  to  be  the  intersection  of  with,  say,  a  ball  of  sufficiently  large  radius,  and  accounting  for  the 
results  of  Proposition  2.1  concerning  the  uniform  convergence  of  the  pn  and  their  partial  derivatives  to  p  and 
its  partial  derivatives  in  the  exterior  of  the  ball,  we  see  that  the  pn  and  p  lie  in  the  Banach  space  €''^■^(0+), 
and  lim„^oo  Pn  =  p  in  the  norm  of  that  space.  ■ 


Returning  to  the  general  setting,  we  shall  study  the  purely  interior  problem  by  setting  up  and  analyzing 
a  weak  version  of  it;  this  will  also  lead  back  in  a  natural  manner  to  weak  formulations  for  the  original 
problems  (Po)  and  (P).  Regularity  results  for  solutions  of  the  weakly  formulated  interior  problem  would 
then  allow  us  subsequently  to  return  to  the  strong  setting.  As  usual,  the  weak  formulati,  n  arises  from  (2.9), 
written  in  the  form 


■f  -b  fo}  •  vdAs I  {£:a,^(u,  v)  -  •  v)  dAa 

Jiio  Jn„ 

=  /  T""[u]  v|r„(iAr„  +  /  fo  -vdAg  for  u  e  and  v  6  H'  (2.18) 
Jr„  Jcio 
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(with  f(,  e  H°).  Now,  if  u  is  a  strong  solution  of  the  purely  interior  problem  for  (Po)>  as  in  Definition  2.2, 
then  we  must  have 


/  QoUi^\i  'v}dX3  + Q+c^.K^  (>4K^(u|r+ ■  n+))v|i  ^  •  n+ cfAr+ 

JOo  Jr+ 

~  Jr  t-  v|r_dAr_ 

+  I  fo  'VrfAs  whenever  V  G  H' .  (2.19) 

•/n. 

This  leads  us  to  define  the  sesquilinear  form  (Tu,(-,  ■)  on  according  to 

<^w(u,v).-  /  {^a,m(u.  v)  -  •  v}  dAa  +  e+c^/c^  /  (/l;c+(u|r+ •  n+))v|r^.  •  n+ ciAr^., 

foru,vGH',  (2.20) 


and,  corresponding  to  selected  functions  fo,  to,  and  yj‘  as  in  Terminology  2.1,  the  conjugate-linear  functional 
Ao;  on  H’  by 

AwV  :=  ^  {t+ -  (7+ -  iu;){^‘|r+ -  •  v|r+dAr+ -I- ^  t_  v|r_dAr_ 

+  1  fo  -TdAa,  for  each  vGH’.  (2.21) 
dfio 

Then  (2.19)  is  written  concisely  as  <ru.(u,v)  =  A„v  for  v  €  H‘;  the  weak  formulation  of  the  purely  interior 
problem  corresponding  to  (Pq)  shall  be  set  up  on  the  basis  of  this  equality. 


It  is  also  convenient  to  identify  a  linear  manifold  comprising  complex  amplitudes  of  weal-  nonradiating 
modes,  proceeding  from  (2.18).  If  Uo,  G  Mq  then  (2.18)  shows  that 


-  ■  v)  <1X3  =  0 


for  every  v  G  H' , 


(2.22) 


while  also  lies  in  the  subspace  Ho„^  given  by 


:={uGHMu|r+n+  =0) 


(2.23) 


(the  closure  in  H'  of  the  collection  of  all  u  G  C'(0<,)  such  that  u|r+  •  n^.  =  0).  With  this,  we  define  the 
linear  manifold  A/^^  in  H'  by  setting 

:=  {  G  I  (2.22)  holds  } .  (2.24) 

Then  C  A/^  At  least  under  sufficiently  strong  hypotheses  of  smoothness  for  To,  we  shall  also  have 

a^A„,.  c  whence  it  shall  follow  that  Wn„  u,  =  A^n„w  example,  in  the  extreme  case  in  which  Do 
is  of  type  we  shall  have  A/p^  ^  C  C°°(n<,);  cf.  Ref.  14.  Proceeding  as  we  did  for  "’e  denote  by 

(A/n„  ^he  orthogonal  complement  of  taken  with  respect  to  H“,  and  introduce  ^  C  H'  by 

n  (2.25) 

u  ,  is  a  closed  linear  manifold  (subspace)  in  H*.  We  shall  show  later  that  AA^  ^,  finite-dimensional, 
and  so  also  closed  in  both  H°  and  H’,  Once  this  is  known,  just  as  before  (cf.  (2.5)  and  (2.6))  it  will  follow 
that  the  direct-sum  decompositions 

H«  =  .VA„,,„©(.VA„,j'-^^<’  (2.26) 
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and 

H' with  =  (2.27) 

are  valid  (of  which  the  first  is  orthogonal).  We  still  have  the  interior  regularity  result  asserting  that  the 
elements  of  „  have  real  and  imaginary  parts  that  are  real-analytic  in  fio,  so  (1.27)  holds  in  the  classical 
sense  in  Qo  for  each  114^,  e  „ . 

As  already  tacitly  noted  in  Proposition  2.1,  for  the  construction  of  the  unique  tp  €  k+)  satisfying 

(2.13),  we  need  know  only  that  u  6  H^.  This  affords  a  simple  means  for  providing  weak  formulations  for  both 
(Po)  and  (P),  based  upon  that  for  the  purely  interior  problem  for  (Po).  Precisely,  the  weak  formulations 
are  given  in 

Definition  2.3.  Recall  Terminology  2.1.  The  function  u  is  a  weak  solution  of  the  purely  interior  problem 
corresponding  to  (Po(fo,  t<,,  ip‘))  iff  u  G  Ad^^  ^  and 

v)  =  Ao/V  for  every  V  €  (2.28) 

with  the  sesquilinear  form  tr^  and  the  conjugate-linear  functional  A^  defined  by  (2.20)  and  (2.21),  respec¬ 
tively.  If  u  is  a  weak  solution  of  the  purely  interior  problem  corresponding  to  (Po(fD,  t^,  yj‘))  and  p  is  the 
corresponding  unique  element  of  1^(04.; /c+)  such  that  (2.13)  holds,  then  the  pair  (u,ip)  is  termed  a  weak 
solution  of  (Pn(fo,  to,  <p‘)) .  Suppose  further  that  the  integral  in  (2.7)  exists  for  every  Uu;  G  A^^^,  and  it 
is  known  that  .Vq  ^  has  finite  dimension  >  0;  if  >  1,  let  be  an  H°-orthonormal  basis  for 

Then  (UjV?)  is  termed  a  weak  solution  of  (P(fo,  to,  FJ,  TJ))  iff  u  can  be  written  as  in  (2.8)  and 
{uo,<p)  is  a  weak  solution  of  (Po(foi  to,  v?‘))- 

Observe  here  that,  for  (Po),  the  requirement  u  G  AdJj^  constitutes  a  weak  form  of  (135),  while,  for 
(P),  the  condition  that  u  have  the  form  in  (2.8)  provides  a  weak  realization  of  (1.36). 

According  to  Terminology  2.1,  we  have  been  supposing  that  the  piece  of  data  to  is  in  L2(ro);  one 
can  relax  this  requirement  by  positing  instead  the  inclusion  to  G  H“>(ro),  provided  that  appropriate 
modification  is  made  in  the  definition  of  A^.  in  (2.21),  viz.,  that  the  sum  of  the  terms  involving  t+  and  t_  be 
replaced  by  the  more  general  expression  (to,v(p^)i,  involving  the  duality  pairing  for  H3(ro).  The  resultant 
conjugate-linear  functional  A^  is  again  bounded  on  and  the  analysis  of  §4  can  be  carried  through,  mutatis 
mutandis.  For  simplicity,  we  shall  continue  to  suppose  that  to  G  L2(ro). 

It  is  clear  that  a  weak  solution  of  the  purely  interior  problem  corresponding  to  (Po)  that  is  also  known  to 
possess  some  additional  regularity  will  be  a  strong  solution  of  the  purely  interior  problem  for  (Po).  For,  ifu 
satisfies  the  requirements  placed  on  a  weak  solution  of  the  purely  interior  problem  for  (Po)  in  Definition  2.3 
and  also  lies  in  then  u  G  A4^  ^  (since  the  inclusion  A/^  u/  C  A/^  implies  that  (A/^  ^„)^'’'^°  C 
(A/q^  while  (2.18)  and  (2.28)  (i.e.,  (2.19))  together  imply  that 

-/  vdAg-h  /  T"°[u]  vlr.dAr, 

Juo 

-  /  {-e+4'«+(^«+(«lr+  ■“+))"+  +t+  -  (t+  -  *w){¥5‘|r+  -  }n+}  •  v|r^dAr+ 

+  J  t_  •  v|r_  dAr_  for  every  v  G  H*.  (2.29) 

Since  the  equality  in  (2.29)  must  hold,  in  particular,  for  every  v  G  Co°(no),  (1-24)  follows  (a.e.  in  fio);  with 
the  simplified  form  then  taken  on  by  (2.29),  we  conclude  that  (2.12)  and  (126)  are  true  (a.e.  on  F^.  and  r_, 
respectively).  Thus,  such  a  u  has  all  of  the  properties  demanded  of  a  strong  solution  of  the  purely  interior 
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problem  in  Definition  2.2,  and,  moreover,  leads  directly  to  the  strong  solution  for  (Po)  itself  (as  described 
in  the  remarks  following  Definition  2.2)  and  to  the  strong  solution  for  (P)  (as  in  Definition  2.1).  Observe 
that  we  do  not  need  to  know  here  whether  ^  and  coincide. 

3.  VARIATIONAL  PROBLEMS  INVOLVING  A  DEFINITE  FORM  WITH  A  COMPACT 
PERTURBATION;  THE  GALERKIN  METHOD 

For  the  analysis  of  the  weak  formulation  of  the  purely  interior  problem,  we  require  some  facts  about 
problems  posed  “in  variational  form;”  in  particular,  for  the  construction  of  a  convergent  approximation 
scheme,  we  shall  use  results  of  Hildebrandt  and  Wienholtz  [6]  concerning  the  Galerkin  method. 

Throughout  this  section,  (H,(-,-)ff)  denotes  a  separable  complex  Hilbert  space,  cr(-,  )  a  bounded 
sesquilinear  form  on  H,  and  A  a  bounded  conjugate-linear  functional  on  H .  We  are  interested  in  studying 
the  solvability  of  the  problem 

find  f  £  H  satisfying  er(f,h)  —  Ah  for  all  h  £  H ,  (31) 


and  in  establishing  a  viable  and  convergent  scheme  for  the  approximation  of  a  solution  (when  at  least  one 
exists),  under  certain  hypotheses  on  the  form  a.  Of  course,  we  know  that  there  are  uniquely  determined  a 
bounded  linear  operator  To  on  H  and  an  element  of  H  such  that 


er(g,h)  =  {Tog,h)„ 

for  all  g,h  £  H 

(3.2) 

and 

=  (?A>  ^)h 

for  all  h£H, 

(3.3) 

so  that  (3.1)  is  equivalent  to  the  operator  problem 

find  f  £  H  satisfying  T^f  =  g/^^. 

(3.4) 

The  latter  equivalence  is  useful  for  examining  the  solvability  of  (3.1),  but  we  prefer  not  to  base  the  com¬ 
putations  of  a  numerical  scheme  upon  the  availability  of  explicit  analytical  forms  for  T„  and  since  such 
formulae  may  be  difficult  to  deduce. 

By  Af(L)  we  mean  the  null  space  of  the  operator  or  functional  L.  Af\{cr)  is  the  subspace  of  H  defined 
by 

Afi(cr)  {g  £  H  \  a(g,h)  =  0  for  all  h  G // } 

(throughout,  a  subspace  of  a  Hilbert  space  is  defined  to  be  a  linear  manifold  that  is  closed  in  that  Hilbert 
space).  The  (bounded)  sesquilinear  form  cr*(-,-)  adjoint  to  a  is  defined  by 

'■=  [or  g,h£  II \ 

tj’  is  represented  by  T*,  the  (Hilbert-space)  adjoint  of  To,  in  the  same  way  that  cr  is  represented  by  To- 
Corresponding  to  the  adjoint  form,  we  have 

M(o-*):={<?€//|<7*(!7,h)  =  0  forallh€//}  =  {<lG//  |<7(h,(/)  =  0  for  all  h  G // }. 

With  (3.2),  it  is  easy  to  check  that  A/’i(<r)  =  Af(To)  and  A^i(<t*)  =  Af{T*). 

Invariably,  we  assume  that  (t(-,  •)  is  of  the  form  6(-,  ■)  +  k(-,  •),  in  which  the  sesquilinear  forms  S  and  k 
are  bounded  on  II,  with  6  definite  and  «  compact.  By  the  requirement  of  definiteness  for  S,  we  mean  that 
there  exists  e<  >  0  such  that 


l'5(if.3)l  >  CillPlIw  for  all  3  G//; 


(3.5) 
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by  the  requirement  of  compactness  for  k,  we  demand  that  whenever  g  and  h  are  elements  of  H  and  (<in)“=o 
and  (hn)“=o  sequences  in  H  converging  weakly  to  g  and  h,  respectively,  then  {K{g„,hn))^_g  converges 
in  C  to  K{g,h).  Of  course,  a  compact  form  is  necessarily  bounded.  For  the  uniquely  determined  bounded 
linear  operators  7i  and  Tk  on  H  such  that 

^(l7. h)  =  h)„  and  K{g, h)  =  {T^g,  h)u  for  all  g,he  H,  (3.6) 

we  know  then  that  Ts  is  definite  {\(^tg,g)H\  ^  for  9  ^  and  Tk  is  compact  (since  ||T,/|(^  = 

for  /  €  /f  and  the  boundedness  of  certainly  implies  that  it  maps  weakly  convergent  sequences 
into  weakly  convergent  sequences).  Under  these  hypotheses,  the  solvability  of  the  problem  (3.1)  is  completely 
characterized  by  the  Fredholm  Theorems  that  are  available  for  it  in  the  form  of  the  following  statement  (for 
clarity,  it  is  desirable  to  provide  a  proof  of  these  facts  showing  that  their  validity  depends  on  only  the 
form  hypothesized  for  a,  and  is  independent  of  the  considerations  involved  in  establishing  an  approximation 
scheme;  cf.  Ref.  6). 

Theorem  3.1.  Let  the  sesquilinear  form  <t  on  H  be  the  sum  6  +  k  of  bounded  sesquiiinear  forms,  with  6 
definite  and  k  compact. 

(i)  The  dimensions  of  Afi(a)  and  are  finite  and  equal. 

(ii)  If  A  is  a  bounded  conjugate-linear  functional  on  H,  then  there  exists  some  fo  £  H  such  that 

<r(/o,  h)  =  AA  for  all  h  &  H  (3.7) 

iff  the  (“solvability”)  condition 

^rl(<r•)C^r(A)  (3.8) 

obtains,  in  which  case  the  collection  of  all  solutions  of  (3.1)  is  given  by  /o  +  A/]('r).  Consequently,  there 
exists  a  unique  solution  of  (3.1)  iffA/’i(<T)  =  {0}  (or,  equivalently,  iff  A/i(<t*)  =  {0}). 

(iii)  Let  M  be  a  subspace  of  H  that  is  complementary  to  Afi(«r),  i.e.,  such  that 

Ari(<T)nAd  =  {0}  and  H  ~Afi(<T)  +  M.  (3.9) 

Then,  whenever  A  is  a  bounded  conjugate-linear  functional  on  H  such  that  (3.8)  holds,  there  exists 
precisely  one  corresponding  £  Ad  for  which 


,h)  =  Ah  for  all  he  H\ 
moreover,  there  exists  M  >  0  with  the  property  that,  for  all  such  A 

WUmWh  <  m\\a\\ 


(3.10) 


(3.11) 


(with  ||A||  denoting  the  norm  of  A). 

Proof  '.  Adhering  to  the  notation  already  introduced,  we  have  T),  =  Tf  -1-Tk.,  with  Ts  definite  and  compact. 
It  is  easy  to  see  that  the  bounded  operator  Ts  is  bijective.  Indeed,  the  inequality  1{75<7,^)^|  >  q||^||J,, 
holding  for  all  g  £  H,  gives  \\Tsg\\H  >  cs\\g\\n  for  all  g,  whence  we  conclude  that  Ts  is  injective  and  has 
closed  range  TliTs)  (as  well  as  bounded  inverse  on  TUfTs)  into  H).  Clearly,  the  same  reasoning  can  be  effected 
for  Tf.  Thus,  TliTs)  =  A/'(T;)-‘-  =  H  (and  HiTf)  =  Af(Ti)-‘-  =  H),  proving  the  claim.  Now  we  can  write, 
e.g.,  Ta  =  Ts(I  ■\-Tf^TK.)  and  argue  in  the  well-known  manner,  from  the  validity  of  the  Fredholm  Theorems 
for  a  perturbation  of  the  identity  by  a  compact  operator,  to  conclude  that  those  same  statements  are  valid 
for  T„.  Having  secured  this,  all  of  the  assertions  of  the  theorem  follow  from  (3.2)  and  the  equivalence  of 
(3.1)  and  (3.4).  Thus,  (j)  holds,  in  view  of  the  equalities  A^i(<t)  =  Af{Ta)  and  A/i(it*)  =  M{Tf).  Further, 
there  exists  fo  £  H  such  that  (3.7)  is  true  iff  Tafo  =  9\  (cf-  (3-3)),  iff  9\  ^  —  A/i(«t*)-‘',  iff 
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A/»  =  {9\,h)fj  =  0  for  all  h  G  iff  (3.8)  holds.  This  proves  (li),  the  final  assertions  concerning  the 

collection  of  all  solutions  of  (3.1)  and  the  existence  of  a  unique  solution  now  being  obviously  true.  Finally, 
for  (»»),  let  fo  €  H  satisfy  (3.7),  wherein  A  is  a  bounded  conjugate-linear  functional  on  H  satisfying  (3.8) 
(such  an  /o  existing  by  (»)).  Writing  /o  =  /i  -t-  ^  M(<’’)  and  €  M,  then  (3.10)  clearly 

holds,  and  is  the  unique  element  of  M  satisfying  (3.10)  since  M  nA/’i(<r)  =  {0}.  This  reasoning  implies 
that  the  bounded  operator  T<r(wM  (the  restriction  of  Ta  to  M)  maps  the  subspace  M  bijectively  onto  the 
subspace  .Afi ((t* )-'■  =  AA(T*)-*-,  and  so  has  bounded  inverse  defined  on  the  latter.  With  M  denoting  the  norm 
of  this  inverse,  since  inequality  (3.11)  follows  upon  noting  that  i 

For  the  construction  of  approximations  to  solutions  of  the  problem  (3.1),  wc  first  consider  the  sp'^cial 
case  in  which  M\{<t)  is  the  trivial  subspace,  then  treat  the  general  case  in  which  may  be  nontrivial, 

under  an  additional  hypothesis  on  <r  (that  will  in  fact  obtain  in  our  subsequent  application  of  the  Galerkin 
method). 

Theorem  3.2.  Let  the  sesquiiinear  form  c  on  H  be  the  sum  S  +  k  of  bounded  sesquilinear  forms,  with  6 
definite  and  k  compact.  Suppose  further  that  Mi(cr)  =  {0}.  Let  A  be  a  bounded  conjugate-linear  functional 
on  H ,  and  denote  the  associated  unique  solution  of  (3.1)  by  f^.  Then  t>e  convergently  approximated 

in  the  norm  of  H  by  the  Galerkin  method  for  (3. 1 ),  by  using  as  coordinate  functions  any  family  that  is  linearly 
independent  and  complete  in  H .  That  is,  suppose  that  is  linearly  independent  and  complete  in  H . 

Then  there  exists  a  positive  integer  N  such  that  for  n>  N  the  system 

n 

53<^(hfc,/i,)Cfc  =  Ah,,  /=l,...,n,  (3.12) 

fc=i 

possesses  a  unique  solution  (4fc"^)2_i  and  the  resultant  sequence  (/„  :=  converges  to 

in  the  norm  of  H .  In  fact,  there  exists  a  positive  c  that  is  independent  of  A  and  such  that 

ll/A-/nllH<Sinf  ||A-/||„  forn>N,  (3.13) 

/€  Wn 


in  which  H„  is  the  linear  span  of  {hifc}J_j. 

Proof:  Retrace  the  reasoning  adduced  in  the  proofs  of  Theorem  1  and  Remark  3®  of  Ref.  6,  mutatis  mutandis, 
replacing  the  “bilinear”  form  B  and  linear  functional  L  of  Ref.  6  with  the  sesquiiinear  form  a  and  conjugate- 
linear  functional  A,  and  accounting  for  the  minor  difference  between  the  form  of  the  variational  problem  given 
in  Eq.  (1)  of  Ref.  6  and  that  of  (3.1).  The  error  inequality  (3.13)  results  from  an  inspection  of  the  proof  of 
Theorem  1  of  Ref.  6.  We  have  written  “bilinear”  here  in  referring  to  the  form  B  of  Ref.  6  because  there  seems 
to  be  some  confusion  in  terminology,  so  one  must  take  some  care  in  working  through  the  developments  of 
Ref.  6.  Evidently  a  complex  Hilbert  space  is  considered  there,  while  the  form  B,  hypothesized  as  “bilinear,” 
seems  to  be  treated  as  sesquiiinear  in  all  of  the  reaisoning  subsequent  to  the  proof  of  Theorem  1  of  Ref.  6.  It 
is  possible  that  the  term  “bilinearity”  in  Ref.  6  is  actually  used  to  indicate  the  property  that  is  more  properly 
designated  “sesquilinearity” ;  if  this  is  so,  then  there  has  been  omitted  a  sign  of  complex  conjugation  in  Eq.  (3) 
of  Ref.  6.  With  the  latter  change,  one  can  read  Ref.  6  with  the  term  “bilinear”  replaced  by  “sesquiiinear" 
throughout.  ■ 

When  M(<^)  is  nontrivial  and  the  solvability  condition  (3.8)  is  fulfilled,  Hildebrandt  and  Wienholtz  [6] 
base  their  method  for  approximating  a  solution  of  (3.1)  upon  application  of  the  Galerkin  method  to  the 
“symmetrized”  problem 


find  f  £  H  satisfying  5-(/',fi)  =  A/i  for  all  h  G  H, 


(3.14) 
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in  which  cr  is  the  (Hermitian-symmetric)  form  corresponding  to  TgT’,  given  by  ff{g,h)  :=  {TgTgg,h)„  = 
^{T^g,  h).  They  show  that  the  Galerkin  method  can  be  applied  to  the  problem  (3.14)  to  produce  approxima¬ 
tions  converging  to  the  unique  solution  of  this  problem  that  also  lies  in  the  subspace  M{T*)^  =  M'{TgTg)^-, 
if  f'^  is  this  solution,  then  /  :=  T*/)(  is  a  solution  of  the  original  problem  (3.1),  unique  to  within  an  element 
oi X{Tg).  Although  this  scheme  may  be  the  best  one  available  in  the  general  case,  it  has  some  unattractive 
features,  t.g.,  one  must  know  how  to  determine  the  action  of  T*  on  the  selected  coordinate  functions.  Of 
course,  a  method  based  upon  computation  of  values  of  the  original  form  a  would  be  much  preferred;  in  at 
least  one  special  case,  viz.,  when  A/i(<r)  =  A/i(<T*),  such  a  method  can  be  employed  (note  that  this  condition 
is  automatically  satisfied  if  a  is  Hermitian  symmetric).  Still,  for  a  complete  resolution  of  the  problem  one 
can  evidently  not  escape  the  necessity  for  determining  Mi(cr),  although  some  information  can  be  salvaged 
without  addressing  the  latter  task;  cf.  the  remarks  following  the  proof  of  the  next  theorem. 

As  a  preliminary,  let  us  recall  some  terminology.  Suppose  that  Ad  is  a  subspace  of  H  that  is  comple¬ 
mentary  to  the  subspace  A/i(<t),  i.e.,  such  that  (3.9)  holds.  Then  the  linear  operator  Pj^  :  H  H  defined 
by 

Pj^h  :=h2,  wherein  h  =  hi  +  with /ii  6  ^((r),  ^2  6  A4,  for  each /i  G  ff 

is  bounded  and  also  a  projection  (i.e.,  satisfies  =  PmI  called  the  operator  of  projection  onto  M  along 
A/'i(«r).  Obviously,  /^(Pj^)  =  A/’i(£r),  while  I  —  Pj,^  is  also  a  bounded  projection  operator,  carrying  H  onto 
(“along  Ad”),  with  A(/  —  P/^)  =  Af.  Of  course,  Pj^  is  self-adjoint,  i.e.,  an  operator  of  orthogonal 
projection,  iff  Ad  =  Ai(<r)-‘-. 

Theorem  3.3.  Let  the  sesquilinear  form  a  on  H  be  the  sum  6  +  k  of  bounded  sesquilinear  forms,  with  6 
definite  and  k  compact,  and  suppose  that  the  equality 

M(<t)  =  M(<T*)  (3.15) 

is  known  to  hold.  Let  M  be  a  subspace  of  H  that  is  complementary  to  Ai(cr),  as  in  (3.9),  and  denote  by  Pj,^ 
the  operator  of  projection  onto  Ad  along  Ai(cr).  Suppose  that  A  is  a  bounded  conjugate-linear  functional 
on  H  for  which  the  solvability  condition  (3.8)  is  valid.  Then  the  unique  G  M  satisfying  (3.10)  can  be 
convergently  approximated  in  the  norm  cf  H  by  use  of  the  Galerkin  method  applied  to  the  form  cr.  More 
precisely,  let  be  a  family  in  H  such  that  is  linearly  independent  and  complete  in  the 

Hilbert  space  (Ad,  (■,  ■)^).  Then  there  exists  N  such  that  for  n  >  N  the  system  (3.12)  possesses  a  unique 
solution  (4"^)lfc=i,  ^nd  the  resultant  sequence  (/„  :=  converges  to  in  the  norm 

of  H .  In  fact,  there  exists  a  positive  c  independent  of  A  and  such  that 

II/aa<  “ /nil// <  c  inf  \\f\M~f\\H  forn>N,  (3.16) 

with  Hn  denoting  the  linear  span  of 

Proof  :  As  noted  in  the  statement  of  the  Theorem,  (Af,  (■,  •);/)  is  a  Hilbert  space,  since  Ad  is  assumed  to  be 
closed  in  H.  With  and  denoting,  respectively,  the  restrictions  of  <t  to  Ad  x  Ad  and  A  to  Ad,  consider 

M(<7jv^)  :=  {  S'  e  Ad  I  <Tji^(g,  h)  =  0  for  all  /i  G  Ad  }. 

Suppose  that  g  G  Ai((T^):  then  we  compute 


(T{g,h)  =  <T{g,  P^/i  -f  (/  -  P^)/i)  =  <T{g,{I  -  PmW  =  0  for  ab  h&H, 

since  (/  -  Pj^)h  G  Ai(<r)  =  Ai((T*),  by  (3.15).  Thus,  g  £  M  O  =  {0},  and  so  Ai((t^)  is  the 

trivial  subspace.  Consequently,  it  is  clear  that  all  of  the  required  hypotheses  are  fulfilled  so  that  we  can 
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apply  Theorem  3.1(tt)  and  Theorem  3  2,  with  H,  it,  and  A  replaced  there  by  (Ad,(-,  )^),  and  A^, 
respectively.  We  conclude  first  that  there  is  a  unique  solution  of  the  problem 

find  f&M  such  that  <»'x(/ih)  =  A^h  forallh^M; 

this  unique  solution  clearly  coincides  with  the  unique  element  €  Af  satisfying  (3.10),  since  the  latter 
holds,  in  particular,  for  every  h  £  M.  Further,  according  to  Theorem  3.2  and  the  given  properties  of  the 
family  there  exists  N  such  that  for  n>  N,  the  system 

n 

P l=l,...,n,  (3.17) 

l:  =  l 

has  a  unique  solution  (^t"^)fc=i,  and  the  resultant  sequence  (/„  :=  52t=i  converges  to  f/^j^ 

in  the  norm  of  Ad,  i.e.,  in  the  norm  of  H.  But,  because  of  (3-15)  and  (3.8),  we  find  that 

=  <{J-  PM)f^k  +  PMhk,{I-  PM)f^l  +  PMf^t))=  <^(*4.  */),  *,  /  =  1,  2,  3,  . .  . , 


and 

f^MPM^i  =  A((7  -  P^)/i,  +  Pj^h,)  =  Ah,,  /  =  1,2, 3, . . . , 
so  the  system  (3.17)  coincides  with  the  system  (3.12).  Finally,  the  error  estimate  (3.16)  follows  from  (3.13).  ■ 

We  conclude  this  section  with  a  remark  that  will  be  important  for  the  application  of  Theorem  3.3 
that  we  have  in  mind,  if  we  should  content  ourselves  in  the  original  fluid-elastic  interaction  problem  with 
the  approximation  of  the  fluid  field  alone.  Maintaining  the  setting  and  notation  of  this  Theorem,  suppose 
that  ii  is  some  normed  linear  space,  L  :  H  —*  H  is  a  bounded  linear  operator,  and  we  are  interested  in 
approximating  in  the  norm  of  If  it  should  be  the  case  that 

A7i(<t)  C  A(L),  (3.18) 

then  obviously 

=  forn>A,  (3.19) 

ifc=i  Jfc=i 

while  the  sequence  (Lfn)^-iv  converges  to  in  the  norm  of  H.  In  our  application,  we  shall  have 

H  =  ,  <T  =  (Tu,  H  =  /f  a(r+),  and  L  given  by  u  i-»  u|r^.  •  n+,  for  u  €  (cf.  Lemma  2.1),  choices  for 

which  (3.15)  and  (3.18)  shall  be  shown  to  hold.  The  point  here  is  that  when  (3.15)  and  (3.18)  are  fulfilled,  we 
can  convergently  approximate  with  very  little  explicit  information  about  the  subspace  A/, (it),  clearly 

a  desirable  feature.  Indeed,  we  need  only  ensure  that  our  choice  of  coordinate  family  {hn}^i  is  such  that 
{ linearly  independent  (i.e.,  has  the  property  that  the  only  linear  combination  of  a  finite  number 
of  elements  of  {/»n}“=i  that  belongs  to  A'i(it)  is  the  trivial  linear  combination,  with  all  coefficients  zero)  and 
complete  in  (Ad,  (•,  O^,);  in  the  later  application,  we  will  even  know  that  the  required  completeness  property 
follows  once  we  have  selected  a  family  {hn}“=i  that  is  complete  in  H  itself  (cf.  Lemma  4.3,  infra).  Then, 
according  to  (3.12)  and  (3.19),  the  computations  necessary  to  construct  approximations  converging  in  H  to 
Lf/^M  1*®  accomplished  without  further  regard  for  the  exact  nature  of  A/,  (it).  In  general,  however,  for  the 
computation  of  approximations  converging  in  H  to  /^x  itself  (which  will  correspond  in  our  application  to 
the  approximation  of  the  elastic-field  portion  of  the  solution  of  the  interaction  problem),  evidently  we  cannot 
avoid  confronting  the  difficult  and  expensive  eigenvalue  problem  that  must  be  solved  to  identify  A/i(cr)  and 
so  determine  how  to  compute  the  action  of  the  projection  operator  Pj,^.  Of  course,  if  we  were  to  know  by 
some  other  reasoning  that  Ai  (<t)  =  {0}  (which  will  correspond  in  the  upcoming  application  to  A/n„,t„  =  {0}), 
then  these  difficulties  disappear. 
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4.  SOLVABILITY  AND  APPROXIMATION-OF-SOLUTION  RESULTS  FOR  THE  WEAK 
PROBLEMS 

In  §2,  we  specified  how  a  solution  of  the  weak  version  of  (P)  is  to  be  obtained  from  a  weak  solution 
of  the  purely  interior  problem  corresponding  to  (Po)  with  data  derived  from  that  of  (P).  Thus,  we  study 
exclusively  the  weak  form  of  the  interior  problem  corresponding  to  (Po)-  Our  first  goal  in  this  section  is 
a  description  of  the  state  of  affairs  concerning  the  solvability  of  the  latter  problem  and  the  applicability  of 
the  Galerkin  method  for  the  generation  of  a  sequence  converging  to  the  unique  solution  when  the  problem 
is  uniquely  solvable.  We  are  restricting  our  considerations  here,  sis  in  §2,  to  the  case  in  which  to  G  L2(ro), 
although  the  statement  and  proof  of  the  following  Theorem  4.1  can  be  modified  in  an  obvious  fashion  to 
permit  the  weaker  hypothesis  to  G  H~3(ro). 

Theorem  4.1.  Recall  the  regularity  conditions  placed  on  flo,  the  conditions  ReK:+  >0  and  ImK+  >  0 
fulfilled  by  k+,  and  the  restrictions  X  >  0,  p  >  0  imposed  upon  the  Lame  parameters. 

(i)  The  dimension  of  the  linear  manifold  finite.  In  particular,  u  closed  in  both  H°  and 

and  we  have  the  direct-sum  decompositions  given  in  (2.26)  and  (2.27).  Thus,  there  is  defined  the 
(bounded)  operator  — <•  of  projection  onto  the  subspace  ^long 

(ii)  There  exists  a  weak  solution  of  the  purely  interior  problem  corresponding  to  (Po(foi  toi  ¥’‘))  >  ^  i^ 

Definition  2.3,  iff  the  associated  conjugate-linear  functional  on  vanishes  on  i.e.,  iff 

/  (t+ -  (t+ •  n+)n+)  •Uo,|r+dAr+ +  /  t_  •  i^lr.  <fAr_  +  /  fo  Ua,dA3  =  0 

Jr+  Jr.  Joo 

for  every  6  An, ,« ■  (^ •  1 ) 

When  this  condition  holds,  there  exists  precisely  one  weak  solution  u  (£  Adn„  u)  problem;  the 

w  -norm  of  this  solution  has  the  bound  given  by 

I|u||h>  <  M||A,.||  <  M  {llfollHo  +  lltollMr,)  +  !r+  -  . 

(4.2) 

for  some  positive  M  that  is  independent  of  the  particular  set  of  data  {fo,  to,  chosen  as  in  Terminol¬ 
ogy  2.1  and  satisfying  (4.1);  here,  IjAioH  is  the  norm  of  A^,  as  an  element  of  the  anti-dual  of  H\  while 
IIAk^II  is  the  norm  of  acting  in  L2(r+). 

(in)  Suppose  that  the  solvability  condition  (4.1)  holds.  Then  the  unique  weak  solution  u  of  the  purely 
interior  problem  corresponding  to  (Po(f<,,to,  ¥>*))  can  be  convergently  approximated  in  the  norm  of 
by  using  the  Galerkin  method.  More  precisely,  suppose  that  is  a  family  in  such  that 

is  linearly  independent  and  complete  in  the  subspace  A4n„  ofH';  here,  completeness 
of  in  will  suffice  to  imply  the  completeness  of  {P;niw„}^_j  in  Adn<,,w-  Then  there  exists 

a  positive  integer  N  such  that  for  n  >  N,  the  system 

n 

^o’u,(wifc,W()^it  =  A„W(,  /=l,...,n,  (4.3) 

possesses  a  unique  solution  (^t"^)lt=i>  resultant  sequence  (u„  ;=  Xrt=i )“=N  con¬ 

verges  to  u  in  the  norm  of  H^.  Moreover,  the  sequence  (Un|r+  •  n+)“_;v  converges  to  ulr+  •  n+  in  the 
norm  of  ff  3(r+),'  since 

n 

«n|r+ •«+ =  11+  forn>N,  (4.4) 

t=i 
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if  follows,  In  view  of  the  form  of  the  systems  in  (4.3),  that  this  sequence  can  be  constructed  even  in 
the  absence  of  explicit  information  concerning  the  subspace  „  and  the  operator  Pm^-  Finally,  the 
following  error  bounds  hold  for  n>  N,  c  and  P  denoting  positive  numbers  that  are  independent  of  the 
particular  set  of  data  {fo ,  to ,  ¥>* }  chosen  as  in  Terminology  2. 1  and  satisfying  (4.1): 

||«-«n||H-<c  ||u-v||h>, 

vePj^iW, 


and 

||(u  -  u„)|r+  •  n+ll^  ll«  -  vIIh-  , 

wherein  Wn  denotes  the  linear  span  of 

Concerning  the  operator  let  us  interject  here  an  explanatory  note:  since  ^  C  ,  it  is  clear  that 

the  operator  on  of  orthogonal  projection  onto  (A/^  carries  the  linear  manifold  of  elements  in  H* 

into  itself;  in  view  of  (2.25)  and  (2.27),  the  restriction  of  this  H°-orthogonal  projection  operator  to  the  linear 
manifold  of  elements  in  when  regarded  as  acting  in  H*,  is  just  the  (nonorthogonal  projection)  operator 
that  we  have  denoted  by  Pm^- 


We  shall  make  some  remarks  on  the  statements  of  Theorem  4.1,  supposing  that  it  has  been  proven. 
Assuming  that  the  solvability  condition  (4.1)  is  fulfilled,  it  is  asserted  that  there  is  a  unique  weak  solution 
u  of  the  purely  interior  problem  for  (Po)<  and  a  method  is  provided  for  constructing  a  sequence  converging 
to  this  solution  in  H*.  To  actually  base  a  computation  on  this  method,  we  must  employ  a  family  {w„}^j 
for  which  the  corresponding  family  of  projections  is  linearly  independent  and  complete  in  the 

subspace  of  Once  we  have  secured  such  a  family  the  numerical  solution  of  the  systems 

in  (4.3)  requires  no  further  knowledge  about  the  subspace  A/n„u,)  evidently  we  cannot  avoid  having  to 
find  the  projections  PM'^n  for  n  =  1, 2, 3, . . . ,  if  we  insist  on  approximating  the  weak  solution  u  itself,  be., 
if  we  wish  to  compute  the  Un  for  n  >  A.  However,  a  further  implication  of  the  Theorem  is  that  we  shall  be 
able  to  approximate  convergently  in  /f  a(r+)  the  normal  component  u|r.^  n^.  of  the  restricted  trace  of  u,  i.e., 
that  we  can  compute  the  u„|r+  -n^,  without  troubling  to  find  those  projections.  In  turn,  we  shall  be  able 
to  approximate  the  function  ^  €  W(n+; «+)  satisfying  (2.13),  the  fluid-field  portion  of  the  weak  solution  of 
(P)  itself,  by  using  Proposition  2.1  (cf.  (2.16)).  Consequently,  we  shall  be  able  to  calculate  approximations 
to  the  fluid-field  portion  of  the  fluid-elastic  interaction  with  little  information  about  the  possible  complex 
amplitudes  of  nonradiating  modes  for  and  w.  It  suffices  to  ensure  that  our  chosen  family  is,  say, 

complete  in  H'  and  such  that  is  linearly  independent,  ».e.,  such  that  any  linear  combination 

of  a  finite  number  of  elements  chosen  from  that  lies  in  A/’n^  is  necessarily  the  trivial  combination 

with  all  coefficients  equal  to  zero  (cf.  Lemma  4.4,  infra).  We  shall  remark  later,  in  Proposition  4.1,  on  the 
selection  of  coordinate  functions  fulfilling  this  latter  requirement. 


The  proof  of  Theorem  4.1  will  be  carried  out  by  checking  that  we  can  apply  to  the  present  setting  the 
abstract  developments  of  the  preceding  §3;  this  will  be  facilitated  by  the  establishment  or  citation  of  certain 
preliminary  results,  to  which  we  proceed.  The  first  statement  is  termed  “the  coercivity  of  strains”  by  Necas 
and  Hlavacek  [8],  but  is  more  commonly  known  as  “Korn’s  second  inequality”  (cf.,  e.g.,  Fichera  [14]). 

Lemma  4.1.  Recall  that  Qo  is  assumed  to  be  a  Lipschitz  domain.  There  exists  a  positive  number  c„, 
depending  upon  only  Qo.  such  that 


L 


Uik[y]ejk[y]  -f-  V  •  v}  dAa  >  Co||v| 


i^> 


for  every  v  G 


(4.5) 


Proof  :  The  statement  is  proven  in  Theorem  3.4  of  Ref  8  for  the  case  in  which  the  underlying  Sobolev  space 
comprises  R^-valued  functions.  But  it  is  clear  that  this  also  suffices  to  cover  the  present  situation,  in  which  the 
elements  of  H*  are  C^-valued,  since  we  have,  for  v  G  H*,  the  equalities  £jt[v]£:jit[v]  =  ejifc[Re  v]£;ifc[Re  v]  -|- 
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Sji[Imv]£ji[Imv]  and  ||v||^n  =  ||Rev||^fc  +||Iinv||„ik  for  A:  =  0, 1.  Alternately,  one  can  appeal  to  the 
reasoning  in  Ref.  14,  carried  out  under  the  assumption  that  the  domain  fio  satisfies  merely  a  “restricted  cone 
hypothesis”  (and  again  for  R^-valued  functions).  ■ 

For  the  sesquilinear  form  defined  on  by  (2.20),  the  associated  “linear  manifolds  of  degeneracy” 
M(<^w)  and  vVi  (K)  in  are  defined  as  in  §3.  Recall  the  definition  (2.24)  of  the  linear  manifold  Aq  of 
complex  amplitudes  of  weak  nonradiating  modes  in 

Lemma  4.2. 

Proof  :  Let  u„  €  ^/'\{cru),  so  that  €  H*  and  eTu{nu,\)  —  0  for  every  v  G  In  particular,  (Ti„(u,,;,  u,^)  = 
0,  and  this  implies,  in  turn,  from  (2.20),  that 

Im  (K+A,c+(ua,|r+  •  n+),  u^(r+  •  n+)i,(r+)  =  0.  (4.6) 

By  (2.1),  we  must  have  Uu,lr+  •  n+  =  0,  so  that  u<„  belongs  to  (cf.  (2.23)).  But  then  (2.22)  is  true, 

as  well  (cf.  the  form  of  cr^  in  (2.20)),  whence  the  inclusion  G  results.  Thus,  Afi{cr^)  C 

The  reversed  inclusion  obviously  holds,  so  A/i(<t„)  =  Similarly,  the  assumption  G  M(<7*),  t.e., 

Uw  €  and  <7,.  (v,  u^,)  =  0  for  every  v  G  H^,  leads  again  to  (4.6),  from  which  the  reasoning  can  be  carried 
through  as  before,  with  trivial  modifications,  to  arrive  at  the  equality  A/’i(cr;)  =  A/^  in  view  of  the  fact 
that  the  first  integral  appearing  in  the  definition  of  generates  an  Hermitian-symmetric  form  on  H'.  ■ 

Lemma  4.3.  If  {w„}“_i  is  a  family  complete  in  HS  then  the  corresponding  collection  is 

complete  in  the  Hilbert  space  (•,  •)h>)>  therein  Pj^^\  :  — ►  H*  denotes  the  operator  of  projection 

ontoMli^  ,^  along 

The  statement  of  Lemma  4.3  anticipates  the  fact  that  has  the  direct-sum  decomposition  given  in 
(2.27).  Therefore,  we  choose  to  defer  the  proof  of  this  Lemma  until  after  we  have  established  that  ^  is 
finite-dimensional  in  the  course  of  the  proof  of  Theorem  4.1  (without  using  Lemma  4.3!). 

Proof  of  Theorem  f.l:  At  the  outset,  we  observe  that  the  sesquilinear  form  cr^  and  the  conjugate-linear 
functional  A^,  are  obviously  bounded  on  H'.  We  wish  to  show  that  (T„  can  be  written  as  the  sum  of  a  definite 
form  and  a  compact  form,  i.c.,  that  it  satisfies  the  fundamental  hypothesis  required  in  the  developments  of 
§3.  To  this  end,  let  us  define  the  sesquilinear  forms  6{-,  •)  and  Ku{-,  ■)  on  by  setting 


dXs 


6{\i,v):=  f  {£:;^_^(u,v)-(-2pu  v} 

«a/(u,  v)  :=  -(2p  +  u  vdXs  4-  q+c\k\  J  (A,^(ujr+  •n+))v|r+  n+dXp^  J 


for  u,  V  G 


Clearly,  =  6  +  k^,-  The  definiteness  of  (the  Hermitian-symmetric  form)  6  follows  from  Lemma  4.1,  upon 
recalling  that  A  >  0  and  p  >  0,  for,  we  can  write 


6(v,v) 


+  2pej  j:[v]ejjt[v]  4-  2/iv  •  v)  dXs  >  2p 


4-  V  •  v}  dXs 


>  2/JCo||v||^,  for  V  G  (4.7) 

The  compactness  of  is  a  simple  consequence  of  the  compactness  of  both  the  natural  injection  map  carrying 
into  H°  and  the  operator  v  i— ►  v|p^  ■  n.(.  taking  into  L2(r+)  (for  the  latter,  cf.  Lemma  2.1(i)),  coupled 
with  the  boundedness  of  on  L2(F+).  We  remark  that  the  compactness  of  is  not  needed  here,  but 
could  have  been  used  to  produce  the  same  conclusion  if  we  knew  only  that  v  t— >  v|p^  •  n+  is  bounded  from 
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into  L2{r+).  Thus,  we  can  apply  Theorem  3.1.  By  doing  so,  we  find  first,  from  statement  (i)  of  that 
Theorem,  that 

dimWo^ =  d\mXi{au)  =  dimM(<r*)  <  oo,  (4.8) 

having  also  used  Lemma  4.2  to  get  the  first  equality  here.  Moreover,  by  Theorem  3.1(i*),  there  exists  some 
u  G  satisfying  (2.28)  iff  jVi(o’w)  C  A'’(A„),  or,  again  by  Lemma 4.2,  iff  ^  C  X{S.u),  an  inclusion  that 
is  equivalent  to  the  explicit  condition  given  in  (4.1),  in  view  of  the  definitions  in  (2.21),  (2.23),  and  (2.24). 
Now  we  know,  in  particular,  that  closed  in  (and  in  H*),  so  (2.26)  holds,  whence  (2.27)  follows; 

■^ho.u  subspace  of  that  is  complementary  to  A/’i(a’^)  (=A/n^  y,).  By  the  direct-sum  decomposition  of 
(2.27),  there  is  induced  the  (bounded)  projection  operator  Pm-  of  onto  along  (= 

Taking  M  =  -^no.u/  ***  Theorem  3.1(ui),  we  can  assert,  when  (4.1)  obtains,  that  there  exists  precisely  one 
u  6  that  fulfills  the  requirement  (2.28),  and  so  provides  the  unique  weak  solution  of  the  purely 

interior  problem  corresponding  to  (Pq);  (3.11)  implies  that  the  H^-norm  of  this  element  has  the  bound  given 
in  (4.2)  (with  M  independent  of  the  data  generating  A^,),  since  it  is  easy  to  see  from  (2.21)  that  !|Au,||,  the 
norm  of  as  an  element  of  the  anti-dual  of  H^,  is  less  than  the  number  in  brackets  on  the  right  in  (4.2). 
This  completes  the  proof  of  (i)  and  (it). 

To  establish  (*ii)i  in  which  we  are  given  that  (4.1)  holds,  we  denote  by  u  G  u,  the  unique  weak 
solution  of  the  purely  interior  problem  for  (Po).  We  wish  to  invoke  Theorem  3.3,  making  the  obvious 
identifications  H  =  a  =  AA  ■=  eic;  this  is  permissible  by  the  facts  discovered  in  the  proof  of 

(t)  and  (ii)  and  by  Lemma  4.2  (giving  (3.15)).  With  the  additional  use  of  Lemma  4.3,  all  of  the  conclusions 
of  (iii)  now  follow  immediately.  For  example,  the  convergence  of  the  sequence  (u„lr^  •  n+)“_jv  to  u|r+  •  n+ 
in  the  norm  of  //^(r+)  results  simply  from  the  convergence  of  (u„)^jy  to  u  in  the  norm  of  and  the 
boundedness  of  the  operation  v  i->  •  n+  on  into  i/3(r+)  (Lemma  2.1(ti)).  The  equalities  of  (4.4) 

are  true  because  (I  —  Pjv^i)wi  lies  in  for  each  k,  so  that  wt|r+  •  n+  =  (PM‘Wifc)|r+  •  n+.  Concerning 

the  error  estimates  displayed,  the  first  is  a  direct  result  of  (3.16),  and  the  second  follows  from  the  first  (by 
Lemma  2.1(<>)).  ■ 

Finally,  we  give  the 

Proof  of  Lemma  4-3:  Recall  the  Hermitian-symmetric  sesquilinear  form  8  defined  on  H'  in  the  proof  of 
Theorem  4.1.  Inequality  (4.7)  shows  first  that  ^(•,  •)  is  an  inner  product  on  HS  and  then,  when  combined 
with  the  boundedness  of  5,  that  the  norms  induced  by  ^(•,  •)  and  (■,  )h>  on  the  complex  linear  space  of 
elements  in  are  in  fact  equivalent: 

2pCo|lv||^.  <  6(v,v)  <  c(,||v|l|,,  for  v  G 

for  some  positive  c(,.  By  Hj,  we  denote  the  collection  of  elements  in  equipped  with  the  inner  product 
8(-,  •).  Then  Hj  is  complete  along  with  H^,  and  the  two  Hilbert  spaces  are  isomorphic  qua  Banach  spaces; 
a  family  is  complete  in  iff  it  is  complete  in  Hj.  The  linear  manifold  Ad^^  ,^  is  closed  in  Hj,  since  we 
already  know  that  it  is  closed  in  ;  a  family  is  complete  in  the  Hilbert  space  (A4^^  (•,  •)h>)  iff  that  family 

is  complete  in  the  Hilbert  space  i(-,  ■)). 

Let  us  suppose,  for  the  moment,  that  the  inclusion 

C  (4.9) 

has  been  shown  to  hold,  wherein  (A/n„,^)f-tb  indicates  the  orthogonal  complement  of  A/'j^^  ^  in  Hj,  and 
check  that  the  claim  of  Lemma  4.3  follows  therefrom.  Thus,  suppose  that  {Wn}^i  is  a  family  complete 
in  H*;  then  it  is  also  complete  in  Hj.  Let  u  be  an  element  of  ^  that  is  6-orthogonal  to  the  family 
of  projections: 

(5(u,  P^iWn)  =  0  for  n  =  1,2,3,... 
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Then,  since  (I  —  PM')'^n  €  >  1.  with  (4.9)  we  get 

5(u,w„)  =  0  for  n  =  1,2, 3, . . .  , 


and  so  conclude  that  u  =  0.  This  shows  that  the  collection  is  complete  in  ^(■.  ))) 

whence  it  is  also  complete  in  (•,  •}h»)- 

Therefore,  for  the  completion  of  the  proof,  it  is  sufficient  to  verify  (4.9).  Then,  let  u  6  i.e., 

by  (2.25),  let  u  G  H'  and  dXs  =  0  whenever  u^,  €  As  a  consequence,  according  to  the 

definition  in  (2.24),  we  must  have 

/  €x,^iu^,n)dX3  =  {fA,^(Uu,,u)  -  u)  dAa  =  0  for  every  Uo,  e  •  (4-10) 

Vflo  J{1„ 

Clearly,  then,  (5(u,Uu)  =  0  whenever  is  in  An„,a-,  implying  that  (4.9)  is  correct. 

Before  leaving  the  proof,  we  note  parenthetically  that  equality  actually  obtains  in  (4.9); 


A4n..u,  = 

Indeed,  if  we  suppose  that  u  G  (A/”^^  **  ^ 

/  {fA,^(ua,,u)  +  2pu„  ujdAs  =  0  for  every  G  , 

Jflo 


then  the  latter  equality  together  with  the  second  equality  in  (4.10)  show  that  •  udAs  =  0  for  every 

“u,  €A/'j^„„,  so  uG(A/'j^,,„)(-^)‘>nHi  =Ad^^  „.  ■ 


Let  us  return  to  the  linear-independence  question:  we  wish  to  know  useful  sufficient  conditions  on  a 
family  {wn}^=i  C  under  which  we  shall  be  certain  that  is  linearly  independent.  It  is  quite 

simple  to  prove  that  this  linear  independence  obtains  for  any  such  family  that  is  itself  linearly  independent 
and  comprises  piecewise-polynomials,  eis  in  the  usual  applications  of  the  finite-element  methods.  We  begin 
by  recording  an  obvious  necessary  and  sufficient  criterion  for  the  linear  independence  of  a  collection  of 
projections  {Pai‘W„}^_i. 

Lemma  4.4.  Let  {w„}^i  be  a  family  in  and  Pm‘  the  operator  in  of  projection  onto  AdJj  along 
Then  is  linearly  independent  iff  for  each  positive  integer  N,  the  inclusion  Wn  € 

A^n„  w  some  {a„}^_i  C  C  implies  that  ai  =  ■  ■  •  =  ajv  =  0.  In  particular,  the  linear  independence  of 
{w„}^j  is  necessary  for  the  linear  independence  of 

Proof:  For  the  necessity,  suppose  that  is  linearly  independent,  let  N  he  a  positive  integer, 

and  let  be  a  complex  set  with  23^=1  G  Then  a„Pji^iw„  =  0,  whence  aj  — 

•  ■  •  =  =  0,  by  the  linear  independence  of  the  family  of  projections.  In  particular,  the  vanishing  of 

a„w„  gives  oi  =  •  •  •  =  aw  =  0,  so  is  linearly  independent.  For  the  sufficiency,  let  the 

condition  hold,  and  suppose  that  anPM''^n  =  0,  for  some  positive  integer  N  and  complex  set  {a„}^_2. 

Then  ^  so  we  must  have  aj  =  •  •  •  =  aw  =  0,  implying  that  is  linearly 

independent.  ■ 


As  usual,  by  a  (C^-valued)  polynomial  in  we  mean  a  function  on  of  the  form  x  i— ►  53o<|cr|<K  Oox“. 
for  some  nonnegative  integer  K  and  collection  {Ca}o<)a)<A-  lying  in  C^;  the  standard  notation  x“  := 
Xi'z2^x^^  is  used  here,  for  a  3-index  a  =  (ai,Qf2,ar3)  and  an  x  G  R^-  A  C^-valued  function  w  defined 
As-a.e.  in  fio 's  said  to  be  a  piecewise-polynomial  in  iff  there  exists  a  family  of  pairwise-disjoint 

open  subsets  of  n„  such  that  the  complement  fi<,\Uf=i  Aa-measure  zero  and  each  of  the  restrictions 

vf\Oj,  j  =  1, . . . ,  J,  coincides  in  Oj  with  some  polynomial  (perhaps  depending  upon  j). 
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Proposition  4.1.  Let  the  fctmily  iie  in  H*,  and  suppose  that  each  w„  is  a  piecewise-polynomial 

in  fio-  Then  {Pa4iw„}^j  is  linearly  independent  iff  {w„}^j  is  linearly  independent. 

Proof  :  By  Lemma  4.4,  we  need  only  demonstrate  the  sufficiency.  Then,  let  {w„}^j  be  linearly  independent. 
Let  iV  be  a  positive  integer,  C  C,  and  suppose  that  Yin=i  6  according  to  Lemma  4.4, 

we  must  show  that  ai  =  •  •  ■=:  ajv  =  0.  Now,  as  already  pointed  out  in  §2,  each  element  of  ^  has  real- 
analytic  real  and  imaginary  parts,  and  is  a  classical  solution  of  (127),  so  a„w„  has  these  properties. 

Since  each  w„  is  a  piecewise-polynomial  in  Qo,  it  is  easy  to  see  that  we  can  find  an  open  subset  of  flo 
in  which  coincides  with  a  polynomial.  Therefore,  by  the  uniqueness  theorem  for  real-analytic 

functions  (cf.,  e.g.,  Ref.  16),  a„w„  must  coincide  with  that  polynomial  throughout  Do.  Consequently, 

the  sum  must  vanish  in  Qo,  since  (1.27)  has  no  nontrivial  polynomial  solutions  (because  >  0).  Now  the 
equality  anW„  =  0,  coupled  with  the  linear  independence  of  implies  the  desired  conclusion 

ai  =  ■  •  •  =  UN  =  0.  I 

Lemma  4.3  and  Proposition  4.1  imply  that  the  computation  of  approximations  to  the  fluid-field  portion 
of  the  unique  weak  solution  of  the  fluid-elastic  interaction  problem  (P)  (corresponding  to  data  fulfilling  the 
requisite  solvability  condition)  can  be  carried  out  in  complete  ignorance  of  the  dimension  of  A^j\^  ^  if  we  select 
the  family  of  coordinate  functions  to  comprise  piecewise-polynomials  in  Qo,  taking  care  to  ensure 

only  that  the  particular  such  family  chosen  is  complete  in  and  linearly  independent.  But  the  denseness 
of  in  H'  and  standard  results  from  the  theoretical  foundations  of  finite-element  methods  concerning 
the  interpolation  of  elements  in  by  appropriate  piecewise-polynomials  in  H’  (cf.,  e.g.,  Theorem  8.2.2  of 
Ref.  17)  show  how  to  construct  coordinate  families  of  piecewise-polynomials  in  Qo  that  are  also  complete  in 
H'.  Of  course,  none  of  this  helps  in  overcoming  the  difficulties  involved  in  computing  approximations  to  the 
elastic  part  of  the  weak  solution  of  the  problem  when  either  it  is  known  that  ^  is  nontrivial  or  there  is 
no  a  priori  assurance  that  A^,  ^  =  {O}. 

5.  ON  THE  DERIVATION  OF  THE  INTERIOR  AND  EXTERIOR  LIMITING-AMPLITUDE 
CONDITIONS 

This  section  is  devoted  to  the  promised  motivation  for  accepting  (1.34)  and  (1.36)  as  limiting-amplitude 
conditions,  i.e.,  under  sufficiently  stringent  assumptions  concerning  the  asymptotic  approach  to  a  time- 
harmonic  state,  as  conditions  satisfied  by  the  assumed-to-exist  complex  limiting  amplitudes  u  and  ip  of  the 
(assumed-to-exist)  solution  (U,<I>)  of  the  fundamental  initial-vdue  problem  governing  the  time-dependent 
fluid-elaistic  interaction,  as  described  in  §1.  Actually,  the  interior  and  exterior  eispects  can  be  treated  essen¬ 
tially  separately;  specifically,  we  shall  provide  the  desired  assertions  in  the  form  of  two  theorems,  the  first 
concerning  solutions  of  the  damped-wave  equation  (1.10)  (and  providing  a  basis  for  taking  the  Sommerfeld 
condition  (1.34)  as  the  appropriate  limiting-amplitude  condition  for  any  boundary-value  problem  for  (1.23), 
with  K+  specified),  the  second  involving  solutions  of  the  Navier  equations  (1.11). 

Consistently,  we  shall  denote  by  d+(x)  the  distance  between  a  point  x  6  R^  and  F^.: 

tf+(x)  ;=  dist(x,  r+)  :=  min  |y  -  x|. 

y€r+ 

For  brevity,  in  this  section  we  shall  write  c  and  7  in  place  of  c+  and  7+,  respectively. 
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Theorem  5.1.  Let  f2+  C  be  as  specified  in  §J  and  §2.  Thus,  $7+  is  connected  and  the  complement  of 
the  closure  of  a  bounded  and  connected  open  set  of  type  3?^;  as  before,  set  r+  ;=  dQ+ .  Let  c  >  0  and  y  €  C, 
with  either  7  =  0  or  Re  7  >  0.  Let  $  €  C^(Q+  x  R)  satisfy 


y  1 

-A$  +  ^^.4  +  ^***>44  =  0,  in  fi+  X  R, 


(5.1) 


and  vanish  in  x  (-00,0); 


${-,  t)  =  0  in  for  t  <  0. 


Suppose  that,  for  some  real  positive  u,  can  be  written  in  the  form 


(5.2) 


<^(-,<)  =  4>’'(-,<)  +  ^(Oe  in  Q+  for  each  t  E  R,  (5.3) 

wherein  95  6  C^(fl4.)  and  E  C^(n+  x  R),  with  the  latter  fulfilling  the  following  conditions  of  transience 
and  boundedness: 

litn  4>’^(y,s)  =  0  whenever  y  E  fly  (5-4)i 

J  — ♦OO 

and,  for  certain  positive  numbers  So  and  M, 


lim  'I>[4(y,  s)  =  0  and  lim  grad  $’'(y,  s)  =  0  if  y  E  0+  and  0  <  <i+(y)  <  So, 

5— *00  3— *00 


(5.4)2 


and 


l^’^(y.«)l  <  A/,  |4>[4(y,s)|  <  M,  and  |grad  4>’’(y,  s)|  <  M 

for  y  Efl+  with  0  <  dy(y)  <  So  and  s  >  0.  (5.5) 

Then  <p  has  real  and  imaginary  parts  that  are  real-analytic  in  0+  and  satisfies  both  (1.23)  and  (1.34),  wherein 
the  complex  k+  =  Ky(ui;c,y)  is  computed  as  in  (1-19)  and  (1.20). 


We  remark  that  (5.3)  obviously  implies  that  the  hypothesized  must  have  more  regularity  than  is  indicated 
by  its  inclusion  in  C^(n+  x  R).  Theorem  5.1  is  easy  to  prove  because  the  hypotheses  have  been  adjusted  to 
make  it  so;  our  only  interest  here  lies  in  showing  that  (1.34)  must  hold  when  a  solution  of  (1.10)  evolves  in  a 
“reasonable”  manner  from  a  quiescent  state  to  approach  a  time-harmonic  form  with  the  complex  amplitude 
tp.  A  far  more  difficult  question  concerns  the  formulation  of  conditions  on  the  shape  of  and  the  forcing 
terms  in  the  problem  under  which  this  approach  to  a  time-harmonic  state  is  assured.  Similar  remarks  are 
pertinent  to  the  next  theorem. 

It  is  very  convenient  to  couch  the  staten.ent  and  proof  of  Theorem  5.2  in  terms  of  Hilbert-space-valued 
functions  on  R;  for  the  basic  definitions  and  results  concerning  such  maps,  one  can  consult,  e.p.,  Refs.  18 
and  19.  With  a  separable  Hilbert  space  H ,  we  associate  the  collection  C(R;  H)  of  all  functions  on  R  into 
H  that  are  strongly  continuous,  i.e.,  that  are  continuous  when  H  is  equipped  with  the  topology  induced 
by  its  inner-product  structure.  For  a  positive  integer  m,  by  C’"(R;  H)  we  denote  the  family  of  all  elements 
U  E  C(R;  H)  possessing  m  strong  derivatives  =  U', . . . ,  that  are  also  in  C(R;  H),  “strong”  here 
signifying  that  the  derivatives  are  defined  with  respect  to  the  inner-product  structure  of  H .  For  integration 
of  such  functions  over  subsets  of  R,  we  use  the  Bochner  theory  (cf.  Refs.  18  and  9),  although  the  strong 
Riemann  integral  (cf.  Ref.  19)  will  suffice,  in  view  of  the  continuity  hypotheses  that  we  pose  in  the  upcoming 
assertions  (in  the  interest  of  simplicity).  If  we  were  to  weaken  these  hypotheses  to  allow  for  the  existence  of 
derivatives  in  a  generalized  sense,  then  use  of  the  Bochner  theory  would  be  essential. 
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Theorem  5.2.  With  setting  and  notation  as  set  forth  in  §I  and  §2,  let  the  functions  Fj,  ;  R  — ►  H°, 
To  ;  R  — ►  L2(ro),  and  ^9+  :  R  -+  /li(r+)  vanish  in  (— oo,0).  Let  Fo  belong  to  C(R;H°)  and  To  be 
contained  in  C(R;L2(ro)).  Suppose  further  that,  for  someu  >  0, 


Fo(t)  =  F;(t)  +  'I 

>  for  each  t  G  R, 

T„(t)  =  T:(0  +  toe-‘“'  J 


(5.6) 


wherein  fo  €  H°  and  to  G  L2(ro),  while 


^lim  ||F:(<)||ho  =  ^Hm  ||T;(<)  -  T^t)  ■  n+ii+||Lo(r^)  =  ||T:(0||Lo(r_)  =  0.  (5.7) 

Let  U  :  R  — ►  vanish  in  (— oo,  0),  belong  to  C^(R;  H°)  when  regarded  as  taking  values  in  H°,  and  satisfy 


eoU"(t)  -  A5,^U(0  =  Fo(<)  for  each  t  G  R, 


along  with 


T“+[U(<)]  =  T+(<)  +  ^^+(<)n+ 


1 


for  each  t  G  R. 


T"-[U(<)]  =  T_(t) 

Suppose  that  U  has  the  form 

U(<)  =  U^(<)  +  ue-"‘,  t  G  R, 

for  some  u  €  and  U’’  G  C^(R;  H°),  the  latter  satisfying  the  conditions  of  transience 


lim  ||U’’^''^(<)1|h'>  =  0  for  j  =  0, 1,  and  2. 


»oo 


Then,  for  every  G  u> 


and 


/  to  Uwlr.  +  /  fo  Uo,dA3  =  0, 

Jr,  Jn, 

[Jr  '^ois)-^\r,dXr,  +  Fl{s) dX,'^  e'-‘ ds. 


(5.8) 

(5.9) 

(5.10) 

(5.11) 

(5.12) 

(5.13) 


The  proof  of  Theorem  5.1  can  be  carried  out  in  a  straightforward  manner  with  the  use  of  an  extension 
of  the  classical  Kirchhoff  representation  for  certain  solutions  of  the  usual  wave  equation  ((5.1)  with  7  =  0; 
cf.,  e.g.,  Ref.  20)  to  cover  the  more  general  case  of  the  damped-wave  equation  (5.1).  This  extension  can  be 
established  with  the  aid  of  a  causal  fundamental  solution  for  the  hyperbolic  operator  in  question  (cf.,  e.g., 
Ref.  21).  Alternately,  one  can  proceed  as  in  the  proof  of  the  representation  result  that  we  provide  in  the  form 
of  Lemma  A.l  of  the  Appendix;  it  is  convenient  to  introduce  here  some  additional  notation  in  preparation 
for  the  use  of  that  Lemma.  When  x  G  R^,  let  the  corresponding  distance- function  Tx  on  R^  be  defined  by 

'•x(y)  :=  |y  -  x|  for  y  G  R^. 

Suppose  that  ’J  is  a  complex  function  defined  on  a  cylinder  F  x  R,  for  some  F  C  R^-  With  'I',  c  >  0,  and 
any  chosen  x  G  R^  and  <  G  R,  we  associate  a  function  defined  on  F,  denoted  by  [^](x,()  :  F  — ►  C  and  termed 
the  (c-)reiardation  of  '9  with  respect  to  (x,<),  according  to  the  rule 

•=  fy-  <  -  ^  ^  ^  ■ 
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Recall  our  convention  of  §1:  for  an  appropriate  complex  function  'I'  defined  in  Q+  x  R,  by  we  mean  the 
function  on  r+  x  R  obtained  by  taking  the  spatial  normal  derivative  at  fixed  values  of  the  fourth  argument. 
That  is,  St,a+(-,i)  :=  ^(•)0in+  on  r+  for  <  6  R. 

Proof  of  Theorem  5.1 :  We  shall  first  assume  that  R^\  is  of  type  Si* ,  4>  P  x  R),  6  C^(n+  x  R), 

tfi  G  C^(Q+),  (5.4)i  holds  for  every  y  G  n+,  and  (5.4)2  and  (5.5)  are  valid  for  y  G  r+  instead  of  for  the 
points  y  specified  there.  Notice  that  then  (5.3),  as  well  as  the  corresponding  equalities  between  the  extended 
first  partial  derivatives,  must  hold  in  all  of  for  each  <  G  R-  Under  these  conditions,  we  shall  show  that 
ip  must  be  given  in  0+  by 

Then,  returning  to  the  original  hypotheses,  and  setting,  for  £  >  0, 

(fi+),  :=  {  y  G  I  d+(y)  >  £  }, 

we  shall  have  shown  that  (5.14)  holds  with  Q+  replaced  by  (fi+)t  (and,  of  course,  r+  replaced  by  5(n+)f),  for 
any  sufficiently  small  positive  £.  From  this,  clearly  it  shall  follow  that  ip  has  real-analytic  real  and  imaginary 
parts  in  all  of  that  (1.23)  holds,  and  that  therefore  (1.34)  is  true,  the  latter  obtaining  because  Im/c+  >  0 
(as  pointed  out  in  §1).  Thus,  it  suffices  to  prove  (5.14)  under  the  modified  hypotheses  listed;  in  particular, 
it  is  to  be  understood  for  the  remainder  of  the  proof  that  (5.4)  and  (5  5)  are  taken  to  hold  for  each  y  6  r+. 


With  this  agreement  in  force,  we  have  the  representation  for  $  available  by  appropriately  specializing 
in  the  statement  of  Lemma  A.l.  Specifically,  taking  70  =  0  and  71  =  7  there,  we  get 


=  (^) 


e-7(<-T)/2 


|?7  ((<  -  ^.n+(-.  r) 

+  ^(r^),n+4  (^{t  -  rf  -  ^r^^  $(• 


dAr+ 


for  X  G  n+  and  <  G  R.  (5.15) 


Here,  with  denoting  the  modified  Bessel  function  of  the  first  kind  and  order  ly,  the  function  :  R\{0}  — ►  C 
is  defined  to  vanish  for  negative  values  of  its  argument  and  to  be  given  by 

g.^(s)  :=  for  s  >  0;  (5.16) 

S  3 

then  3-,l(0,oo)  G  C°°([0,oo)).  In  particular,  30  =  0,  so  that  the  terms  in  (5.15)  that  contain  g^  vanish  for 
the  ordinary  wave  equation,  and  this  gives  the  classical  Kirchhoff  representation. 

Now,  let  X  be  a  point  chosen  in  n+;  since  (p(x)  =  limt_oo  e'"‘$(x,t).  we  examine  the  product  e*"‘^'(x,() 
for  t  G  R.  Employing  (5.15)  and  observing  that,  by  (5  3)  and  the  additional  regularity  assumptions  noted, 
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[^](x.t)(y)  =  [‘^’■](x,t)(y)  +  ¥’(y)e  [^,4](,,,)(y)  =  [^r4](x,o(y)  -  ^^(y)«  and 

(^,  n+](,_,)(y)  =  [^rnj(,  ,)(y)  +  V?.  n+(y)c"‘"*c‘“’'«(y)/'  hold  for  y  e  r+,  we  find 


-yr,l2c 


rx 

-K^) 


f--yr. 


'x,n+[^r4](x,«) 


for  t  6  R,  (5.17) 


having  set 


a  =  a(u>,y) 


(5.18) 


In  the  second  integral  over  r+  on  the  right  in  (5.17),  a  simplification  has  resulted  from  the  combination 

'e-7r,/2c\  ,  /^-orr./c\  /p-orr./i:' 


on  Fj 


We  write  (5.17)  in  the  form 

e'"‘  1 

e-‘«>(x,t)  =  I,(x,t)  -  ^2:2(x,t), 

Ii(x,t)  denoting  the  first  integral  over  r+  appearing  on  the  right,  Z2(x, <)  denoting  the  second.  We  intend 
to  show  that 

limli(x,t)  =  0  (5.19) 

(-♦OO 

and  /  •  s  V 

^Ijm  l2(x,t)  =  dAr„  (5.20) 

whence  (5.14)  shall  follow,  and  the  proof  will  be  effectively  complete,  as  outlined.  The  equalities  (5.19)  and 
(5.20)  are  established  by  considering  in  turn  the  various  terms  appearing  in  Ji(x,  t)  and  Z2(x.  i),  ide.ntifying 
for  each  integrand  its  pointwise  limit  as  <  — »  oo  and  an  estimate  valid  for  all  sufficiently  large  t,  in  such  a 
way  that  the  Dominated-Convergence  Theorem  of  Lebesgue  will  be  applicable  to  give  the  limiting  value  of 
that  term.  Introducing  the  number 

d'''(x)  ;=  max  |y  —  x|, 
yer+ 

henceforth,  unless  otherwise  specified,  we  suppose  (at  least)  that  t  >  d‘*'(x)/c,  so  that  t  —  rx(y)/e  is  positive 
for  all  y  6  r+. 
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With  this  plan  in  mind,  we  examine  first  the  integral 

r  /e-7r./2c\ 

2u(x,0  := 

Since  t)(y)  =  ^rn+(y,t  —  »'x(y)/c)  for  every  y  G  r+,  it  is  clear  from  (the  modified  form  of)  (5.4)2 

that  the  integrand  here  converges  pointwise  on  r+  to  zero  as  t  -+  oo.  Moreover,  it  is  just  as  clear  that  the 
integrand  is  majorized  in  modulus  by 

for  all  t  >  </+(x)/c,  the  inequality  here  following  from  (the  modified  form  of)  (5.5).  Consequently,  the 
Dominated-Convergence  Theorem  allows  the  assertion  that  limt^oo  Iii(x,  t)  =  0.  The  integrals  correspond¬ 
ing  to  the  second  and  third  terms  in  the  integrand  of  Ii(x,t)  can  be  shown  also  to  have  limit  zero,  in  a 
similar  manner,  by  using  the  other  (modified)  hypotheses  in  (5.4)  and  (5.5).  Now  (5.19)  has  been  verified  if 
7  =  0,  so  we  assume  in  the  further  reasoning  concerning  Ii(x,t)  that  Re^  >  0. 

Next,  let 


Ii4(x,<)  ;=^  ^  ,.)drdAr+, 


which  can  be  wiitten 


wherein  ho  denotes  the  Heaviside  function: 


L  /  \ _ f  0  if  5  <  0, 

'■»<** -ti  iis>o: 


Once  again,  with  (5.4)2  it  is  evident  that  the  integrand  in  (5.21)  converges  to  zero  as  t  — ►  oo,  pointwise  on 
r^.  X  (0,oo).  To  derive  an  acceptable  uniform  estimate  for  the  modulus  of  that  integrand,  we  need  some 
information  about  the  behavior  of  h  when  its  argument  lies  in  the  right  half-plane  and  has  large  modulus.  A 
standard  asymptotic  result  for  Ii,{z)  when  the  principal  argument  argz  is  restricted  by  |  argz|  <  (ir/2)  —  S, 
with  6  fixed  in  (0, 7r/2)  (cf.,  e.g.,  Eq.  (5.11.10)  of  Ref.  22),  coupled  with  the  fact  that  now  Rej  >  0,  shows 
that  there  exist  positive  numbers  M^iy)  and  s^(y)  such  that 


3R€7/2  r  M  t-vl  1 

|7,,(tsV2)|  < 1+ — for  s>  5^7).  (5.22) 

(xl7l)3s«  L  S3  J 

Thus,  for  y  G  r^.,  r  >  \/(d+(x}/c)^  +  Si(y)^,  and  any  t,  the  modulus  of  the  integrand  in  (5.21)  at  (y,7;<) 
is  bounded  by 


„-TRe7/2. 


p(T’-(r.(y)/c)=’)iRe7/2 
(7r|7|)3(r2  -  (rx(y)/c)2): 


■l<J'rn+(y,<  -  r)ho{t  -  r)|  < 


M 


(x|7|)i(r2-(d+(x)/c)2)r 


Meanwhile,  in  the  remaining  (bounded)  portion  r+  x  (O,  y(d+(x)/c)2  -f-  81(7)2]  of  r+  x  (0,oo),  and  for, 
say,  <  >  0,  it  is  easy  to  see  that  the  integrand  in  (5.21)  is  uniformly  bounded  in  modulus,  in  view  of 
the  boundedness  of  gy  on  bounded  subsets  of  (0,oo)  and  of  on  r+  x  (0,oo).  Consequently,  there  is  a 
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nonnegative  element  of  Li  (r+  x  (0,  oo))  dominating  the  modulus  of  the  integrand  in  (5.21)  for  all  sufficiently 
large  t,  whence  limt_oo  2’i4(x,  t)  =  0. 


Denoting  by  the  final  term  of  Zi(x,t),  an  integration  by  parts  produces 

((«  -  rf  -  ♦'(  ,  r)dr  JAr, 

+  «W . r)+ (I  +  »'(  ..)) 

^’’x)  (5.23) 

with  57(0+)  :=  lim,_o+  gy(s)  =  7/4.  Obviously,  the  integral  corresponding  to  the  first  term  in  the  integrand 
on  the  right  in  (5.23)  can  be  treated  by  essentially  the  same  reasoning  applied  to  Iii(x,  •)  and  so  vanishes 
in  the  limit  as  f  — ♦  00.  For  the  second  term  in  the  integrand,  with  the  aid  of  (5.22)  we  derive  the  estimates 


l(»•x).n+(y)^-j — ^^(y,0)ffy  (^t^-  ^r^(y)"j  I 

g-tRey/2  gRey(t^-(r^(y)/ef)i/2  f 

<  2d+(x) - ^ - 3-  { 

‘  (jr|7|)3(f2_(r,(y)/c)2)^ 


1  + 


M.(7) 


-  (r*(y)/c)2) 


TH«^^(y,o)| 


4d+(x) 


(x|7|)3<(<2  -  (d+(x)/c)2) 


rnapc|$^(z,0)| 


21?  «€r+ 


for  all  y  €  r+  and  all  sufficiently  large  t  >  d+(x)/c,  (5.24) 

showing  at  once  that  the  corresponding  integral  of  this  term  over  r+  has  limit  zero  as  t  00.  The  same 
assertion  can  be  proven  for  the  final  (iterated-integral)  term  on  the  right  in  (5.23)  by  proceeding  along  the 
lines  of  the  analysis  carried  out  for  Ii4(x,  ■);  only  straightforward  modifications  are  required,  so  we  omit  the 
details.  The  limit  equality  in  (5.19)  can  now  be  regarded  as  established. 


Turning  to  T2(x,  •),  it  is  first  of  all  obvious  that  (5.20)  is  true  if  7  =  0,  since  ^0  =  0  and  — a(w,0)/c  = 
iui/c  =  c,  0).  Thus,  we  suppose  in  the  remainder  of  the  proof  that  Re7  >  0.  At  the  outset,  we  perform 

integrations  by  parts  in  the  inner  integral  appearing  in  l2(x,  <),  once  for  the  term  involving  and  twice  for 
the  term  containing  g'^.  In  the  former,  we  base  the  computation  on  the  observation 

^ffy(0-rf-rx(y)/c^)  =  for  0  <  r  <  t  -  rx(y)/c,  with  y  6  r+  and  <  >  d+(x)/c,  (5.25) 

wherein,  for  such  r,  y,  and  t, 

Gy(r;  y,  <)  ;=  lo  Q  ((<  -  r)^  -  ^r^(y))^^  I 


in  the  term  involving  g'y,  a  first  integration  by  parts  results  in  the  appearance  of  g^,  whereupon  (5.25)  is 
used  once  again,  to  effect  a  second  integration  by  parts.  With  A,,  denoting  Macdonald’s  function  of  order  1/ 
(cf.  Ref.  22),  sometimes  also  called  the  modified  Bessel  function  of  the  third  kind  and  order  1/  (cf.  Ref.  23), 
we  use  the  equality 


for  z  ^  0  and  |  arg  z|  <  tt. 


(5.26) 
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to  replace  the  appearance  of  the  exponential  function  in  favor  of  K^,  subsequently  applying  the  well-known 
formula  asserting  that  the  derivative  of  the  function  z  i— ►  z~''  K„{z)  is  simply  the  function  z  i— ►  —z~‘' K„+i{z). 
In  this  manner,  we  find 


r«- Jr,(y) 


^  °  ^  ^  ^'>g\  (^{t  -  rf  -  dr 

"~16c2  r,(y)  "^4^3  t 

+  _i!_\  e-o.r.(y)/c 

\  (orx(y))^  (arx(y))^  / 

(i^ '  (i 

each  holding  for  all  y  €  r+  and  t  >  d+(x)/e.  Inserting  these  results  into  the  expression  for  J2(x,  t)  and 
performing  some  rearrangements  and  simplifications  (under  which  all  terms  involving  the  exponential  e""''”/' 
disappear),  we  come  to 


1  /  e-°‘  h{'r(t^  -  rl/c^)^  12) 
4c^  I  t  (<2  _  r^/c^)^ 


0 (H''  ■  *)) 

(K»  -  -  i-o‘) 


for  t  >  d+{x)/c.  (5.27) 


There  is  no  difficulty  in  showing  that  the  first  integral  on  the  right  in  (5.27)  approaches  zero  ast-^oo, 
since  we  have  available  for  each  of  the  terms  in  the  integrand  an  estimate  of  the  sort  appearing  in  (5.24) 
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(used  to  analyze  the  second  term  in  the  integral  over  r+  on  the  right  in  (5.23)).  One  should  note  here  that 
Rea  =  Re (7/2),  while  both  y|r+  and  ¥>,n+  are  now  in  C(r+),  and  so  are  bounded. 


Finally,  we  must  deal  with  the  second  and  third  integrals  taken  over  r+  on  the  right  in  (5.27).  To 
evzduate  the  pointwise  limits  of  the  integrands,  t.e.,  of  the  inner  integrals,  we  first  write 


^  =  Uj  Tl  (2 

Ki^a(T^ +  r^(y)/c^)^')  ^  r-y  \  ^ 

V’tJ  c  Jo  (r2  +  r2(y)/c2)‘  V2  / 


dr 


(5.28) 


and 


/2\^at  ,  r*-i’'*(.y)  Ks.(a(i  -  t))  fy/  1  \j 

n,(,,(;x):=(-)  _(r.),„,(y)^  _,J(y,) 

=  b  5;j('-x).n+(y)  / 

\^J  2c^  Jo  (r2  +  ra(y)/c2)  *  ''2  J 


dr 


(5.29) 


for  y  €  r4.  and  t  >  d^{x)/c. 

Now,  for  the  evaluation  of  the  limits  as  <—►  00  of  the  integrals  on  the  right  in  (5.28)  and  (5.29),  we  have 
prepared  Lemma  A.2  of  the  Appendix.  In  the  statement  of  that  Lemma,  we  wish  to  take  a  =  rx(y)/c  (with 
y  chosen  in  r+),  /?  =  *7/2  =  — (Im7/2)  +  i(Re7/2),  and  z  =  a  =(Re7/2)  +  i((Im7/2)  -  w) ;  in  particular, 
we  shall  then  have 

JaiPr)  =  Jo  (*|^)  =  h  for  r  >  0 

(with  Jo  denoting,  of  course,  the  Bessel  function  of  the  first  kind  and  order  zero)  and  Re^  =  Im/?  = 
Rey/2  >  0,  and  consequently  must  use  equality  (A. 19)  for  the  evaluation  of  the  limits  in  question.  With 
these  identifications, 

z^  +  =  -U)(u  +  ij)  =  ]u{u)  +  i7)|e*('’"’^^  ^  0, 

in  which  we  have  denoted  the  principal  argument  of  u)(u  + 17)  by  d.  Since  Re7  >  0  (and  w  >  0),  lies  in 
(0,5r),  whence  0  —  t  lies  in  (— •)r,0),  and  we  compute,  recalling  (1.19), 

\/z^~+^  —  \/ —w{<jj  +  17)  =  |w(w  4-  i7)|e‘'’^^e~”'^^  =  —i\/u){u  +  i-y)  =  — ic/c+. 

Provided  that  Re^  >  1,  (A.  19)  then  gives 

/■C^(a(r2  +  r2(y)/c2)^) 

- T-^/o  (Jr)  rdr 

(r2  +  r2(y)/c2)^ 

By  first  taking  ^  =  3/2  in  (5.30)  and  recalling  (5.26),  it  is  easy  to  check  that 

g«<t+r,(y) 

lim  222(y,<;x)  =  - 7—^  for  each  y  e  r+ ;  (5.31) 

<—00  rx(y) 


I  rx(y)/c  J 


Af  - 1  (rx  (y )  n/zM^/c) 
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by  setting  ^  =  5/2  in  (5.30)  and  using  (5.26)  with  the  recursion  relation  for  to  compute 


one  also  finds 


lim  Iisiy, t\x)  =  r«(y)r*,„  (y)  (  — I - ^ - y  + - ^ - y 

\'x(y)J  \(-i„yr.(y))*  (-■«+r.(y))S 

/gt/c+r,  \ 

=  -  I  — j; — )  ,n+(y)  for  each  y  G  r+. 


p**+’',(y) 


(5.32) 


Now,  to  prove  that  the  pointwise  limits  in  (5.31)  and  (5.32)  imply  (5.20),  we  have  only  to  produce  appropriate 
estimates  for  the  integrands  of  the  second  and  third  integrals  taken  over  r+  in  (5.27).  To  this  end,  we  cite 
the  bound 


(5.33) 


wherein  M((a)  and  l{(a)  are  certain  positive  numbers  depending  upon  only  the  indicated  parameters, 
following  from  an  asymptotic  formula  for  K^{z)  that  is  valid  when  |arga:|  <  ir  -  6,  with  6  fixed  in  (0,5r) 
(cf.  Eq.  (5.11.9)  of  Ref.  22).  With  this,  we  find  that  the  integrands  in  question  are  uniformly  bounded  for 
y  G  r+  and  all  sufficiently  large  t.  For  example,  consider  the  second  term  on  the  right  in  (5.27):  if  y  G  r+ 
and  t  >  to  >  d+(x)/c, 


dr 


<  max 
“  »er+ 


dr  >  . 


By  using  (5.33)  and  recalling  (5.22),  the  second  integral  on  the  right  here  is,  for  a  fixed,  sufficiently  large  tg, 

4  /■<  e-(H.eot)T  g(R.e7/2)(T»-rJ(y)/c*)^  4  1 

~  {2\a-y\)i  Jto  1x2  -  r2fvl/c2'l‘  ~(2|a7|)W,„  r(r2  -  (d+fxl/cl^)  * 


{2\Qy\)  2  Jt,  T  (r2 -r2(y)/c2)‘ 

and,  for  that  <<,,  the  first  integral  on  the  right  is 


It  fy  \\  f  j 

<  max  /olxs)  /  — ^-r - dr, 


inequalities  providing  uniform  bounds  as  claimed.  The  third  term  on  the  right  in  (5.27)  can  be  handled  in  a 
similar  manner,  whence  we  conclude  by  (5.31),  (5.32),  and  the  Dominated-Convergence  Theorem  that  (5.20) 
is  indeed  correct.  This  completes  the  proof  of  Theorem  5.1.  ■ 

Proof  of  Theorem  5.2:  Fix  u^,  G  corresponding  to  this  choice,  set 

Uu/(<)  :=  Uu,e~’"‘  for  each  f  G  R, 
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to  get  an  element  of  C°°(R;  H^)  with  the  properties  A*  ^Uu,(t)  +  =  0,  T“‘’[Ui^(<)]  =  0,  and 

Ua,(t)|r^.  •  n+  =  0  for  each  <  €  R;  obviously,  any  one  of  the  derivatives  of  Uu,,  say,  U/1,,  possesses  the  same 
properties.  Consequently,  the  equality 

f  {U:(s)  .  At,^U(s)  -  U(s)  •  dAg 

J(io 

=  f  {U:(s)|r,  •  T“»[U(s)]  -  U(s)|p„  .  T”‘’[u:(s)]}  dAr„, 
dr, 

holding  for  each  s  G  R  and  resulting  from  two  applications  of  (2.9),  leads  to 

go  [  {  U"(s)  •  U:(s)  +  w=U(s)  •  U:(s) }  dXa 
Jtlo 

=  iwc‘“*  To(s)  ■  dAr,  +  J  Fo(s)  •  dA3|  for  each  s  G  R,  (5.34) 

in  view  of  (5.8)  and  (5.9).  We  note  here  that  the  function  need  be  considered  no  further.  By  using  (5.6) 
and  (5.10),  in  particular,  by  computing  U"(s)  from  the  latter  equality  as  U’'"(s)  —  ui^ue“‘“*  for  each  s  G  R, 
a  first  application  of  (5.34)  yields 


go  I  {U^"(s)+a^^U^(s)}  .  U:(s)dA3 


=  iwe*"* 


{Jr 

{/r 


Uu,  |r.  dAp  + 


r, 

to  •  Uto  Ip  dAp  + 


jf  F:(s)  u,odA3| 

J  fo(s)  •  Uu,  dA3| 


for  each  s  G  R. 


Directly  from  the  latter  equality,  (5.12)  results  by  employing  the  conditions  of  transience  (5.7)  and  (5.11) 
(for  j  =  0  and  2),  the  vanishing  of  u^^jp^  •  n+,  and  the  Cauchy-Schwarz  inequality. 


For  the  proof  of  (5.13),  we  return  to  (5.34).  By  noting  that 

U"  ■  U:  +  •  Ui  =  (U'  -  U:,)'  +  (U  •  U,,)' 

holds  by  virtue  of  the  relation  U"  =  — and  again  using  (5.6),  but  now  taking  into  account  (5.12),  we 
are  led  this  time  from  (5.34)  to 

Qo  f  (U'  •  UL  +  ■  U„)'  (s)  dAg  =  rwe'^  (  /  T^s)  ■  u,olr,  dAp^  +  /  FJ(s)  •  dA3 

dfi,  Ur,  dn. 


for  each  s  G  R. 


Choose  any  <  >  0  and  integrate  this  equality  over  (0,<);  applications  of  the  Cauchy-Schwarz  inequality  and 
accounting  for  the  continuity  of  the  functions  s  t-*  ||U(-'^(s)||j,o  serve  to  verify  that  Fubini’s  Theorem  can 
be  invoked  to  justify  reversing  the  order  of  integration  in  the  resultant  left-hand  side.  By  recalling  that 
U(0)  =  U'(0)  =  0  and  using  (5.10)  to  replace  the  appearance  of  U  and  U',  we  come  to 


2goUJ^  f  uu^dX3  +  gof  {U"'(<)  ■  U^t) -H  ^^3 

dn,  dn, 

=  iu,  T,^(s)  ■  u,,|p,  dAp^  +  Flis)  ■  u„  dA3|  ds, 
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holding  then  for  each  t  >  0.  Appealing  to  (5.11),  this  time  for  j  =  0  and  1,  we  see  that  the  second  term  on 
the  left  here  has  the  limit  0  as  <  — ^  oo.  Thus,  the  limit  of  the  integral  on  the  right  exists,  and  (5.13)  results 
upon  letting  <  oo.  I 


Finally,  supposing  that  the  hypotheses  of  Theorem  5.2  remain  in  force,  we  remark  that  it  is  easy  to 
identify  additional  conditions  on  under  which  u  satisfies  (1.24).  For  example,  suppose  that  U’'(<)  € 
for  each  <  6  R  and  lim(_oo  ||U’^(<)IIh*  =  0.  Then  u  must  lie  in  and  (5.8)  will  give,  with  the  limits  being 
taken  in  the  norm  of  , 


lim  e' 

t— *oo 


A5;,^U^(0  +  eoU^"(t)}  -  =  lim  +  U, 

«r-  If-  t-^OO 


whence  (1.24)  follows.  Similarly,  one  can  formulate  reasonable  conditions  of  this  sort  on  the  functions  U,  <I>, 
To,  Fo,  and  on  the  summands  in  their  forms  assumed  in  §1,  under  which  (1.18),  (1.25),  and  (126)  are 
consequences  of  (1.7),  (1.13),  and  (1.14),  respectively. 
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Appendix 


AUXILIARY  RESULTS 


We  have  collected  here  the  formulations  and  fairly  detailed  proofs  of  two  statements  that  are  used  in 
the  proof  of  Theorem  5.1. 

The  first  result  provides  an  integral  representation  for  sufficiently  regular  and  (for  simplicity)  initially 
quiescent  solutions  of  the  so-called  damped-wave  equation  (110)  or  (5.1)  in  the  exterior  of  a  cylinder  in 
space-time.  The  representation  coincides  with  that  resulting  from  the  use  of  Green’s  Theorem  in  conjunction 
with  a  causal  fundamental  solution  (in  the  distributional  sense)  for  the  hyperbolic  operator  involved  (cf., 
t.g.,  Ref.  A.l),  and  reduces  to  the  classical  Kirchhoff  expression  (cf.,  e.g.,  Ref.  A. 2)  in  the  case  of  the  ordinary 
wave  equation.  A  full  and  precise  statement  of  the  representation  does  not  seem  to  be  readily  available,  and 
so  we  have  chosen  to  outline  here  the  main  points  in  a  derivation;  the  present  proof  does  not  employ  the 
theory  of  distributions  but  is  more  in  the  spirit  of  “elliptic”  reasoning,  essentially  consisting  in  an  extension 
of  the  development  of  the  Kirchhoff  result  as  presented  in  Ref.  A. 2.  Actually,  we  shall  consider  a  more  general 
hyperbolic  operator,  since  but  little  additional  work  is  thereby  required. 


It  is  convenient  to  introduce  some  notation,  to  prepare  for  the  statement  and  proof  of  the  first  Lemma. 
When  X  €  R^,  let  the  corresponding  distance-function  r*  on  R®  be  defined  by 

'•x(y)  :=  |y-x|  foryeR®. 


Suppose  that  is  a  complex  function  defined  on  a  cylinder  F  x  R,  for  some  F  C  R®.  With  $,  c  >  0,  and 
any  chosen  x  €  R^  and  <  €  R,  we  associate  a  function  defined  on  F,  denoted  by  ^y.  F  C  and  termed 
the  (c-)retardation  of '6  with  respect  to  (x,t),  according  to  the  rule 


[^](x,o(y)  •=  (y<^-  y  e  f. 

Recall  our  convention  of  §1:  for  an  appropriate  complex  function  ^  defined  in  0+  x  R,  by  we  mean  the 
function  on  r+  x  R  obtained  by  taking  the  spatial  normal  derivative  at  fixed  values  of  the  fourth  argument. 
That  is,  ^,„^(',0  :=  ^('.0.0+  1’+  <  €  R.  We  shall  write  Ba{x)  :=  {  y  G  R^  |  |y  -  x[  <  a }  for  the 

open  ball  in  R^  of  radius  a  >  0  and  centered  at  x  G  R^.  Finally,  we  use  the  standard  notation  F  for  the 
modified  Bessel  function  of  the  first  kind  and  order  w,  cf.  the  definition  in  Eq.  (5.7.1)  of  Ref.  A. 3. 

Lemma  A.l.  Let  the  open  set  in  R^  be  connected  and  the  complemeni  of  the  closure  of  a  bounded 
connected  Lipschitz  domain;  write  r+  :=  <9f2+,  and  let  n+  denote  the  unit-normal  field  that  is  defined  at  all 
appropriate  points  of  r+  and  “directed  into  P4..”  Suppose  that  c  is  a  real  positive  number,  while  jo  ^nd  71 
are  complex.  Let  $  G  C^(fl+  x  R)  satisfy 


-Ad>  -I-  -I-  ^<^,4  q-  -4«^,44  =  0,  in  n+  X  R, 


(A.l) 


and  vanish  in  x  (— oo,0); 


$(•,<)  =  0  in  124.  for  <  <  0. 


(A.2) 
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Then  ^  has  the  representation  given  by 


<^(x,<)  =  - 


_1_ 

4w 


g-yir,f2c 


+  ^(»‘x).n+/{  ((<  -  rf  -  r)|  dr^ 


I 


for  X  6  !i2+  and  <  £  R,  (A. 3) 


wherein 

^ -=71-470,  (A. 4) 

and  -s/^'  denotes  a  selected  square  root  of  corresponding  to  which  the  function  :  R\  {0}  — ►  C  is  defined 
by 

if  X  <  0  . 

/((^)  ■-  I  l^(,/^>xi/2)/xi  ifx  >  0 

(so  that  fo  is  the  zero-function). 

Proof:  A  preliminary  transformation  serves  to  simplify  the  form  of  the  underlying  partial  differential  equation: 
specifically,  setting 

^(y,  s)  :=  ^(y,  s)e'’'‘’^^  for  y  £  0+ and  s  £  R,  (A. 6) 

we  get  a  function  ^  £  C^{fl+  x  R),  vanishing  in  x  (— oo,0)  and  satisfying,  with  ^  defined  as  in  (A. 4), 


-A'P  -  4- ^4',44  =  0,  in  X  R,  (A. 7) 

as  one  can  readily  check.  In  the  remainder  of  the  proof,  let  x  be  fixed  in  17+;  until  further  notice,  let  t  be  fixed 
and  positive.  Then  [^](x  t)  and  0  are  defined  in  0+  and  lie  in  (7(0+),  while  the  supports  of 

and  ['I',4](x  I)  are  compact  in  0+,  in  fact,  are  contained  in  the  closure  of  the  set  0+  n5ct(x),  since 
and  [’I',4](x  <)(y)  vanish  when  t  -  (l/c)rx(y)  <  0.  Further,  from  the  chain-rule  computations  yielding 


W(x,<)o(y)  =  ['^->](x,<)(>')  -  T''xv(y)f'^-4l(x,«/y) 


for  y  £  0+  \  {x],  j  =  1,2,3, 


and  the  further  expressions  giving  and  ['I',4](x.j)d  \  ^  2,  3,  it  is  apparent  that 

W(x  <)  ['^’4](x  ()  partial  derivatives  at  x,  but  their  restrictions  to  0+  \  {x}  do  lie  in 

C^(0+  \  {x})  and  (7^(0+  \  {x}),  respectively,  with  their  first  partial  derivatives  remaining  bounded  in  that 
set.  A  bit  of  computation,  along  with  an  appeal  to  the  obvious  fact  that  the  retardation  with  respect  to 
(x,  ()  of  the  function  appearing  on  the  left  in  (A. 7)  vanishes  throughout  f7+,  will  verify  that  the  Laplacian 
of  satisfies  (cf.  Ref.  A. 2,  for  the  case  ^  =  0) 


lx.tpjj 


in  a+  \  {x}. 


(A.8) 


For  the  moment,  let  /  denote  a  complex  function  defined  in  R  \  {0},  with  f{s)  =  0  for  s  <  0  and 
the  restriction  /|(0,oo)  £  C^([0,oo));  we  shall  presently  choose  /  judiciously.  With  /,  x,  t,  and  any  ^  £ 
C(ft+  X  R),  we  associate  the  function  {/  ★^}(x,o  ^  ^  according  to 

ft-ir,(y)  /  \\- 

{/*'*'}(x,t)(y)  :=  ^»-*(y)j fory£n+. 
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The  argument  of  /  appearing  here  is 

(t  _  rf  -  ^rl(y)  =  {  (<  -  ^’•«(y))  “’■}{(*  +  ’ 

which  is  positive  when  0  <  r  <  <  — rx(y)/c  and  negative  when  t  —  r^(y)/c  <  r  <  0.  In  particular,  {f*^}(x,t) 
vanishes  outside  the  set  n  Bct(x),  since  f{s)  =  0  for  s  <  0;  actually,  the  values  of  /  for  negative  values  of 
its  argument  will  be  irrelevant  for  us,  since  we  shall  be  interested  exclusively  in  the  cases  when  ^  is  either  9 
or  each  of  which  vanishes  in  x  (— oo,0).  The  regularity  properties  of  {/*’4'}(x,t)  can  be  discerned  a 
bit  more  easily  by  making  a  change  of  variable:  setting,  for  definiteness,  /(O)  :=  /(O'*')  ;=  lim,_o+  /(«),  we 
can  write 


=  (t-  “'■*(y))^  /((^  -  “'■^(y))  (1  “  ((<  +  “  (*  “  *))  ^  (y’(' "  *) 

whenever  y  €  f2+  (by  first  supposing  that  i  —  rx(y)/c  /  0,  and  then  noting  that  the  result  is  true  also  when 
t  -  rx(y)/c  =  0).  From  the  latter  form,  it  is  a  simple  matter  to  check  that  {/  ★  'i'}(x,()  is  in  C(n+),  and, 
moreover,  that  the  functions  {/★^}(x,«)  and  {/*(’*',;  )}(x,0  of  particular  interest  are  at  least  in  C^(n+\{x}) 
and  C^(n+  \  {x}),  respectively,  with  their  first  partial  derivatives  bounded  in  0+  \  {x}.  The  only  difficulty 
at  all  lies  in  examining  the  behavior  of  these  functions  near  x  and  r+  and  on  the  set  0+  n(9Bc<(x).  Consider, 
for  example,  f'’om  the  second  form,  it  is  easily  seen  that  the  first  partial  derivatives  of  {/■*-^}(x,*) 

exist  and  equal  zero  at  any  y  €  12+  with  i  —  rx(y)/c  =  0,  owing  to  the  regularity  hypothesized  for  /  and 
^  and  the  vanishing  of  ®  in  12+  x  (-oo,0).  At  a  point  y  G  f2+,  y  ^  x,  and  contained  in  the  ball  Bct(x), 
the  usual  formula  for  differentiation  of  a  parameter-dependent  integral  can  be  applied  to  compute  the  first 
partial  derivatives  of  {/*^}(x,«)-  Upon  doing  so,  an  inspection  shows  that  the  resultant  expression  vanishes 
outside  Bet(x),  and  so  is  valid  in  all  of  12+  \  {x};  that  is,  we  have 

{/★^}(x,()u(y)  =  -^/(O'^)''x,j(y)['i’](x,o(y) 

/(-ir.(y)  I-  /  ,  \ 

+ 1  |/((<-^)^-^^x(y)j%(y>r) 

-^(»-*)u(y)/'  -  Tf  -  ^’■x(y))  ^(y-T)| 

for  y  6  12+ \  {x}.  (A. 9) 

Directly  from  (A. 9),  it  is  obvious  that  the  first  partial  derivatives  of  {/★^}(x,t)  can  be  extended  continuously 
to  all  points  of  F^.  and  are  bounded  in  12+  \  {x}.  The  second  partial  derivatives  of  {/★^}(x,t)  can  be  shown 
to  be  in  C(12+  \  {x})  (but  not  necessarily  bounded  there)  in  a  similar  manner,  starting  from  (A. 9).  The 
study  of  the  first  derivatives  of  {/ *(^,j)}(x,r)  is  entirely  analogous. 

Now,  by  using  (A. 7)  and  the  vanishing  of  ^  in  12+  x  (— oo,0),  an  elementary,  albeit  rather  tedious, 
computation  reveals  that  the  Laplacian  of  {/  *^}(x,t)  in  12+  \  {x}  satisfies 
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W(x,o(y) 


+2/' 


-  -  ^’'x(y))  f"  ((‘  -  -  ^’’x(y)) 

-  rf  -  ^rl(y)j  -  ((t  -  rf  -  |  ^(y.  r)  dr, 


for  y  G  n+ \  {x}.  (A. 10) 


Evidently,  we  should  select  /|(0,oo)  in  such  a  way  that  this  restriction  lies  in  C^([0,oo))  and  the  bracketed 
factor  in  the  integrand  on  the  right  in  (A. 10)  vanishes  for  0  <  r  <  <  —  rx(y)/c,  for  each  y  G  0+  PI  Bct{x); 
since  (t  —  r)^  —  r^(y)/c^  is  positive  for  all  such  r  and  y,  we  demand  that 

sf'is)  +  2/'(s)  -  ~f(s)  =  0  for  all  s  >  0. 

It  follows  that  we  must  take  /l(0,oo)  to  be  given  by  a  multiple  of  f(,  wherein 

/:  (v^'sV2) 

/^(s)  := - 1 -  for  s  >  0,  (A. 11) 

with  -y/?'  indicating  a  fixed  choice  of  square  root  of  ^  (in  which  case  we  shall  even  have  /|(0,oo)  G 
C“([0,oo))).  For  a  reason  that  will  become  apparent  shortly,  we  choose  the  multiple  so  that  /(O'*")  =  ^/8c; 
since  we  find  /{(O'*")  =  \/?V4  from  (A. 11),  this  requires  that  we  take 

/(«)  := 

Now,  with  /  so  specified,  (A. 10)  gives 


{/  ★  ®}(x,()';;  —  2{/  *  {®)i)}(*,0>i 


in  n+  \  {x}; 


combining  this  with  (A. 8),  we  come  to 


-m 


(x.O’JJ 


{/ ■*■  ^}(x,()i  j;  —  ~  ^  ^  ^  Tx,,- [’t',4](x  t)^  u  —  2{/ ★  (^,j)}(x,t),j ,  inn-|.\{x}.  (A. 12) 


Now  let  the  positive  number  R  be  greater  than  ct  and  so  large  that  the  ball  f?ft(x)  contains  for 

all  sufficiently  small  positive  f,  denote  by  the  bounded  open  set  {f2+  PI  Sr(x)}  \  S((x),  with  boundary 
-  r+  U(95fl(x)U(95f(x).  For  any  such  c,  according  to  the  preparation  laid  down  in  the  discussions  of 
the  regularity  properties  of  [^](x,t).  ['J'.4](x,()’  {/*'*')(x.O’  {/*('*'>;  ))(x,0’  ^PP’y  Green’s  Theorem 

to  each  of  the  terms  in  the  integral  of  the  left-hand  side  of  (A.  12)  over  fi”’'  (recognizing  that  the  Laplaciansof 
both  the  function  l/r^  and  the  function  having  everywhere  the  value  1  vanish  in  fi^'*),  while  simultaneously 
using  the  Divergence  Theorem  to  transform  each  of  the  terms  in  the  integral  of  the  right-hand  side  over  the 
same  open  set.  By  proceeding  in  this  manner  and  using  the  support  properties  of  the  integrands  along  with 
the  condition  t  -  R/c  <  0  to  conclude  that  all  resultant  integrals  taken  over  the  portion  dBf{{x)  of 
must  vanish,  we  get  (with  n+j ,  j  =  1,2,3,  denoting  the  Cartesian  components  of  n+) 


52 


NRL  REPORT  9172 


Ir  {  (tx)  (  Tx  {/*®}(x.«)>n+|  ‘^•^r+ 

~  idB  (x)***’^  {{^}  t^l(x.<)’.»  “  (tv) ~  I  <^^dB,(x) 

~  X"  *  {(^)  ~  <^^3 

^  ii"'  {~c  (;r'‘’'’it'^’4](x,,))  -2{/*('^,>)}(x,o|,>‘^A3 
~  ^  (  (cr  )  ^*’“+^^’4l(x,0  {/ *  (®>i)}(x,t)ii+i  ^  <^-^r+ 

so,  noting  that 

[^in+](x,t)  =  W(x,0’"+  ~’‘*»«+l^»4](x,()lr+>  on  r+, 

we  have 

/b.(x)  {(^)  ^  (ii)  ~ 

■*■  (^)  t^>4](x,«)  +  2{/*($,j)}(x,()  7x!j|  <i^dBt(x) 

~  ~  Ir+{  (^)  ’"+f'^J(x.‘)  ■*■  (^)  ’■*’"+ [^’4l(x,0 

+  2{/*(^,j)}(x,t)n+j  -  {/★^}(x,<)>ii+|  tiAr+.  (A. 13) 


The  expression  on  the  left  in  (A.  13)  is  consequently  independent  of  the  sufficiently  small  positive  e;  its  limit 
for  €  — ►  O"*"  is  easily  proven  to  be  4x$(x,<),  by  noting  that  the  first  partial  derivatives  of  r^,  [^](x()>  ^^d 
{/*^}(x,()  are  bounded  in  Q+  \  {x},  appealing  to  the  continuity  of  [^](x,t).  ['^.4](x  t)>  {/*(’^.j)}(x,t)  in 

n+,  and  observing  that  [^](x  <)(*)  —  Thus,  by  inserting  the  definition  of  {/ * {'®',>)}(x,<)  and  using 

(A. 9),  from  (A. 13)  we  derive  the  equality 


=  - 


8^2 '■*'"+)  W(x,«) 

^••x)  "l^.n+C-.r) 


+  ^(r2),„^/'  ^(<  -r)2  -  't(-,r)|  drj  dAr+. 


(A.14) 


Upon  using  (A. 6)  to  replace  each  appearance  of  'I'  in  (A.14)  by  the  appropriate  expression  involving  <I>,  it  is 
a  simple  mailer  to  verify  that  (A. 3)  results  after  some  rearrangement.  Finally,  we  have  assumed  here  that 
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t  >  0,  but  it  is  obvious  that  (A. 3)  holds  for  t  <  0,  as  well,  since  the  expression  on  the  right  in  that  equality 
vanishes  for  any  such  t.  ■ 

Next,  we  turn  to  the  explicit  evaluation  of  an  integral  of  Sonine-Gegenbauer  type,  for  which  it  is  first 
necessary  to  interject  some  remarks  about  complex  algebra.  With  Log(  )  denoting  the  principal  branch  of 
the  logarithm- function,  given  as  usual  by 

Log  C  :=  In  Id  -f  i  arg  C  for  C  #  0, 

wherein  In  a;  is  the  natural  logarithm  of  the  positive  number  x  and  argC  £  (— 5r,7r]  indicates  the  principal 
argument  of  the  complex  number  C  0.  throughout  C  C“  signifies  the  principal  branch  of  the  general 
power  function: 


<“  :=  exp(aLogC)  =  exp(aln  |d  iaargC)  for  <  6  C,  C  7^  0,  and  a  e  C. 
When  Q  =  1/2,  we  may  write  in  place  of  Then  we  have  the  following  rules  of  calculation: 

=  and  r“  =  i/c“, 

but,  if  a  5^  0, 

(CiC2r=CrC2  and  (d/C2r  =Cr/C2“  iff  -’r<argC,+argC2<^; 
moreover,  in  general  we  have  only 

_  ^a0g2rin(<:,a)/} 

in  which  n(dQ:)  denotes  the  unique  integer  such  that 

arg(e‘‘"’^‘‘^°*''i)  =  Im{aLogC)  +  2n(C,  a)7r. 


Consequently,  if  ^  0,  then 

(C“)^  =  iff  —  T  <  Imaln  Id  +  ReaargC  <  TT. 


For  example,  the  equality  (C“)^  =  holds  if  both  and  a  are  real  or  if  only  a  is  real  and  —  1  <  a  <  1. 
Later,  we  shall  employ  these  facts  without  explicit  mention. 

Now,  with  A'f  denoting  Macdonald’s  function  of  order  (,  Ji,  the  Bessel  function  of  the  first  kind  and 
order  u,  and  z,  a,  and  /?  complex  numbers,  let  us  formally  set 


1(2  I  a,  /?;  d  v 


K({z  +  a^) 
{Vx^  + 


dx. 


(A. 15) 


The  standard  works  provide  an  explicit  evaluation  oiJ{z\a,P;^,iy)  for  2  and  0  real  and  positive,  a  complex 
with  its  principal  argument  arg  a  restricted  by  |arga|  <  7r/2,  and  and  u  complex,  with  Rei/  >  —1;  under 
these  conditions, 

(. - .. 

AV.-,(a^/I^T^),  (A. 16) 

as  shown,  e.g.,  by  Watson  [A. 4,  p.  416,  Eq.  13.47(2)].  However,  our  interest  necessarily  lies  in  evaluating 
T(2  I  a, /3;  ^, //>)  for  certain  nonreal  values  of  2  and  0,  although  only  for  a  real  and  positive.  Watson  (loc.  cit.) 
notes  merely  that  (A. 16)  is  valid  for  complex  2  and  0  “with  certain  limitations,”  while  Erdelyi  ct  al.  [A. 5, 
p,  94,  Eq.  7.14.2(46)]  provide  the  formula  with  no  indication  of  the  restrictions  on  2  and  0.  Evidently,  the 
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evaluation  has  not  been  carried  out  for  values  of  the  parameters  pertinent  to  our  setting,  and  so  we  extend 
the  reasoning  of  Refs.  A. 3  or  A. 4  to  develop  the  required  result  in 

Lemma  A. 2.  Let  (  and  v  be  complex,  with  Rei/  >  —1.  Let  a  be  real  and  nonzero,  and  suppose  that  p  is  a 
nonzero  complex  number  with  arg/3  ^  ir.  Introduce  the  half-plane  C  C  by 

C+  :=  {z€C  I  Rez  >  llmpl). 


(i)  The  integral  appearing  on  the  right  in  (A. 15)  exists  for  each  z  6  Cj ;  the  resultant  function  I(  ■  \a, 
is  analytic  in  Cj  and  is  given  by 

(-v-i  _ 

/f{_v_i(|a|\/z2  for  Re  z  >  (A. 17) 


(ii)  Suppose  it  also  true  that  Re((  —  i/  —  1)  >  0.  Then  the  integred  on  the  right  in  (A. 15)  exists  for  each 
nonzero  z  lying  on  the  boundary  of  Cj,  i.e.,  for  z  0  and  Rez  =  |Im/?|  (which  is  all  of  the  boundary 
of  Cj  unless  ImP  =  0);  moreover,  in  this  case,  the  resultant  values  provide  a  continuous  extension  of 
the  function  I(  •  |  a,  i/)  to  such  boundary  points,  so  that 


f 


-—j===—J„(px)x+'-dx=  hm 


<ec+ 


r  (s/eTW\ 

V  |al  ) 


(-^-1 


for  Rez  =  |Im/?|  and  z  ^  0.  (A. 18) 


Explicitly,  (A.  18)  appears  as 


f 


K({zy/x^  +  a^) 


J„{Px)x‘''^^  dx  =  \ 


(5( 


Vz^  + 


A:f_,_:(|a|v/z2+/?2)  ifz^-hP^^O 


pu2(-' 


-2 


ifz2+/?2  =  0 


for  Im/?  ^  0  and  Rez  =  |Im/?|,  (A. 19) 


and 


f 


K({z\/x^  +  a2) 


(t6)<  I  |a| 


i-u-\ 


ffe_,_i(|a|v//?2-62)  if\b\<p 


(■v/a;2  +  a2)f 


Jv{Px)x'''^^  dx  —  I 


if  161  >  P  and  6^0 


(±»)fy3f-''a2«-2‘'-2 


z^r(^-..-i) 


if  6  =  ±p 


for  P  real  and  positive,  and  z  =  ib  with  b  real  and  nonzero. 
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Before  giving  the  proof,  we  remark  that  we  need  only  the  result  given  in  (A. 19),  for  Im^  /  0,  in  the 
proof  of  Theorem  5.1.  The  form  of  (A. 18)  for  the  case  in  which  /?  is  real  and  positive  and  z  —  ib,  with  b 
real  and  nonzero,  appearing  here  as  (A. 20),  has  been  presented  only  for  the  sake  of  completeness.  If  desired, 
(A. 20)  can  be  recast  into  a  form  involving  the  Hankel  functions  of  the  first  and  second  kinds  rather  than  the 
Macdonald  functions,  by  using  the  familiar  relations  between  these  functions  given  in  Ref.  A. 3,  Eqs.  (5.7.5) 
and  (5.7.6);  in  this  manner,  one  can  derive,  e.g.,  the  formulate  be  found  in  Ref.  A. 5,  p.94,  Eq. 7.14.2(48). 

Proof  :  We  begin  by  studying  the  z-dependence  of  the  convergence  properties  of  the  integral  appearing  in 
(A. 15),  with  the  remaining  parameters  fixed  and  restricted  eis  in  the  hypotheses.  Throughout,  it  is  to  be 
understood  that  x  takes  only  real  values,  while  ci,  Ca,  etc.,  indicate  positive  numbers  that  may  change 
from  estimate  to  estimate,  and  may  depend  upon  the  secondary  parameters  in  a  given  setting,  but  shall  be 
constant  with  respect  to  the  primary  variable(s);  we  may  indicate  the  explicit  dependence  of  one  or  more  of 
these  numbers  upon  certain  of  the  secondary  parameters. 


First,  the  restriction  Rei/  >  —1  ensures  that  the  integrand  in  (A. 15)  is  sufficiently  well-behaved  in  the 
neighborhood  of  the  lower  limit.  For,  given  any  positive  M,  from  the  expansion  of  JuiC)  if  ascending  powers 
of  C  we  have  an  inequality 

U.(C)I  <  ^i(A/)Kr*‘'  for  0  <  Id  <  and  |argd  <  ^r, 
so  that,  if  also  0  <  <62,  there  exists  a  c^iM,  61,62)  such  that 

<C2(A/,5l,i2)x2'^^‘'  +  ' 

for  0  <  X  <  M,  0  <  61  <  I2I  <  ^2,  and  [argzj  <  n,  (A. 21) 
since  it  is  easy  to  check  that  is  bounded  in  any  slit  annulus  {2€C|0<63<|2|<64,  |  argd  <  ir }. 


(Vx2  -1-  a2)f 


Turning  to  the  examination  of  the  integrand  in  (A. 15)  for  large  values  of  x,  we  recall  the  asymptotic 
formulas 


as  Id  00  with  I  argd  £  ’i'  — 


(A.22) 


(Ref.A.3,  Eq.  (5.11.9))  and 


4(0  =  y^{cos  (c  -  -b  (^1^)  }  as  Id  —  00  with  I  argd  <  tt  -  60  (A.23) 

(following  from  Ref.A.3,  Eq.  (5.11.6)),  wherein  5o  is  fixed  in  (0,  tt),  from  which  it  is  clear  that,  given  any 
positive  61,  there  can  be  found  corresponding  positive  numbers  ci,  C2,  and  M'{6i)  such  that 

_ _  g-(Re 

\K^{z\/x^  +  a'^)\  <  Cl - 1 —  for  X  >  A/'(5i)  and  0  <  <  |2|  with  I  arg  2|  <  TT  -  (So  (A. 24) 

(x^  -1-  a2)< 


l4(/?x)|  <  C2 - j —  for  X  >  Af'(6i), 

X3 


(A.25) 


the  latter  inequality  following  since  cos<(  and  sinC  are  always  majorized  in  modulus  by  el''"'-!.  Consequently, 
for  X  and  2  as  in  (A. 24), 


„-(Rei)>/x2+a3  ^Rei/  +  1 

(x2  +  a2)^  (x2  +  a2)^ 


^  2(\/r’  +  a’-r)g-(Re  i-llm  /3|)r 
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which  gives  an  estimate  companion  to  (A. 21), 


+  g^) 


for  X  >  M'{6i)  and  0  <  6i  <  |z|  with  Rez  >  0.  (A.26) 


Now  suppose  that  z  G  C^,  i.e.,  that  Rez  >  |Im/?|:  then,  recalling  that  Rei/  >  —1,  so  2Rei/  +  1  >  — 1, 
inequalities  (A. 21)  and  (A.26)  show  first  that  the  integrand  appearing  in  (A. 15)  is  in  Li{0,oo),  whence 
I{z\a,^;^,i/)  is  defined.  Selecting  any  Zq  €  and  applying  (A. 21)  and  (A.26)  a  second  time,  it  is  clear 
that  the  modulus  of  the  integrand  in  (A. 15)  is  majorized  by  a  nonnegative  function  in  ii(0,  oo),  uniformly 
for  all  z  lying  in  a  closed  disc  contained  in  and  centered  at  zq;  since  the  function  z  i— ►  ^^(z'v/x^  +  a^) 
is  continuous  (even  analytic)  in  for  fixed  x,  Lebesgue’s  Dominated-Convergence  Theorem  implies  that 
I(  •  I  i/)  is  continuous  at  Zq,  and  so  also  in  all  of  C^.  Moreover,  the  same  estimates  imply  that, 

whenever  C  is  a  closed  rectifiable  path  lying  in  C^,  so  that  |z|  and  Rez  —  |Im/?|  are  uniformly  bounded 
below  by  a  positive  number  for  all  z  on  C,  Fubini’s  Theorem  can  be  applied  to  reverse  the  order  of  integration 
in  I{z  \a,p-,^,i>)dz,  whence  Cauchy’s  Theorem  shows  that  the  integral  vanishes,  and  so  J(  •  |  a,  u)  is 
analytic  in  Cj,  by  Morera’s  Theorem.  This  proves  the  first  two  assertions  of  (i). 

In  this  paragraph,  let  Re(^  —  i/)  >  1  hold  (as  well  as  Rei/  >  —1),  and  suppose  that  Rez  = 
Returning  to  (A. 21)  and  (A.26),  it  is  evident  that  the  integrand  in  (A. 15)  is  in  this  case  still  in  ii(0,oo) 
when  Im/?  5^  0  (so  that  z  yi  0),  but  when  Im/?  =  0  we  must  stipulate  that  z  ^  0  to  obtain  the  same  inclusion. 
Thus,  I(-  is  now  defined  by  (A. 15)  at  all  nonzero  boundary  points  of  .  Let  us  show  at  this 

point  that  !( ■  |  a,  /?;^,  u)  is  continuous  in  \  {0},  in  particular,  that  the  same  strategy  employed  to  prove 
the  continuity  in  also  serves  to  verify  that 


J„(Px)x‘''*'^  dx  =  lim 

<ec+\{o) 

for  Rez  =  lIm/?(  and  z  ^  0  when  Re(^  -  i/  -  1)  >  0.  (A. 27) 

For  this,  let  the  nonzero  z  be  fixed,  with  Rez  =  |Im/9|;  obviously,  as  ^  approaches  z  while  remaining  in 
Cp  \  {0},  the  integrand  on  the  right  in  (A. 27)  approaches  that  appearing  on  the  left,  for  each  positive  x. 
Further,  we  again  appeal  to  (A. 21)  and  (A.26),  this  time  to  conclude  that  there  exists  a  nonnegative  element 
of  Li(0,  oo)  that  majorizes  the  modulus  of  the  integrand  on  the  right  in  (A. 27),  uniformly  for  all  ^  lying  in  the 
intersection  of  with  a  closed  disc  centered  at  z  and  of  radius  sufficiently  small  to  exclude  the  origin.  Thus, 
the  Dominated-Convergence  Theorem  can  be  applied  again,  establishing  (A. 27).  Now  the  second  sentence  of 
(ii)  has  been  proven  except  for  the  equality  in  (A. 18),  which  shall  follow  immediately  once  (A. 17)  has  been 
shown  to  be  correct. 

Returning  to  the  general  case,  as  a  preliminary  to  proving  (A. 17)  (as  well  as  (A. 19)  and  (A. 20)),  let  us 
review  the  properties  of  the  function  defined  by  the  expression  appearing  on  the  right  in  that  equality.  Let 
D(3  be  the  subset  of  C  defined  by 

D/3  :=  {  z  G  C  I  arg(z^  +  ^^)  =  n  }  (J  {i/d,  -t/?}  =  {  z  G  C  |  Re(z^  /9^)  <  0  and  Im(z^  +  0^)  =  0  }; 

then  T>0  is  closed,  and  z  »-+  \/z^+^  is  analytic  in  C  \  (recall  that  denotes  the  principal  branch 
of  the  square-root  function,  analytic  in  the  plane  cut  by  removal  of  the  nonpositive  points  on  the  real  axis). 


Jo  (vr'-l-a^)t 


f 


K^{zVx^  +  g^) 
{Vx^  +  a2)« 
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Now,  the  set  Dp  is  found  as  the  collection  of  all  z  =  a:  +  iy  (x,  y  G  R)  lying  in  the  closed  set  described  by 
the  inequality 

x2-i/2<(ini/?)2_(Re^)2 

and  on  the  hyperbola  determined  by  the  condition 

xy  —  -Re,dlm/? 

(which  degenerates  to  coincide  with  the  real  and  imaginary  axes  when  Re/?Im;d  =  0).  Thus,  when  Im/?  =  0, 
so  that  (3  is  real  and  positive,  lies  on  the  imaginary  axis  (just  the  boundary  of  Cj  in  that  case),  comprising 
those  points  ih  with  16|  >  /3;  if  Im/3  0,  is  contained  within  the  closed  strip  {  z  |  |Rezl  <  jlm/?! }  and 
meets  the  boundary  of  that  strip  only  at  the  two  points  ±i0,  while  touching  the  boundary  of  at  a  single 
point  {viz.,  -ij3  if  Im/?  >  0,  i/?  if  Im/?  <  0).  In  either  case,  we  find  that  C  C  \  D^.  Consequently,  the 
restriction  to  of  the  function  z  i— ►  \/ z^  +  j3^  is  analytic  throughout  Cp,  and  it  is  easy  to  see  that  the 
restriction  can  be  extended  by  continuity  to  all  of  C^;  in  fact,  if  Im/d  ^  0,  the  values  of  the  continuous 
extension  are  again  given  simply  by  i/(-)2  +  (3'^,  while  for  ImP  =  0  these  values  are  found  to  be,  recalling 
that  then  0  is  assumed  positive,  and  writing  z  —  ib,  with  6  real, 

'  v/;32  -  52  if  j6|  <  p 

<  (A.28) 

-is/b"^  -  /?2  if  6  <  —j3. 


Next,  since  is  analytic  in  the  open  set  obtained  by  deleting  from  C  the  nonpositive  portion  of  the 

real  axis,  while  \/z^  +  0^  ^  0  and  |  arg  \/z^  +  /?2|  <  jr/2  for  z  €  C  \  D/j,  by  setting 


K{z)  =  K{z\a,0-,^,i.):=  /^^_„_l(lal^/^M^)  forzeC\D;,, 


we  obtain  a  function  analytic  in  C  \  Dp.  Moreover,  from  what  has  been  said,  the  restriction  AC+  ;=  /CjCj 
is  analytic  in  and  can  be  extended  continuously  to  all  points  of  the  boundary  of  Cj  with  the  possible 
exception  of  those  z  on  the  boundary  at  which  z^  +  =  0  (there  being  one  such  point  if  Im/?  /  0,  two 

if  Im/?  =  0).  The  value  of  the  continuous  extension  of  IC+  at  a  point  z  on  the  boundary  with  z^  +  ,0^  0 

is  just  IC{z)  when  Im/?  /  0,  and  in  the  contrary  case  when  Imjd  =  0  can  be  easily  written  down  by  use  of 
the  values  supplied  in  (A.28)  for  the  continuous  extension  of  the  restriction  of  C  \/C^  +  to  Cj .  Now, 
if  Rex  >  0  it  is  easy  to  check,  by  using  the  definition  and  series  expansions  for  (cf.,  e.g.,  Ref.  A. 3,  §5.7) 
that 

lim  eK^{0  =  2^-^Tix), 

larg<|<T 


and  so 


when  Re(^  -  «/ -  1)  >  0  (A. 29) 

iec\D» 


(of  course,  with  F  denoting  the  Gamma  function).  Thus,  when  Re  (^  —  «/  —  1)  >  0  we  conclude  that  K.+  can 
be  extended  by  continuity  to  all  of  Cj,  the  value  of  the  extension  at  the  point  (if  Im0  ^  0)  or  points  (if 
Im/?  =  0)  z  of  the  boundary  at  which  z^  +  /?^  vanishes  being  found  from  (A. 29).  Finally,  the  modifications  in 
these  statements  necessary  to  describe  the  properties  of  the  function  z  (/?*'/z^)A+(z  |  a,  /?;^,  t/),  z  6  C^, 
are  obvious;  in  particular,  the  origin  must  be  excluded  from  the  boundary  points  of  extended-continuity  of 
the  latter  function  if  Im/?  =  0  and  Re^  >  0.  In  fact,  when  Re(^  —  —  1)  >  0,  we  do  have  Re^  >  0  (since 
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Rei/  >  —1),  and  we  find  that  the  values  of  the  continuous  extension  of  this  function  to  the  nonzero  points 
of  the  boundary  of  are  given  by  the  expressions  appearing  on  the  right-hand  sides  of  (A. 19)  and  (A.20), 
in  the  respective  cases  Imy9  ^  0  and  Im/?  =  0. 


According  to  what  has  been  shown,  to  establish  (A.17)  it  is  sufficient  to  derive  the  equality  for,  say,  z 
real  and  greater  than  |Im/3|;  actually,  we  shall  suppose  that  z  lies  in  the  leirger  set  C  specified  by 

Cj+  :=  {C  I  Re«^)  >  (Im^)2  and  ReC  >  0} 

and  prove  that  J(z  |  v)  is  then  given  by  {0‘'fz^)K{z  |  u).  For  this,  we  require  the  evaluation  of 

two  auxiliary  integrals.  First,  from  Ref.  A.6,  p.  146,  Eq.  (29),  we  have 

e-Cir-(c,/(4r))^f-i  ^  2  for  R«Ci  >  0,  R^Cz  >  0,  and  |  €  C.  (A.30) 

Also,  in  Ref.  A.4,  p.394,  Eq.  13.3(4),  we  find 

J  «fi/(Ci7)e~<®’'*r*'+^dr  =  for  Ci  0,  ReCz  >  0,  and  Rei/ > -1.  (A.31) 

Select  and  fix  any  z  €  Cj"*".  Then  we  can  apply  (A.30)  with  Ci  =  1.  Cz  =  +  a^)  (for  x  6  R),  and 

f  ^  to  get  (by  recalling  that  K~(  =  and  noting  that  y/^  =  z,  (z^(x^  -I-  =  zy/x^  a^,  and 

(z^(x^  -I-  a’)/4)”^^*  =  2*/(z^(Vx*  +  o^)t)  hold  because  Rez  >  0) 


I 

I 

r 


and  so 


K({zy/x^  -f  a*) 


2f+r 


I  J"  |  jf”  drj  J4^x)x‘'+i  d 

I  (A.32) 

f 

I 

i  in  which  we  have  used  Fubini’s  Theorem  to  reverse  the  order  of  integration,  a  step  that  we  shall  justify 

I  presently.  We  apply  (A.31)  for  the  evaluation  of  the  inner  integral  in  (A.32)  for  each  positive  r  (a  computation 

that  is  permissible  because  here  Re  (z^)  >  0andRei/>  — 1)  eind  obtain  (2‘"*"^/?‘'/z^*''''^)t‘'+^c“^*’’/'^  ,  whence 

I(z|a,/?;^,i/)  =  2''-«/9''z«-2''-2  dr.  (A.33) 

Jo 

Now  we  observe  that  the  inclusion  z  €  C^'*'  ensures  not  only  that  z  and  z^  have  positive  real  parts  but  also 
that  the  same  is  true  of  z^+)9^  and  (z^  +  /j^)/z*  (and  so  each  number  has  principal  argument  in  (— ir/2,  Jr/2)). 
Indeed,  we  find 

Re(z*  +  P^)  =  Re(z^)  —  (Im/?)*  +  (Re/?)*  >  0, 

and 

=i^{{Re(z*)-KRe;3)*}{Re(z*)-(Im/?)*}-t-{Im(z*)-hRe^Im/?}*}  >0 

(note  that  the  second  of  these  inequalities  provides  an  independent  check  on  the  existence  of  the  integral  in 
(A.33)).  This  shows,  firstly,  that  to  evaluate  the  integral  on  the  right-hand  side  of  (A.33)  it  is  permissible  to 


59 


ALLAN  G.  DALLAS 


apply  (A.30)  a  second  time  by  taking  Ci  =  C2  =  and  |  =  j/  -  ^  +  1,  and,  secondly,  that 

the  following  computations  are  correct:  =  \alz,  and 


UcJ  “V  4  22  +  /?2y  -22f-2./-2(^  |a| 

In  this  manner,  the  equality  in  (A.17)  clearly  results  from  (A.33)  and  (A.30),  and  so  the  proof  of  (A. 17) 
is  complete,  modulo  the  demonstration  that  the  second  equality  in  (A. 32)  is  correct.  Taking  up  this  latter 
question,  we  first  write 


drdx 


Jo  (a^Rc  {z^)(x^  +  u^) ) 


(A.34) 


having  once  more  applied  (A.30)  to  obtain  the  second  equality  here.  To  show  that  the  nonnegative  integrand 
appearing  in  the  final  integral  in  (A.34)  is  in  Li(0,oo),  we  return  to  (A. 21)  and  (A.26),  replacing  there  ^  by 
Re^  and  2  by  ^Re  (2^),  to  conclude  that  for  any  positive  M  there  can  be  found  c^iM)  such  that 


KR.^(N/Re(2^)(z2  +  a2)) 
(VRe(22)(i2+a2))''"^ 


\M^x)\x^‘'+^ 


<C2(M)x^^’'+' 


for  0  <  X  <  A/, 


(A.35) 


and  that  there  exists  a  sufficiently  large  positive  M'  for  which  we  have  an  estimate  of  the  form 


(A/Re(22)(x2  +  a2)) 


Since  2  €  we  have  a/Rc  (2^)  -  llm/?|  >  0,  whence  it  follows  from  (A.35)  and  (A. 36)  that  the  final 

integral  in  (A.34)  is  finite,  and  this  justifies  the  application  of  Fubini’s  Theorem  upon  which  the  second 
equality  in  (A. 32)  is  based.  This  completes  the  proof  of  (A.17),  and  so  also  of  (i). 

Considering  once  more  the  case  in  which  Re  (^  —  1/  —  1)  >  0,  since  the  limiting  relation  (A. 18)  is 
an  immediate  consequence  of  (A.17)  and  (A. 27),  now  the  explicit  forms  (A. 19)  and  (A. 20)  of  (A. 18)  in 
the  respective  cases  Im/?  ^  0  and  Im^  =  0  follow  from  our  discussion  of  the  properties  of  the  function 
2  !-►  {(3’' /z^)IC+{z  \  a,  2  G  C^,  in  which  we  supplied  the  details  necessary  for  checking  that  the 

expressions  on  the  right  in  (A. 19)  and  (A. 20)  give  the  values  of  the  continuous  extension  of  this  function  to 
the  nonzero  points  of  the  boundary  of  .  ■ 
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